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FOREWORD 
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« 


The  effort  described  in  this  report  is  one  part  of  the  theoretical  development 
necessary  to  make  the  last  version  of  the  Naval  Surface  Warfare  Center  (NAVSWC) 
aeroprediction  code,  published  in  1981,  more  applicable  to  present  and  future  weapon 
concepts.  Hie  particular  theoretical  development  included  in  this  report  addresses 
the  changes  in  skin  temperatures  and  aerodynamics  associated  with  real  gas  effects 
that  occur  at  high  Mach  number  flight. 

Appreciation  is  given  to  Francis  Priolo,  who  ran  the  ZEUS  code  to  provide  the 
exact  computations  with  which  to  compare  the  present  new  engineering  methods. 
Appreciation  is  also  expressed  to  Dr.  Roy  Mclnville  for  consultation  on  heat  transfer 
requirements  and  techniques. 

The  work  described  in  this  report  was  supported  through  the  Office  of  Naval 
Research  (Dave  Siegel)  and,  more  specifically,  the  Surface-launched  Weapons 
Technology  Block  Program  mana#*sd  at  NAVSWC  by  Robin  Staton,  .xppreciation  is 
expressed  to  these  individuals  for  the»r  support  in  this  work. 


.Approved  by: 


RICHARD  W.  DORSEY,  Deputy  Head 
Weapons  Systems  Department 
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ABSTRACT 


New  methods  have  been  developed  to  compute  inviscid  surface  pressures  and 
temperatures  for  both  perfect  and  equilibrium  chemically  reacting  flews  on  both 
pointed  and  blunt  bodies  of  revolution.  These  new  methods  include  an  improved 
Shock-expansion  Theorj',  an  improved  Modified  Newtonian  Theory  (MNT),  and  an 
improved  method  for  angle-of-attack  effects.  Comparison  of  these  approximate 
engineering  techniques  to  exact  inviscid  computations  using  a  full  Euler  code  showed 
that  these  new  methods  gave  very  good  agreement  of  surface  temperature  and 
pressures  as  well  as  forces  and  moments.  Incorporation  of  this  new  technology  into 
the  NAVSWC  aeroprediction  code  will  allow  the  code  to  be  used  for  engineering 
estimates  of  inviscid  surface  temperature  at  high  Mach  numbers.  These  approximate 
temperatures  can  then  be  used  as  inputs  to  perform  heat  transfer  analysis. 
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1.  INTRODUCTION 


For  the  past  20  years,  the  Naval  Surface  Warfare  Center  (NAVSWC)  has  been 
involved  in  developing  codes  to  calculate  aerodynamics  on  tactical  weapons.  These 
codes  have  attempted  to  meet  the  changing  needs  of  the  Tacticai  Weapons 
Community  and  keep  pace  with  aerodynamic  requirements.  A  recent  efforti  was 
undertaken  to  look  at  where  we  have  been,  where  we  are,  and  where  we  need  to  go  in 
the  future  with  respect  to  aerodynamic  codes.  This  effort  identified  three  needs:  (1) 
an  improved  aerodynamic  prediction  code  that  allows  Mach  numbers  to  20,  has 
improved  low  aspect  ratio  wii^  lifting  capabilities,  and  allows  improved  base  drag 
prediction;  (2)  an  improved  full  Euler  code  that  has  real  gas  effects  included, 
improved  low  supersonic  Mach  number  ability,  and  an  integral  boundary  layer  code 
to  compute  viscous  effects;  (3)  a  Navier  Stokes  code  to  provide  aeroheating  and  high- 
angle-of-attack  aerodynamics  on  a  limited  number  of  high-value  projects.  This  report 
is  concerned  with  partially  meeting  the  number  1  need  above;  that  is,  including  real 
gas  effects  into  the  computations  so  that  higher  Mach  numbers  can  be  considered. 

The  latest  version  of  the  aeroprediction  code2'4  calculates  aerodynamics  up  to 
Mach  8.  The  main  reason  for  limiting  the  code  to  this  Mach  number  is  that  real  gas 
effects  start  becoming  important  around  Mach  6,  and  at  Mach  8  can  still  be  neglected 
for  most  applications.  However,  as  Mach  number  increases  substantially  above  Mach 
6,  the  requirement  to  have  real  gas  effects  included  in  the  aeroprediction  code 
increases.  Over  the  past  5  years,  tactical  weapon  concepts  have  been  investigated 
that  fly  at  Mach  numbers  up  to  20.  For  this  reason,  it  is  believed  that  real  gas  effects 
need  to  be  included  in  the  aeroprediction  code. 

The  theoretical  methodology  iised  to  predict  aerodynamics  above  about  Mach  3 
is  Second-order  Shock-espansion  Theory  (SOSET)  combined  with  Modified 
Newtonian  Theory  (MNT).  A  brief  search  of  the  literature  revealed  no  present 
second-order  accurate  theoretical  techniques  for  including  real  gas  effects  luLo 
aerodynamic  computations.  The  only  approaches  available  were  either  incorporation 
of  real  gas  effects  into  the  full  inviscid  Euler  equations  of  motion  or  first-order  local 
slope  techniques.  The  first  approach  is  beyond  the  computational  complexity  desired 
for  an  approximate  engineering  code,  whereas  the  latter  approach  does  not  yield 
accurate  axial  force  calculations.  Hence,  extending  SOSET  to  include  real  gas  effects 
will  meet  an  immediate  need  for  the  aeroprediction  code  and  will  also  fill  a  gap  in  the 
external  literature. 
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2.  ANALYSIS 


This  section  is  written  at  a  level  that  covers  the  high-speed  fluid  flow 
phenomena  discussed  in  the  Introduction  from  a  first  principles  viewpoint.  The 
intent  is  to  acquaint  the  reader  who  has  no  prior  background  in  Real  Gas  effects  with 
the  importance  and  the  physics  and  chemistry  of  these  effects  as  well  as  how  one 
mathematically  goes  about  solving  for  the  flow  field. 

2.1  PERFECT  AND  REAL  GAS  BACKGROUND 

2.1.1  Definition 

Usually  when  one  thinks  of  a  definition,  one  thinks  in  terms  of  a  clear  and 
concise  answer.  Such  is  not  necessarily  the  case  with  real  and  perfect  or  ideal  gases. 
This  is  because  there  are  several  definitions  of  perfect  gases  and  real  gases 
depending  on  which  reference  one  uses^tQ  and  whether  one  uses  the  chemist’s  or 
aerodynamicist’s  definition.  To  define  a  real  gas,  we  must  first  define  a  perfect  gas. 
Once  we  have  defined  a  perfect  gas,  a  real  gas  will  be  defined  to  cover  all  ofeer  cases. 

Some  references  define  a  perfect  gas  as  one  that  follows  the  equation  of  state 
identically 


p  =  pRT  (1) 

and  has  constant  specific  heats 
C  C  =  cona 

P.  '• 

(NOTE:  Here  a  bar  is  placed  over  the  Cp,  Cv  to  distinguish  the  specific 
heats  from  pressure  coefficient  Cp,  used  throughout  the  report.) 

Others  break  down  the  definition  of  pierfect  gas  into  a  calorically-perfect  gas,  which  is 
defined  by  Equations  (1)  and  (2),  and  a  thermally  perfect  gas  where  Equation  (1)  is 
still  satisfied  but  the  specific  heats  are  functions  of  temperature.  That  is 


C  =  fiT» 

p 

(3) 

C  =ftT'- 

V 

Finally,  one  author^*  allowed  a  perfect  gas  to  include  cases  where  the  specific  heats 
were  functions  of  pressure  and  temperature,  i.e. 


C 

p 


f'Xp) 


C  =f(T,  D> 

V 


(4) 
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and  Equat  ion  (1)  was  not  required  to  be  satisfied.  The  only  requirement  for  a  perfect 
gas  v/as  that  the  intermolecular  attractive  or  repulsi'-’e  forces  between  atoms  or 
molecules  be  neglected  (this  is  the  original  definition  by  the  chemist). 

For  purc.jses  of  this  report,  the  most  strict  definition  for  a  perfect  gas  is  taken 
(which  is  sometimes  referred  to  as  a  caioricaliy -perfect  gas).  That  is.  Equations  (1) 
and  (2)  must  be  satisfied.  A  real  gas  is,  therefore,  one  where  either  inter- 
molecular  forces  cannot  be  neglected  (low  temperature  and  high  pressure)  or 
one  w'here  the  temperatare  is  high  enough  to  excite  the  internal  energy  modes 
of  vibration,  electronic  excitation,  or  cause  chemical  reactions  and 
dissociation  of  molecules  or  ionizations  of  atoms,  in  some  cases,  Elquation  (3)  is 
valid  (computation  of  additional  internal  energy  from  vibration)  whereas,  in  other 
cases,  chemical  reactions  are  ocemring  or  the  intermolecular  forces  are  not  small  and 
Equation  (4)  applies  and  Equation  (i)  is  no  longer  valid. 

in  summary,  a  perfect  gas  is  defined  here  as  one  that  obeys  Equations  (i)  and  (2) 
exactly,  and  a  real  gas  is  defined  as  one  that  violates  either  Equation  (I)  or  (2)  or 
both.  This  definition  is  consistent  with  most  of  the  aerodynamics  literature. 

2.1.2  Importance  of  Real  Gas  Effects 

One  may  ask:  Why  be  concerned  about  real  gas  effects  for  tactical  weapons 
(since  they  have  been  neglected  up  to  this  point  in  time,  to  the  author’s  knowledge,  for 
aerodynamic  compulations}?  A  good  illustration  for  the  importance  is  given  by 
Anderson.^5  He  shov/ed  the  results  of  static  temperature  behind  a  normal  shock 
where  the  uerfect  gas  assumption  was  made  and  then  these  results  were  compared  to 
a  real  gas.  For  a  perfect  gas,  the  static  temperature  ratio  across  a  normal  shock  is 


»> 

(>•  r  ir 


I. 


(5) 


where  T?  is  temperature  beliind  the  shock  and  subscript  <»>  represents  conditions  in 
the  freestream  ahead  o*'  the  shock.  Figure  1  plots  these  results  as  a  functicn  of  Mach 
number  for  an  altitude  of  170,000  ft.  At  this  altitude,  the  speed  of  sound  is  approx¬ 
imately  1100  Il/sec  and  the  freestream  air  temperature  is  approximately  283‘’K.  The 
normal  shock  would  occur  in  the  vicinity  immediately  ahead  of  the  blunted  portion  of 
a  seeker  or  the  missile  nose.  Note  that  the  temperatures  of  interest  to  tactical 
weapons  aerodynamicists  can  be  very  high  for  high  Mach  number  conditions 
as  ;uming  a  perfec*  gas.  Also  shewn  on  the  figure  are  the  real  gas  results  taken  from 
Rc  'brcnce  8.  Note,  in  particular,  the  plot  of  Tr/Tp,  the  ratio  of  ♦he  real  gas  to  perfect 
gas  temperature.  For  Mach  numbers  of  6  or  less,  this  ratio  is  unity  or  near  unity. 
This  is  the  rea.scn  that  aerodynamic  computations  below  M,  =  6  could  neglect  real 
gas  eiTects  with  little  error.  However,  as  one  goes  above  =  6,  the  error  in 
temperature  using  the  real  gas  assumption  can  be  quite  large.  This  is  of  particular 
importance  to  materials  and  structures  engineers  who  must  design  the  system  to 
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withstand  these  temperatures.  Also  shown  in  Figure  1  is  the  melting  point  of  typical 
structural  materials  used  in  presentrday  missile  design.  The  actual-use  temperature 
is  less  than  the  melting-point  temperature.  For  missiles  that  fly  at  any  appreciable 
time  above  the  maximum-use  temperature  of  a  given  material,  some  form  of  active 
cooling  or  insulation  would  be  required.  This  means  additional  dead  weight  and, 
hence,  less  performance  for  the  missile,  it  is  therefore  obvious  that  a  reasonably 
accurate  estimate  of  temperature  is  essential  for  the  design  of  the  seeker  and  the 
structure  of  the  weapon.  The  reasons  for  this  drastic  effect  on  the  temperature  will  be 
discussed  in  the  nest  section  on  hindamentai  phenomena. 


STATIC 

TEMPERATURE 
BEHIND 
NORMAL 
SHOCK  T2 
{kx10-3) 


FIGURE  1 .  TEMPERATURE  BEHIND  A  NORMAL  SHOCK  .43  A  FUNCTION  OF 
FREESTRE.4MM.4CK  NUMBER  tH  =  iTO.kft; 


Before  we  finish  the  discussion  on  the  importance  of  real  gases,  it  should  be  noted 
that  pressure  is  also  afiected  but  to  a  much  les-s  extent  than  temperature.  This  leads 
one  to  tne  conclusion  that  if  temperature  is  not  of  concern  and  aerodynamic  forces  and 
moments  are  the  major  interest,  the  assumption  of  a  perfect  gas  can  be  made  with 
reasonable  validity  for  Mach  numbers  much  higher  than  6.  The  reason  for  this  is 
that  pressure  is  mainly  governed  by  the  kinetic  energy  of  the  flow  and  not  the 
potential  energ>'.  On  the  other  hand,  the  temperature  is  a  function  of  both  and, 
therefore,  varies  considerably  due  to  the  real  gas  effects.  In  fact,  many  of  the 
hypersonic  flow  assumptions  are  based  on  the  assumption  of  a  thin  shock  layer  and 
>  >  1 ,  and  reasonable  results  for  aerodynamic  forces  can  be  obtained. 
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2.1.3  Physical  f  lienomana  That  Cause  Real  Gas  Effects 

A  gas.  like  any  other  subjrtance,  is  composed  of  atomic  structures  either  in 
singular  form  (monatomic)  or  joined  together  (diatomic  or  polyatomic).  Air  at  room 
temperature  is  composed  primarily  of  nitrogen  (N2)  and  oxygen  (O2)  with  minute 
traces  of  caibon  dioxide  (COzK  argon  (Ar),  and  nitrogen  oxide  (NO).  The  total 
amount  of  these  trace  elements  is  less  than  2  percent  of  the  mass  of  a  given  volume  of 
air,  wheieas  the  N2  is  about  78  percent  and  O2  about  20  percent.  Hence,  for 
engineermg  purposes,  air  at  low  to  moderate  temperatures  (T  <  2500“K)  can  be 
assumed  to  be  composed  of  N2  and  O2- 

These  diatomic  molecules  are  in  continuous  motion  and  can  translate  in  three 
directions  (x,  y,  2)  and  can  rotate  in  two  directions  about  axes  A  and  B  as  shown  in 
Figure  2.  (Note  that  rotational  energy  about  an  axis  going  through  both  N  atoms  can 
be  neglected  in  comparison  to  that  about  the  transverse  axis.) 

B  A 


TRANSLATION  VELOCITY  =  v'Vi2  +  V-2  +  V,2 
ROTATION  VELOCITY  =  ua.  ws 


FIGURE  2.  .\'2  .MOLECULE  UNDERGOLN’G  TR.ANSL.4T10S*  .4.\*D  R0T.4T10.N 


So  as  air  siaj's  beiow  about  600“K.  the  only  energy  it  contains  is  the  rotation 
and  translation  of  the  diatomic  N2  and  O2-  This  is  what  allows  us  ixi  make  the 
assumptions  of  a  perfect  gas  that  must  follow  Equations  (1)  and  (2).  However,  as  the 
air  is  heated  further,  an  additional  energy  mode  enters  into  the  physical  model  of  the 
N2,  O2  molecules.  This  mode  is  vibration.  That  is,  referring  to  Figure  2,  one  could 
imagine  a  spring-mass  system  with  a  spring  that  allows  additional  internal  energy  to 
be  generated  between  the  nitrogen  or  oxygen  atoms.  At  temperatures  above  2500‘‘K, 
the  molecule  has  acquired  about  90  percent  of  the  internal  energy  possible  from 
vibration  of  the  molecules. 

As  the  air  is  heated  above  2000‘'K,  Uie  diatomic  molecules  of  Nv,  O2  begin  to 
dissociate  and  chemically  react.  Air  then  begins  to  contain  significant  amounts  of  O2, 
0,  N2.  N,  and  NO  as  the  temperature  is  increased.  Finally,  as  the  air  is  healed  above 
8000‘’K,  significant  ionization  occurs  tan  electron  leaving  one  atom  and  being  picked 
up  by  another  atom),  which  gives  rise  to  an  additional  mode  of  internal  energy  (the 
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gain  or  loss  of  electrons)  as  well  as  additional  chemical  reactions.  An  example  of  the 
species  composition  of  air  as  a  fl  notion  of  temperature  is  shown  in  Figure  3,  which  is 
taken  from  Reference  1 1  Figure  3  shows  how  the  O2,  N2  go  to  N,  0,  NO,  N  ,  O  -f- , 
and  e  (free  electrons)  as  temperature  is  increased. 


T.  1000  degK 

FIGURE  3.  EQUILIBRIUM  COMPOSITION  OF  .\IR  .AT  DENSm*  OF  10-2  aun 

•FRO.M  REFERENCE  Hi 


The  fundamental  physical  phenomenon  that  occurs  as  air  is  heated  (or  as  Mach 
number  of  a  vehicle  is  increased)  is  that  the  additional  energy  modes  of  vibration  and 
ionization  are  added  and  the  air  molecules  break  down  and  chemically  react  to  form 
other  species.  It  is  this  transfer  of  energy  to  vibration  and  the  dissociation  process 
that  accounts  for  the  fact  that  the  temperature  is  much  lower  for  a  real  gas  than  for  a 
perfect  gas  at  high  Mach  numbers.  That  is,  the  translational  and  rotational  energy 
that  accounts  for  the  high  temperatures  in  a  perfect  gas  is  shared  by  other  forms  of 
energy  and  the  chemical  reactions  in  a  real  gas.  It  is  this  sharing  of  energy  across 
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energy  modes  and  the  conservation  of  energy  that  require  the  temperature  for  a  real 
gas  to  be  lower  than  that  for  a  perfect  gas.  How  one  treats  the  computation  of  real  gas 
effects  from  a  computational  standpoint  is  discussed  later  in  the  Real  Gas 
Computational  Process  section. 


2.2  PERFECT  GAS  COMPUTATIONAL  PROCEDURE 

Before  getting  into  the  mathematics  of  real  gas  computations,  it  is  instructive  to 
briefly  discuss  perfect  gas  computations  on  a  missile-like  body  using  a  combination  of 
Newtonian  Theory  and  SOSET.  Here  it  is  assumed  that  the  missile  has  a  blunt  nose 
where  Modified  Newtonian  Theory  (MNT)  would  be  used  for  flov/  computations  and 
SOSET  would  be  used  past  some  match  point  (between  the  pressures  predicted  by  the 
two  theories).  If  the  nose  is  sharp  or  pointed,  SOSET  can  be  used  directly. 

2.2.1  Modified  Newtonian  Theory 

Newtonian  Impact  Theory  assumes  that,  in  the  limit  of  high  Mach  number,  the 
shock  lies  on  the  body.  This  means  that  the  disturbed  fiow  field  lies  in  an  infinitely- 
thin  layer  between  the  shock  and  body.  Applying  the  laws  of  conservation  of  mass 
and  momentum  across  the  shock  yields  the  result  that  density  behind  the  shock 
approaches  infinite  values  and  the  ratio  of  specific  heats  approaches  unity.  The 
pressure  coefficient  on  the  surface  becomes  1 2 

C  =  2  sinti  '6) 

p 

where  6eq  is  the  angle  between  the  velocity  vector  and  a  tangent  to  the  body  at  the 
point  in  question  (see  Appendix  A). 

Leesl2  noted  that  a  much  more  accurate  prediction  of  pressure  on  the  blunt- 
nosed  body  could  be  obtained  by  replacing  the  constant  “2”  in  Equation  (6)  with  the 
stagnation  pressure  coefficient  Cp^  .  Cp^  can  be  found  from 


where  p-i^T,  ^md  pOo^  a*"-  given  byl3 

p,  2yM;-(y-15 
p.,,  1  +  i 


(7) 


(8) 


p 


y  -  1 
2 


r 


r-J 


(9) 


t 
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Also 


(y  -  l)M^  -f-  2 

M^= — —  (10) 

2yM;-(y-l) 

SO  that  Equation  (9)  can  be  related  to  and  y  only  by  substitution  of  Equation  (10) 
into  (9)  to  obtain 


-y 

4yM^  -  2(y  -  1)  ly-i 
Po  (y  +  ! 


(ll) 


Equations  (8)  and  (10)  hold  only  for  a  perfect  gas.  They  are  the  relationships  taken 
from  Reference  13  for  flow  across  a  normal  shock  wave.  Combining  Equations  (11) 
and  (8)  into  Equation  (6),  we  obtain 


2  \ 

{y+  DM^I 

y-l 

yri  lr-1  A 

vM;1 

2  ! 

2yM^-(y-l>i 

(12) 


and 


C„  =  Cp  sin^S 
P  Pq  «5 


(13) 


Equation  (13)  allows  the  calculation  of  the  pressure  coefficient  ail  along  the  blunt 
surface  of  a  missile  nose  or  wing  leading  edge  for  a  perfect  gas  where  Cp^  is  given  by 
Equation  (12)  and  sin  6eq  from  Appendix  A, 


We  also  must  derive  an  expression  for  static  temperatiu'e  along  the  surface.  To 
do  this  recall  that,  for  an  adiabatic  flow,  total  temperature  is  constant.  This  means 
that 


=  T.(l 


(14) 


Also  recall  that  the  total  pressure  and  density  relationships  are  constant  only  for 
isentropic  flow.  This  means  that  these  relationships  do  not  apply  across  a  shock  but 
they  are  constant  along  a  streamline  for  a  perfect  gas,  frozen  flow,  or  equilibrium 
flow.  Since  the  stagnation  streamline  wets  the  body,  we  can  apply  these  relation¬ 
ships  along  the  body  tr  compute  local  Mach  number  and  then  apply  Equation  (14)  to 
compute  local  temperature.  The  total  pressure  is  given  by  Elquation  {9}  and  the  total 
densit3'  is 


i 


(15) 
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Now  from  Eqiiation  (13) 


Also 


(16) 


Po,  =  Pl 


at  the  stagnation  point.  Thus 
1  .2 

p_  =  —  p  V  C  T  p 
‘'o  2  *  *  Pn  “ 


(17) 


Elquation  (17)  gives  as  a  fiinction  of  freestream  conditions  only  so  it  is  a  known 
number.  Since  Equation  (16)  is  also  known  all  along  the  surface,  then  from  Equation 
(9)  local  Mach  number  can  be  computed  by 


1 

'v=!(rT)l(^F-‘'' 


(18) 


Then  from  Equation  (14) 


-1 


(19) 


Also  from  the  equation  of  state  for  a  perfect  gas 


RT, 


(20) 


Other  properties  such  as  internal  energy,  enthalpy,  and  entropy  could  also  be 
computed,  if  desired,  at  each  point.  That  is 


=  C  T  =f-^  }t 
V  y  -  1  /  *■ 


(21) 
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The  change  in  entropy  across  the  shock  is 


AS  (\\ 

^  =  -(Y-l)ln  — 

C  ^Po/ 


(22) 

(2S) 


(24) 


Since  the  flow  is  isentropic  along  a  streamline  (which  can  include  the  body  siurface), 
AS  along  the  body  is  zero. 


2.2.2  First-order  Shock-expansion  Theory 

First-order  Expansion  Theory  was  first  proposed  by  Eggers  et  al.l4-i7  for  bodies 
of  revolution  flying  at  high  supersonic  speeds.  Basically,  the  Shock-expaiision 
TTieory  computes  the  flow  parameters  at  the  leading  edge  of  a  two-dimensional  (2-D) 
surface  with  the  oblique  shock  wave  relations  and  with  the  solution  for  a  cone  at  the 
tip  of  a  three-dimensional  {3-D)  body.  Standard  Prandtl-Meyer  Expansion  (PME)  is 
then  applied  along  the  surface  behind  the  leading  edge  or  tip  solution  to  get  the 
complete  pressure  distribution  over  the  body  surface.  Referring  to  Figure  4,  this 
theory  inherently  assumes  that  the  expansion  waves  created  by  the  change  in 
curvature  around  the  body  are  entirely  absorbed  by  the  shock  and  do  not  reflect  back 
to  the  body  surface.  Since  the  theory  assumes  constant  press\n-e  along  one  of  the 
conical  tangent  elements  of  the  surface,  fairly  slender  siu-faces  must  be  assumed  or 
many  points  along  the  surface  assumed  to  obtain  a  fairly  accurate  pressure 
distribution.  Another  way  of  stating  this  is  to  minimize  the  strength  of  the 
disturbance  created  by  Mach  waves  emanating  from  the  expansion  comer  and 
intersecting  the  shock,  the  degree  of  turn  should  be  small. 


For  the  2-D  starting  solution,  the  equations  for  a  perfect  gas  across  an  oblique 
shock  areis 


o  o 

p,  2  yM‘ sin'o  -  ly  —  1) 

p,,  r  + 1 


(25) 


o  o 

p,  ^  -r  I'M'an'o 
P=  (y  —  1 )  an*  o  2 


(26) 
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FIGURE  4.  APPROXLMATION  OF  TRUE  BODY  BY  ONE  COMPOSED  OF 
STRAIGHT  UNE  SEGMENTS  TANGENT  TO  SURFACE 
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Here  the  weak  shock  solution  is  taken  as  physically  more  correct  (intermediate  value 
of  sin2  o).  For  a  wedge  flow,  conditions  behind  the  shock  are  constant  so  Equations 
(25)  through  (28)  determine  the  conditions  at  the  surface  as  well.  Total  conditions 
can  be  determined  by  Equations  (9),  (14),  and  (15)  and  energy,  enthalpy,  and  entropy 
change  by  Equations  (21)  through  (24).  The  2-D  starting  solution  would  be  used  for 
the  leading  edge  portion  of  the  lifting  surface  on  missile  configurations. 

For  the  3-D  starting  solution,  the  tip  is  approximated  by  a  cone  w'here  the  initial 
cone  angle  is  defined  by  a  tangent  to  the  surface  at  the  tip  as  shown  in  Figure  4. 
Properties  immediately  behind  the  shock  are  also  calculated  by  the  oblique  shock 
relations  given  by  Equations  (25)  through  (28).  However,  the  flow  between  the  shock 
and  body  is  not  constant  throughout  the  flow  field,  as  in  the  case  with  the  v/edge 
solution,  but  is  constant  along  rays  emanating  from  the  cone  tip.  This  requires 
solution  of  a  differential  equation  to  define  the  shock  wave  angle  and  properties  at 
the  body  surface.  This  was  done  by  Taylor  and  McCall. Without  going  through  the 
derivation  of  the  differential  equation,  it  is  repeated  here  for  reference. 


where 


V‘  =  2  H  =  Const 

Refer  to  Figure  5  for  the  nomenclature. 


FIGURE  5.  NOMENXL.4TURE  FOR  C0N1C.4L  SOLUTION 
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To  solve  this  eqiiation,  we  assume  a  shock  angle  o.  With  this,  the  conditions  behind 
the  shock  can  be  computed  for  a  given  freestream  Mach  number  by  Equations  (25) 
through  (28),  and  the  initial  velocity  components  Vr  and  Vq  behind  the  shock  can  be 
computed.  Equation  (30)  is  integrated  until  Va  =  0.  If  the  angle  ©c  is  not  correct,  a 
new  value  of  shock  angle  is  guessed  and  an  iteration  occurs  until  the  correct  value  of 
©c  is  found  for  Ve  =  0. 

The  Mach  number  at  the  surface  is  then  related  to  the  total  enthalpy  as  follows: 


or 


.r2 


=  1  + 


1 


(y  -  DM" 

Solving  for  M^,  there  is  obtained 


m2  = 


(V/V 


i-(v/v  r 

max 


1 


(y-  D 


(31) 


(32) 


The  temperature  at  the  surface  can  then  be  determined  from  Equation  (14)  since  Tq 
is  constant  throughout  the  flow.  Knowing  the  total  pressure  and  densify  behind  the 
shock  (where  the  fiow  field  is  isentropic)  also  allows  one  to  then  compute  the  local  p 
and  p  throughout  the  flow  and,  in  particular,  at  the  surface  from  Equations  (11)  and 
(15). 


Once  the  starting  solution  and  the  various  fiow  properties  are  known  at  the 
surface  of  a  2-D  or  3  -D  body,  the  PME  Theory  is  applied  for  points  downstream. 

For  a  perfect  gas,  one  can  write ' 


v{M)  = 


/m 

y-  I 


tan 


Y-  1 
y+  1 


(M“  —  I )  —  tan 


-1 


(ad) 


If  A©  is  the  change  in  the  local  surface  slope  in  going  firom  one  tangent  segment  of  the 
body  or  airfoil  to  another  tangent  surface  ^see  Figure  4),  then 


V21M2)  =  A0  +  v.(Mj) 


(34) 
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This  means  that  one  calculates  the  Prandt-Meyer  angle  for  the  first  surface  with 
Equation  {33),  using  from  Equation  (32),  then  must  solve  Equation  (33) 
numerically  for  Mg  based  on  some  value  of  Vg  from  Equation  (34).  Knowing  Mg,  we 
can  compute  Tg,  Pg,  Pg  from  Equations  (11),  (14),  and  (15).  Thus  for  «  =  0,  sharp¬ 
nosed  bodies,  or  airfoils  in  a  perfect  gas,  we  now  have  expressions  for  pressure  and 
temperature  all  along  the  surface  of  the  body  or  airfoil. 


2.2.3  Second-order  Shock-expansion  Theory  (SOSET) 

Syvertsoa  (et  al.)  extended  the  generalized  Shock-expansion  Theory  on  pointed 
bodies  and  sharp  airfoils  to  what  he  called  a  second-order  theory.  19  He  defined  the 
pressure  along  a  conical  frustum  by 

p  ~  p(.  -  (P(3  -  Pg)^"'’  (35) 

instead  of  a  constant  on  each  segment  as  was  the  case  in  the  generalized  theory.  Here 
Pj,  is  the  pressure  on  a  cone  with  the  given  cone  half  angle  equal  to  the  slope  of  the 
conical  segm  mt  with  respect  to  the  axis  of  symmetry.  Pg  is  the  pressure  just  aft  of  a 
conical  segment  (see  Figure  6),  which  is  calculated  from  the  PME,  Equations  (33)  and 
(34). 


FIGURE  6.  FLOW  ABOUT  A  FRUSTU.M  ELEMENT 


Also 


n  = 


(36) 


Thus  examining  p  from  the  equation,  it  can  be  seen,  for  example,  on  the  frustrum 
element  2-3  that  the  pressure  varies  from  the  pressure  of  the  generalized  theory  at 
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point  2  to  that  of  a  cone  of  angle  6^  and  Mach  number  M2  as  s  gets  large.  Reference  19 
approximated  the  pressure  gradient  as 


where 


^1.2  = 


ypMj, 


=  — 
^  M, 
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y  - 1  2 

1  + 


y+l 


y-rl 

ay-n 


(37) 


Finally,  for  negative  angles  such  as  would  occur  on  a  boattailed  configuration,  p^,  was 
replaced  by  p„.  No  discussion  was  given  for  blunt  bodies.  It  should  be  noted  that  if  q 
of  Equation  (36)  becomes  negative,  the  SOSET  reverts  to  the  generalized  or  First- 
order  Shock-expansion  Theory  of  section  2.2.2.  This  is  because  (35)  will  not  give  the 
correct  asymptotic  cone  solution  for  negative  values  of  q. 


Jackson  et  al.20  combined  SOSET  with  MNT  to  treat  blunt-nosed  configurations 
with  or  without  flares.  They,  like  the  authors  of  Reference  18,  assumed  that  the 
lifting  properties  could  be  predicted  by  assuming  that  the  original  body  is  made  up  of 
several  equivalent  bodies  of  revolution  represented  by  the  various  meridians  (see 
Figure  7).  They  assumed  the  match  point  between  the  MNT  and  second-order  shock 
pressure  prediction  to  be  the  angle  that  corresponds  to  shock  detachment  on  a  wedge 
with  the  given  freestream  Mach  number. 

DeJarnette  et  al.2l-23  made  significant  improvements  to  the  work  of  Jackson  et 
al.20  and  Syvertson.i9  These  new  improvements  included  the  following; 

1.  A  new  empirical  equation  for  calculating  pressures  on  blunt  noses. 

2.  A  new  matching  point  to  match  the  blunt-nose  pressures  with  those  of 
SOSET. 

3.  An  exact  (as  opposed  to  an  approximate)  expression  for  the  pressure 
gradient  downstream  of  a  corner. 

4.  A  new  expression  for  pointed-cone  pressures  at  angle  of  attack  which 
improves  the  initial  pressure  prediction  over  that  of  tangent  cone  theory. 

5.  A  new  technique  for  calculating  pressures  on  bodies  at  incidence. 
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FIGURE  7.  TYPICAL  EQUIVALENT  BODY  SHAPES  USED  FOR  COMPUTING 
LIFTING  PROPERTIES  WITH  SECOND-ORDER  SHOCK  EXPANSION 
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The  pressure  gradient  after  a  corner  was  calculated  exactly  through  a  numer¬ 
ical  integration  of  the  equation 


/  8  G  \  y  +  1  Q  d  I 

\dlJ~2~  (m2_1)  '  d6 
Integration  of  this  equation  gives  dG/ds,  and  dp/3s  is  then  computed  from 

8s  \  ds  8s  y 


(38) 


(39) 


where 


2yp 

A  — 

sin  2ji 

and  d6/ds  is  the  curvature  of  the  body  at  a  given  point. 

The  exact  pressure  prediction  gives  much  more  accurate  results  than  does  the 
approximate  method  given  by  Equation  (37),  particularly  for  large  jumps  A0.  As 
noted  by  DeJamette,  it  is  important  to  get  an  accurate  initial  pressure  to  accurately 
predict  the  pressures  over  the  rest  of  the  body  profile.  Similar  expressions  were 
derived  for 


8p 

—  wheno  >  0 
ds 

but  are  not  repeated  here. 

Finally,  the  expression  for  the  pointed-cone  pressure  at  a  >  0  was  estimated  by 
combining  Slender  Body  Theory,  Newtonian  Theory,  and  an  approximate  expression 

=  Cj,  +  ACp  (40) 

a  =  0 


where 


AC 


H 


—  san  2a  sin  20cos4>  +  sin'o  cos* G 


/'  M 

/  2  \  o 

2 - (1  -  tan'O)  - 

1  2  +  —  ).sin“4> 

V  liJ 

{40a) 


sin'O 


1  + 


(y  +  1  ?  K“  +  2  /  y  -r  1 

- In  (  - 

(y-l.K--2  ^  2 


(40b) 
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=  <40c) 

Note  that  4>=0  is  the  leeward  plane  in  this  report  versus  the  windward  plane  in 
Reference  23. 

As  shown  in  Reference  23,  the  addition  of  these  improvements  gave  a 
significant  improvement  in  pressure  and  force  and  moment  prediction  over  a  limited 
range  of  configurations  for  which  computations  were  performed.  Reference  2,  which 
used  the  techniques  of  Reference  23  to  conduct  more  extensive  comparisons  with 
experiment,  found  that  improvements  in  MNT  were  minimal  and,  therefore,  the 
standard  MNT  was  used  with  a  match  point  occurring  where  the  local  Mach  number 
was  1.1.  Hence,  the  latest  version  of  the  NAVSWC  Aeroprediction  Code  (References 
3  and  4)  contains  the  last  three  improvements  by  DeJamette  et  al.  but  does  not 
contain  the  first  two. 

In  more  recent,  unpublished  research,  DeJarnette  was  able  to  derive 
expressions  for  improving  the  Modified  Newtonian  Theory  on  blimt-nosed  bodies 
compared  to  present  usage.  The  derivation  of  this  new  methodology  is  given  in 
Appendix  B.  This  new  methodology  is  given  by 


C  =  C  -kcos“6  cos  5  -cos(5  ) 

P  «!  m 


where  m  =  2.78,  (5eq)^  =  25.95“,  and 


(41) 


k  = 


1 


^  '6  )  <xs“  (6  ) 

eq  m  eq  m 


Here 


(42) 


Thus,  the  theory  that  will  be  extended  to  M»  =  20  and  to  calculate  the  effects  of 
equilibrium  chemically  reacting  flow  behind  the  shock  is  given  by  Equations  (38) 
through  (42). 
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Note  that  making  the  above  extensions  and  modifications  will  require  several 
changes  to  the  existing  Shock-expansicn  Theory  combined  with  MNT  for  blunt 
bodies.  These  include  the  following: 

1.  Calculation  of  real  gas  effects  behind  a  normal  shock. 

2.  Calculation  of  real  gas  effects  across  an  oblique  shock  attached  to  a  cone 
and  wedge. 

3.  Calculation  of  PME  for  chemical  equilibrium  flows. 

4.  Extension  of  the  exact  pressure  gradient  formula  for  a  =  0  and  a  s:  0  to 
include  frozen  and  chemical  equilibrium  flows. 

5.  Modified  Newtonian  Theory  extended  to  real  gases. 

6.  Definition  of  match  point  up  to  =  20. 

7.  Computation  of  temperature  along  body  surface  (in  addition  to  pressure). 

8.  Modification  of  any  empirical  formulas  as  necessary  for  fi^ozen  and  equi¬ 
librium  flows. 

These  changes  will  be  discussed  later  in  this  report.  It  should  be  noted  that,  while  the 
first  three  changes  are  state  of  the  art,24-27  the  last  five  have  never  been  presented 
before  to  the  author’s  knowledge. 


2.3  REAL  GAS  COMPUTATIONAL  PROCESS 


2.3.1  Summary  of  Procedure  for  General  Chemical  Species  Composition 

To  put  real  gas  computations  in  proper  perspective,  consider  the  mathematical 
model  of  an  inviscid  flow  field.  These  equations  are  standard  in  the  literature^lO  and 
are  sometimes  referred  to  as  Euler’s  equations.  They  will  be  stated  here  in  vector 
form  without  derivation. 


Continuity  Equation 

dp 

—  +  V  -  p  V  =  0 

at 

(one  equation,  four  unknowns) 


Momentum  Equation 


(three  equations,  one  additional  unknown) 


(43) 


(44) 
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Energy  Equation 


D{h  +  V^/2)  dp  ^  -:w 

p - - - =  —  +  pq  +  p(f-V) 

Dt  at 

(one  equation,  three  additional  unknowns) 


Equation  of  State 

h  =  h{p,p) 

(one  equation,  zero  additional  unknowns) 


(45) 


(46) 


These  equations  have  few  assumptions.  They  assume  inviscid  flow  with  no 
soiurces  or  sinks  present  and  the  potential  energy  of  the  air  mass  due  to  gravity  is 
negligible.  There  is  a  total  of  six  equations  and  eight  unknowns.  If  we  make  the 
assumption  of  no  body  forces  (i.e.,  f  =  0)  and  no  heat  added  to  the  ^stem  firom  the 
outside  (radiation,  etc.  so  that  q  =  0),  then  the  system  reduces  to  six  equations  and 
six  unknowns.  For  a  perfect  gas.  Equation  (46)  is 


h  =  C  T  (47) 

p 

and  from  the  equation  of  state  for  a  perfect  gas: 


T  = 


P 

pR 


(48) 


Here  Cp  is  a  constant.  Equation  (48)  alci^  with  Equations  (43)  through  (47)  form  a 
deterministic  mathematical  model  for  general  3-D  computation  of  inviscid  flows  over 
configurations  in  a  Mach  number  range  where  these  equations  are  valid.  They  are 
solved  in  finite  difference  form  such  as  References  28  and  29  or  in  various  approx¬ 
imate  forms  such  as  References  3, 4,  or  30-32  (among  others).  If  viscous  flows  are  of 
interest,  the  viscous  stress  terms  can  be  added  to  Equation  (44),  and  the  set  of 
equations  then  becomes  the  Navier  Stokes  equations.  Also,  they  can  be  solved 
numerically  or  by  approximations  of  a  thin  layer  near  the  body  (boimdary  layer). 


The  question  we  must  answer  is  how  we  calculate  the  enthalpy  (or  internal 
energy  since  h  =  e  +  p/p)  when  the  gas  is  at  a  temperature  such  that  the  perfect  gas 
law  is  no  longer  valid.  Numerous  efforts  are  reported  in  the  literature  going  back  to 
the  1940s  and  1950s  when  high-speed  flight  was  really  becoming  a  practical  problem. 
Again,  a  rederivation  of  results  will  not  be  repeated  here — only  a  summary  com¬ 
putation  process.  The  process  draws  heavily  on  the  kinetic  theory  of  gases  and 
Quantum  Mechanics. 
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The  internal  energy  of  an  atom  (O,  N,  0 + ,  N + )  or  molecule  {NO,  O2,  N2)  is 
comprised  of 


e  =  +  e..  +  -f  (49) 

where  subscripts  t,  r,  v,  e,  d  stand  for  translation,  rotation,  vibration,  electronic 
excitation,  and  dissociation.  For  a  single  atom,  only  et  and  eg  are  important  since  the 
rotational  energy  is  small  compared  to  translational  energj'.  From  the  kinetic  theory 
ofgases^ 


for  both  atoms  and  molecules.  Also 


(50) 


e  =0  for  atoms 
r 


e  =  RT  for  molecules 
r 


(51) 


The  vibration  energy  is  somewhat  more  difficult  to  compute.  The  vibration  energy  is 

e  =  0  for  atoms 


(52) 


where  Gy  is  2270°K  for  O2, 3390‘'K  for  N2,  or  approximately  3160®K  for  air  up  to  the 
point  of  dissociation.  Also,  exp  has  been  used  to  denote  exponential  to  distinguish  the 
normal  terminology  of  e  for  exp  from  the  total  internal  energy  e. 

Combining  Equations  (49)  through  (52),  for  air  up  to  approximately  2000°K 


e  =  e  +  e  +  e 

t  r  V 


Re 


e=  -  RT  + 
2 


r'\ 

exp  5  —  ?  — 


V  T  / 


(53) 


Equation  (53)  is  for  molecules  only  as  we  assume  air  is  about  80  percent  N2  and  20 
percent  O2  and,  therefore,  no  atomic  species  are  present.  Note  that  Equation  (53)  is  a 
function  of  several  fixed  constants  and  temperature  only.  Hence,  Equation  (1)  or  (46) 
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holds  identically,  and  one  can  replace  Equation  (48)  with  Equation  (53)  for 
temperatures  in  the  range  of  BOO^K  to  2000®K.  Figure  8  is  a  plot  of  the  specific  heats 
and  the  ratio  of  specific  heats  where  only  the  rotation,  translation,  and  vibration 
modes  of  energy  are  accounted  for  (see  References  6  and  8). 


0.01  .6  1  2  3  4  5  6  7  8 


LOCAL  TEMPERATURE  X 10-3  (K) 

HGURE  8.  SPECinC  HEATS  FOR  AIR  (NO  DISSOCIATION  OR  IONIZATION) 


As  air  is  heated  above  2000‘’K,  it  begins  to  dissociate  into  N,  0,  and  NO  atoms 
and  molecules  in  addition  to  the  N2, 02  already  prevalent.  Initially,  the  O2  begins  to 
dissociate,  but  at  temperatures  of  about  4000°K,  the  N2  dissociation  also  becomes 
appreciable.  Finally,  above  temperatures  of  9000°K,  the  oxygen  and  nitrogen  atoms 
and  molecules  begin  to  give  up  electrons  (ionize).  This  dissociation  and  ionization  of 
the  air  molecules  creates  additional  internal  energy  that  must  be  accounted  for  to 
accurately  predict  the  thermodynamic  and  transport  properties  in  the  flow.  To 
compute  these  thermodynamic  properties,  we  first  must  determine  the  species 
concentrations  at  the  given  T  and  p  of  the  mixture.  Then  enthalpy,  internal  energy. 
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and  values  of  specific  heats,  ratio  of  specific  heats,  and  gas  constant  can  be  computed. 
Some  typical  chemical  reactions  for  air  are 


02^5  20 

(54) 

N2±;2N 

(55) 

N  +  0±5N0 

(56) 

N  +  Oi^NO^  +  e~ 

(57) 

Reference  11  shows  there  is  a  total  of  22  chemical  reactions  that  occm-  for  temper¬ 
atures  up  to  15000“K,  resulting  in  a  set  of  28  nonlinear  equations.  Each  of  these 
reactions  will  have  an  equilibrium  constant.  These  equilibrium  constants  are  known 
functions  of  temperature  fi’om  either  measurements  or  statistical  mechanics 
computations.  These  equilibrium  constants  are  related  to  the  partial  pressures  of  the 
individual  constituents  by  the  Law  of  Mass  Action.5  Here  the  partial  pressures  can 
be  thought  of  (for  example,  P(j)  as  the  pressure  that  would  result  if  a  container  of 
given  volume  were  filled  with  only  oxygen  atoms.  As  an  example,  for  Equations  (54) 
through  (57),  the  equilibrium  constants  are  related  to  the  partial  pressures  of 
individual  constituents  as  follows: 


(58) 


(59) 


Pn^o 


(60) 


PnoP# 

PnPo 


=  K  (T) 

P.NO"^ 


(61) 


There,  of  course,  would  be  a  single  equation  corresponding  to  each  significant 
chemical  reaction.  For  temperatures  below  about  9000°K,  Equations  (58)  through 
(61)  are  the  most  important  reactions.  In  addition  to  Equations  (58)  through  (61),  we 
have  Dalton’s  Law  of  Partial  Pressures,  which  states  that  the  total  pressure  of  the 
mixture  is  the  sum  of  the  partial  pressures.  Mathematically,  for  the  constituents  of 
Equations  (58)  through  (61),  this  can  be  wTitten  as  follows: 

P  Po  +  Po  Pn,  Pn  Pno  +  P„„^  -^-P  -  ^®2) 

2  2  NO  e 
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Of  course,  if  the  temperature  is  higher  than  9000'’K  and  other  chemical  reactions 
occur,  the  other  species  must  be  added  to  Biquation  (62),  just  as  the  chemical  reaction 
equations  are  included  in  the  set  (58)  through  (61).  Note  that  there  are  seven 
unknowns  in  Equation  (62)  and  five  equations  so  far.  The  remaining  two  equations 
come  from  a  chemical  balance  of  the  number  of  O  and  N  atoms  and  the  fact  that 
electric  charge  must  be  conserved.  These  facts,  in  equation  form,  may  be  written  as 


0  2  ^  "  "no 

—  =  —  =  0.25  = - 

0.8  2p^  ^Pn-Pno+Pno 


(63) 


and 


P  =  P  _  (64) 

.so*  e 

Again,  if  more  reaction  equations  are  included,  then  Equations  (63)  and  (64)  will 
change  to  include  other  species.  The  fundamental  principle  of  correct  balance  of  O2, 
N2  and  electric  charge  will  not  change,  however.  For  a  given  temperature  and 
pressure.  Equations  (58)  through  (64)  give  a  unique  set  of  seven  partial  pressures  for 
the  seven  species  (Reference  5,  of  course,  has  22  species).  This  is  a  nonlinear  set  of 
algebraic  equations  that  must  be  solved  for  each  given  T,  p  to  determine  the  species 
content  once  temperatures  are  high  enough  for  dissociation  to  occur.  Knowing  the 
partial  pressures  of  the  species  present  allows  one  to  uniquely  determine  the  other 
properties  of  the  system.  The  details  will  not  be  repeated  here  but  can  be  found  in 
References  5, 9,  or  33.  Figure  3  is  an  example  of  an  equilibrium  composition  of  air  as 
a  function  of  temperature  at  a  given  density. 

2.3.2  Simplified  Procedures  for  .^ir 

This  discussion,  thus  far,  has  focused  on  how  to  compute  the  thermodynamic 
and  flow  field  properties  of  a  chemically  reacting  gas.  For  general  chemical 
reactions,  this  is  the  process  that  must  be  u.sed.  However,  for  air  other  alternatives 
are  available.  For  example,  Reference  11  has  produced  tables  of  the  properties  of  air 
as  a  function  of  p  and  p.  These  could  be  used  in  a  table  look-up  mode  in  the 
computation  process.  A  more  simplified  approach  was  produced  by  Tannehill  and 
Mugge34  and  later  extended  to  other  thermodynamic  variables  by  Srinivasan  et  al.35 
They  produced  curve  fits  from  algebraic  equations  for  p  =  p(e,  p),  a  =  a(e,  p),  T  =  T(e, 
p),  h  =  h(p,  p),  and  T  =  T(p,  p).  These  fits  are  valid  up  to  temperatures  of  25000°K. 
Since  this  is  by  far  the  most  computationally  efficient  method  and  since  we  will  only 
be  dealing  with  air,  this  is  the  technique  used  in  this  report.  Since  it  is  the  technique 
used,  a  brief  discussion  of  the  equations  relevant  to  the  work  herein  will  be  given. 

Appendix  B  shows  the  iteration  process  to  determine  properties  behind  a 
normal  or  oblique  shock  wave,  which  involves  assuming  a  value  of  pi/p2  behind  the 
shock,  computing  p2  and  h2  from  this  assumed  value,  and  then  recomputing  a  new 
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value  of  p2  based  on  these  values  of  p2  and  h2  from  the  thermodynamic  properties  of 
equilibrium  air.  An  example  of  the  algebraic  equation  contained  in  Reference  34  is 
shown  here.  From  Reference  34,  an  effective  y  as  a  function  of  p  and  p  is  defined  by 


y  =  C ,  +  C.,Y  +  az  +  C  YZ  + 

i-  ^  J  4 


C.  +  C  Y+C,Z  +  C,YZ 

3  6  I  6 

l  +  expi(Cg{X  +  CjgY  +  C,j)) 


(65) 


where 


‘“^10(1.292; 


X  =  log 


I - 5 - ') 

1.013x10=*^ 


Z  =  X- Y 

The  coefficients  Ci.  C2,  C3,  C4,  C5,  Cs,  C7,  Cs,  C9,  Cio,  Cn  are  tabulated  in  Table  1  of 
Reference  34.  This  allows  one  to  compute  entiialpy  from  the  relation 


h  = 


P 

P 


and  energy  from 

e  =  h-^ 

P 

Also,  from  values  of  p  and  p,  the  temperature  can  be  determined  from 


(66) 


(67) 


=  dj  +  d2Y  +  d3Z  +  d_,YZ  +  d.Z^ 
d.  +  d.  Y  +  d^Z  +  d„  Y  Z  +  d,„Z^ 

6  4  b  9  10 

1  +  expldjj(Z  +  d^  J 


(68) 


Here 


Y  =  loK,„ 

^10 


X  =  log 


10 


'  1.0134  X  10 


i) 
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Z  =  X- Y 

and  di,  d2,  da,  d4,  ds,  de,  d?,  dg,  dg,  dio,  dn,  di2  are  given  in  Table  2  of  Reference  34. 
Knowing  p,  p,  T,  h,  y,  e,  the  only  quantities  remaining  are  velocity  and 
compressibility  factor.  These  can  be  determined  from  the  constantcy  of  total 
enthalpy  relation  and  the  equation  of  state 

uj  vl 

H  =  Const  =  h,  +  —  =h^+ —  (69) 


Z  = 


P 

pRT 


(70) 


Knowing  H  from  freestream  conditions  and  a  new  h2  allows  one  to  compute  the  local 
U2.  It  should  be  noted  here  that  the  relationship 


y  p 

p 

does  not  yield  the  true  speed  of  soimd  as  in  the  case  of  a  perfect  gas.  This  is  because  y 
is  h/e  as  defined  by  Reference  34  and  y  =  y  only  if  no  dissociation  of  air  molecules 
occurs.  Tannehill34  gave  an  expression  for  the  approximate  speed  of  sound  if  e  and  p 
were  known.  That  is 


a  = 


dy 


This  is  also  provided  in  curve  fit  form  by  Tannehill 


(71) 


a  =  a  (e,  p) 

The  correct  local  Mach  number  is  then 


(72) 


Srinivasan35  produced  additional  curve  fits  of  various  thermodynamic 
variables.  These  included 

p  =  p(h,p}  ^ 

P  =  p{p,T) 


S  =  S  (e,  p) 
p  =  p  (p,  S) 
a  =  a  (p.  S) 
e  =  e(p,S) 


J 


(73) 
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Having  curve  fits  available  for  thermodynamic  variables  means  that  the 
computational  time  for  real  gas  flow  fields  has  been  reduced  considerably  over  what 
it  would  be  if  the  chemistry  of  the  gas  had  to  be  computed  at  each  point  in  the  flow 
field. 

Some  results  of  the  Reference  34  curve  fit  procedure  are  given  in  Figures  9 
through  12  in  terms  of  various  thermodynamic  properties.  In  each  of  the  figures,  the 
thermod\'namic  properties  are  referenced  to  seme  set  of  standard  conditions  that 
could  be  the  sea  level  standard  atmosphere  for  example. 


PRESSURE  RATIO,  p/po 


FIGURE  9.  CURVE  FIT  EQU.ATIOXS  tL\TERN.4L  ENERGY  RATIO* 
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FIGURE  11.  CURVE  FIT  EQUATIONS  (TEMPERATURE  RATIO) 
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2.3.3  Chemical  Reaction  Rates 

In  the  preceding  discussion,  it  was  assumed  that  the  reaction  rates  of  the 
chemical  processes  were  infinitely  fast  so  that  the  flow  could  be  assumed  to  be  in 
chemical  equilibrium  at  all  points  in  the  flow  field.  If  one  defines  the  characteristic 
length  of  the  flow  field  to  be  the  body  length,  then  the  time  it  takes  for  a  fluid  element 
to  traverse  the  shock  wave  and  full  body  length  is 


Defining  the  time  it  takes  for  a  chemical  reaction  to  come  into  equilibrium  as  to,  then 
if 


t,>>t2  (75) 

the  flow  is  chemical  equilibrium  for  all  practical  purposes.  On  the  other 

hand,  if 


t.  >  >  t. 


(76) 
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we  may  assume  that  the  flow  is  “frozen”  at  the  chemistry  that  exists  just  behind  tn? 
shock.  That  is,  we  can  compute  the  flow  as  if  it  were  a  perfect  gas  except  we  have  a 
different  value  of  y  =  yj. behind  the  shock.  Finally,  if 

=  (77) 

we  have  nonequiiibrium  flow,  which  means  that  the  chemical  processes  are  occurring 
and  changing  significantly  along  the  body.  This  results  in  an  order  of  magnitude 
more  difficulty  in  making  real  gas  computations.  What  can  be  done,  however,  from  a 
more  practical  standpoint  is  to  compute  the  Iccal  conditions  assuming  “frozen”  and 
“equilibrium”  flows  and  assume  that  the  nonequilibrium  case  falls  in  between  these 
two.  That  is,  generally  the  equilibrium  and  frozen  flow  cases  present  boundaries  for 
flow  properties  and  the  nonequilibrium  case  fails  in  between. 

Since  the  focus  of  the  present  work  is  to  incorporate  real  gas  effects  into  SOSET, 
only  equilibrium  and  frozen  flows  will  be  considered.  The  computations  involved 
v-ith  a  nonequilibrium  flow  are  inconsistent  with  an  approximate  mathematical 
model  of  the  flow  field  and  are  more  consistent  with  solution  of  the  full  Euler  or 
Navier  Stokes  equations. 


2.3.4  Modified  Newtonian  Theory 


Recall  from  Section  2.2.1.  that  the  MNT  pressure  coefficient  is  given  by 


C 

? 


=  C 


•  2c 

Sin  5 


*^0 


e<j 


(78) 


where  5eq  is  evaluated  according  to  Appendix  A.  (Note  that  even  though  Equations 
(41)  and  (42)  are  the  relations  used  to  calculate  pressures  around  the  blunt  portion, 
they  are  both  functions  of  modified  Newtonian  pressures.)  For  the  perfect  gas 
computations,  Cp^  was  evaluated  according  to  Equation  1 12 However,  Equation 
(12)  no  longer  he  Ids  for  real  gas  computations  across  normal  or  oblique  sh.Kk  waves 
(see  Appendix  C).  The  conditions  behind  the  shock  wave  can  no  longer  be  solved  for 
in  closed  form,  as  in  Section  2.2.1,  because  the  specific  heats  and  compressibility 
factor  are  no  longer  constant.  This  means  that  an  iteration  takes  place  to  determine 
the  properties  immediately  behind  the  shock.  Appendix  C  thus  is  the  process  for 
defining  conditions  behind  the  shock,  i.e. 


>■■>•  P?  Py  ^2-  ^2’  ^2 

The  question  that  must  be  addressed  is  twofold; 

1.  How  do  these  parameters  vary  along  the  stagnation  streamline  between 
the  body  and  shuck? 

2.  How  do  these  properties  vary  from  the  stagnation  poir. :  around  the  body  to 
the  point  where  shock  expansion  ih^jry  is  applied? 
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This  section  of  the  report  will  address  these  questions  for  both  frozen  flow  and 
equilibrium  chemically  reacting  flows  where  MNT  is  used  on  the  blunt  portion  of  a 
configxu'ation  for  pressure,  temperature,  force,  and  moment  computations. 

2.3.4.1  f  rozen  Flow.  For  frozen  flow,  it  is  assumed  that  the  r  and 
compresSibil  jsy  factor  ve  fixed  a^  the  values  thej'  have  just  behind  the  shock  (i.e.,  y 
=  y,,  Z  —  Z|).  It  IS  also  known  that,  along  the  stagnation  streamline,  the  flow  is 
isentropic  equilibrium  flow  conditions.  Furthermore,  it  is  known  from  many 
orevious  Works  (ser,  for  example.  References  10, 13,  and  24)  that  the  Mach  number 
behki^.  a  cormal  shock  is  much  less  than  1  for  high  freestream  Mach  numbers.  Also, 
this  Mach  number  must  go  to  0  at  the  stagnation  point  which  is  only  a  slight  distance 
away  from  the  shock.  This  has  led  to  approximate  theories  in  the  stagnation  region 
based  on  constant  density,  temperature,  and  pressure  which,  in  effect,  is  the 
Newtonian  assumption  in  the  shock  layer.  That  is,  the  shock  is  assumed  to  lie  on  the 
body  and  the  conditions  behind  the  shock  are  the  same  as  those  on  the  body  surface. 
The  physics  of  the  flow  between  the  shock  and  body  is  what  allows  one  to  assume 
and  Zf  are  in  fact  approximately  constant  for  both  equilibrium  and  frozen  flow 
between  the  shock  and  body  and  along  the  stagnation  streamline.  Reference  24 
indicates  a  maximum  2.3  percent  error  in  making  this  assumption,  and  in  most  cases 
the  error  was  much  less. 

A  second  assumption  will  also  be  made.  This  assumption  will  allow  the  use  of 
Equations  (9),  (14),  and  (15)  in  a  local  sense  along  the  body  where  y  =  y,  and  M  is  the 
local  value.  Strictly  speaking,  these  equations  relating  total  to  static  conditions  are 
derived  based  on  perfect  gas  assumptions  and  with  constant  y  =  1.4.  The  accuracy  of 
this  assumption  (combined  with  others)  can  be  assessed  in  comparison  with  exact 
calculations  over  a  blunt-nosed  body  by  solution  of  the  full  Euler  equations.  This  will 
be  done  in  the  Results  section  of  thi.s  report. 

Returning  now  to  Equation  ( 78)  and  recognizing  the  fact  that  all  the  properties 
behind  the  shock  are  known  from  the  real  geis  computations  of  the  flow  across  a 
normal  shock  (Appendix  C),  one  can  write 


along  stagnation  streamline.  Here  pi,  p2:  h2,  V2,  a?,  M2,  >/ ,  Zf  etc.  are  the  values 
calculated  immediately  behind  the  shock  from  Appendix  B.  Also 
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^2 

H  =  h2+Y 


=  T,(l  + 


00,  =  02(1  + 


(81) 

(82) 


(83) 


Now,  since  at  the  stagnation  point  local  velocity  and  Mach  number  are  zero,  the  static 
pressure,  temperature,  density,  and  enthalpy  at  the  stagnation  point  are  simply  the 
total  values  behind  the  shock  given  by  Equations  (80)  through  (83).  Thus,  at  the 
stagnation  point,  for  &ozen  flow 


hs  =  H 

Ps  =  P02 

Ms  =  Vs  =  0 

y  =  Vj.  =  Const 

C  =C  =  Const 
f 

C  -  C  =  Const 

Zj.  =  Zj  =  Const 
5  =  82  =  Const  a!cn^  body 


(84) 


To  continue  the  computations  around  the  blunt  surface,  an  isentropic  expansion  is 
a^umed  and  Equation  (78)  is  applied  at  each  point  on  the  surface  where  8  is  given  by 
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Appendix  A  and  Cp^  is  constant  at  the  value  given  by  Equation  (79).  Knowing  8eq, 
CpQ  and  Cp  allows  one  to  calculate  the  local  pressure  at  a  given  point.  That  is 


Pl  =  p« 


1  + 


(85) 


Since  the  flow  is  frozen,  we  can  treat  it  just  like  a  perfect  gas  except  that  the  values  of 
y  =  yf  and  Rf  =  RZf.  Kence,  knowing  local  pressure  from  Equation  (85),  total 
pressure,  temperature,  density,  and  other  flow  properties  behind  the  normal  shock 
[Equations  (80)  -  (84)1,  the  fact  that  total  conditions  are  constant  along  a  streamline, 
the  local  conditions  at  any  point  on  the  blunt  body  can  be  calculated  by 


2.3.4.2  Equilibrium  Flow.  Equilibrium  chemically  reacting  flow  properties  at 
the  stagnation  point  are  computed  just  like  the  frozen  flow  case.  This  is  based  on  the 
rationale  discxissed  previously  of  only  slight  changes  in  the  specific  heats  and 
compressibility  factor  between  the  body  and  shock.  However,  for  the  isentropic 
expansion  around  the  body,  >•  and  Z  are  allowed  to  vary  from  the  values  at  the 
stagnation  point  in  contrast  to  the  frozen  flow  case  where  these  parameters  are  frozen 
at  values  behind  the  shock.  Hence,  for  equilibrium  chemically  reacting  flows. 
Equation  (84)  defines  the  stagnation  point  conditions.  Also,  Equation  (85)  defines 
the  pressure  at  the  next  point  around  the  body.  To  get  the  remaining  properties, 
recourse  is  once  again  made  to  the  real  gas  curve  fits  of  Section  2.3.2.  To  obtain  p,  we 
know  local  pressure  and  entropy  so  that 
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Pl  =  p(p.S)  (92) 

Then,  knowing  pressure  and  density 

h  =  h  (p,  p),  =  T  (p,  p),  =  a  (p.  p)  (93) 

The  remaining  quantities  can  be  computed  from  thermodynamic  relations  as  fellows: 


VL=^2(H^-h,) 

(94) 

(95) 

(96) 

(97) 

(98) 

This  process  is  continued  around  the  body  surface  to  the  point  where  shock  expansion 
begins.  Again,  Equations  (85),  (86),  and  (92)  through  (98)  provide  initial  conditions 
for  the  shock  expansion  process. 


2.3.5  Second-order  Shock-expansion  Theory  (SOSET) 

2.3.5. 1  Zero  Angle-of- Attack  Solution 

Recall  that  SOSET  was  given  by  Equations  (35)  and  (36),  repeated  here  for 
convenience. 

P  =  Pc  -  (Pc  -  P2’®  '’  ^35) 


n  = 


(36) 


We  desire  to  develop  the  theory  to  allow  computations  to  proceed  around  the  body 
using  modified  versions  of  £)quations  (35)  and  (36)  for  real  gas  effects.  To  do  this 
requires  pressures  on  wedges  and  cones  calculated  based  on  real  gases  (see  Appendix 
C)  and  pressures  behind  an  expansion  comer  computed  by  the  PME  for  real  gases  (see 
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Appendix  D).  Note  some  of  the  results  of  wedge,  cone,  and  expansion  flows  of  real 
versus  perfect  gases  shown  in  Appendixes  C  and  D.  The  pressures  in  compression 
processes  are  almost  identical  for  perfect  and  real  gases,  whereas  the  temperatures 
show  significant  differences  as  the  Mach  number  and  compression  angle  increase.  On 
the  other  hand,  expansion  processes  show  differences  for  both  pressure  and  temp¬ 
erature  of  real  versus  perfect  gas  computations. 

If  one  were  only  interested  in  applying  generalized  Shock-expansion  Theory  for 
equilibrium  or  frozen  flows  over  wing  body  configurations  in  a  local  sense  at  near¬ 
zero  angle  of  attack,  one  would  stop  at  this  point.  Appendixes  C  and  D  give  the 
results  for  pressures  and  temperatures  that  are  used  for  points  aromad  the  2-D  or 
axisymmetric  body.  One  would  simply  compute  the  pressure  and  temperature  on  the 
blunt  portion  of  the  body  using  section  2.3.4,  determine  a  match  point,  use  the  PME  of 
Appendix  D  to  march  around  the  surface  where  the  pressure  and  other  properties  are 
constant  on  each  straight-line  segment  of  Figure  6.  For  a  sharp  conical  tip  or  wedge, 
results  of  Appendix  C  could  be  used  for  the  initial  solution,  and  the  solution  could 
proceed  along  the  body  or  wing  siu’face  again  using  the  PME  results  of  Appendix  D. 
However,  experience  has  shown  that  a  first-order  solution  is  unacceptable  for  pres¬ 
sure  on  most  bodies  of  revolution  at  lower  supersonic  Mach  numbers,  so  a  second- 
order  accurate  pressure  scheme  is  necessary  for  good  force  and  moment  predictions. 
Since  this  is  the  case,  a  comparable  second-order  technique  for  real  gas  pressures  and 
temperatures  at  the  surface  to  be  available  for  inputs  to  heat  transfer  and  force  and 
moment  computations  is  also  desirable. 

The  crux  of  the  problem  in  addressing  second-order  solutions  for  inclusion  of 
real  gas  effects  is  to  successfully  address  the  pressure  gradient  of  Equation  36 


behind  the  comer  and  to  find  simple  but  accurate  ways  of  addressing  angle-of-attack 
computations  for  values  of  Pc  in  Equation  (35)  and  temperature  along  the  surface. 
The  pressure  gradient  change  behind  a  comer  will  be  investigated  first.  To  do  this 
will  require  several  steps: 

1.  Conversion  of  the  equations  of  motion  [Equations  (43  and  (44)1  from 
rectangular  to  streamline  coordinates  and  derivation  of  the  characteristic 
equations. 

2.  Ocrivi  cio.i  of  the  pressure  change  along  a  left,  mnning  characteristic  (or 
Mach  line)  since  this  is  the  mechanism  for  differential  pressure  change 
along  the  surface. 

3.  For  angie-of-attack  compulations,  a  method  for  accurately  computing  Ap 
and  AT  due  to  a. 
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4.  An  algebraic  approximation  as  a  solution  to  an  ordinary  differential 
equation  (ODE). 

While  the  labor  and  math  are  fairly  involved  to  get  to  step  4,  the  solution  of  the 
equation  in  step  4  is  the  only  thing  the  computer  sees  and  it  is  very  fast.  Most  of  the 
details  of  the  derivation  of  steps  1-2  will  be  indicated  in  the  following  discussion  and 
will  be  included  either  in  the  text  or  an  appendix.  Angle-of-attack  effects  are 
discussed  in  section  2.3.6.  The  theoretical  methodology  follows  the  approach  of 
Reference  23  for  perfect  gases  and  modifies  the  mathematical  model  as  necessary  or 
states  assumptions  made  to  allow  real  gas  computations. 

Appendix  E  converts  the  continuity  and  momentum  equations  from  rectangular 
to  streamline  coordinates  and  then  derives  the  characteristic  compatibility  relations. 
These  equations  for  axisymmetric  bodies  at  zero  angle  of  attack ,  repeated  here,,  are 

JL  -  ~  sin6.sinti 

pV2  dC,  dC^  ■  r 


jj  dp  d6  —  sin0sinp 
pV2  dC2  dC^  r 


(IOC) 


Since  we  are  looking  for  dp/ds  downstream  of  an  isentropic  turn,  information  away 
from  the  body  surface  must  be  obtained.  To  do  this,  we  must  relate  the  change  in  flow’ 
variables  along  the  left  running  characteristic  to  the  change  of  these  variables  along 
a  streamline.  In  general,  one  can  write 

.s  =  slCj.Cj) 


w’here  s,  Ci,  and  C2  are  all  curvilinear  coordinates.  By  the  chain  rule 


d 

ds 


dCj  ^  dC.,  ^ 
ds  dCj  ^  ds  dC, 


Also,  the  differential  ds  is 


ds 

ds  =  —  dC  - 
dC.  1 


ds 


(101) 


(102) 


And,  since  ds  =  dCi  cos  p  +  dC2  cos  p  and  along  a  streamline  dCi  =  dC2,  then 


dCj  dC, 


1 


ds  ds  2cosp 


36 


NAVSWC  TR  90-683 


Thiis  Equation  (101)  becomes 


or 


1.  -  ^  ^ 

3s  2cosp  3Cj  dC^ 


=  2  cos  p  — 
ds 


(103) 


(le-y 


Now  adding  Equations  (99)  and  (100),  using  the  relation  (104),  there  is  obtained 
(where  A  =  pV^/p) 


—  an  p  sin  6 
r 


cosp 


dp 

ae 

_  \ _ 

ds 

A 

ds 

(105) 


Likewise,  subtracting  Ekjuation  (99)  from  (100),  the  relationship  for  pressure  along 
the  left  running  characteristic  is  obtained  as 


dp  (  dp  36 

dC^  j  ds  ds  1 


(106) 


Equations  (105)  and  (106)  are  the  two  relations  we  will  xise  to  provide 
information  to  allow  computation  of 


\  ds  4 

To  do  this,  consider  a  streamline  very  near  the  body  surface  that  will  form  a 
streamtube  between  the  body  and  the  streamline  (see  Figure  13).  Referring  to  Figure 
13,  the  distance  along  a  Mach  line  to  the  streamline  is  b,  which  again  is  assumed  to 
be  very  small  compared  to  the  body  radius.  The  thickness  of  the  streamtube  is 
b  J  sin  p,  ahead  of  the  turn  and  b2  sin  p2  after  the  turn. 


MACH  UNES 
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Applying  Equation  (106)  along  the  streamline  near  the  body  between  points  ® 
and  (D,  wMch  correspond  to  the  beginning  and  end  of  the  expansion  turn,  there  is 
obtained 


fPs  dp 

^P  •  4. 

ae 

J  4  A  cosp 

-  QS  + 

3C^ 

1  —  d  s 

J4  ds 

(107) 

The  left  side  of  Equation  (107)  may  be  expanded  as 


p!!p=  f  ilLac  *  f 

ip  A  J4  A  ac,  1  ii  A  J2  A  dC,  1 


(108) 


The  first  and  third  integrals  on  the  right  side  of  Equation  (108)  can  be  approximated 
by 


fi  IL 

\A  aCj 


(109) 


Here  1-4  and  5-2  refer  to  the  average  values  of  these  parameters  between  these 
points.  The  differential  dCi  between  1-4  and  5-2  is  simply  the  distance  along  the  left 
running  characteristics,  which  we  defined  as  bi  and  b2.  Note  that  since  the  first 
integral  goes  from  4  to  1,  the  integral  of  dCi  is  negative  so  it  integrates  to  -bi. 

The  second  integral  on  the  right  side  can  be  integrated  directly  using  the 
Prandtl-Meyer  Expansion.  Referring  to  Appendix  C,  where 


IP  =  + 

<1®  Vm^-i 


then 


Substituting  Eqiiations  (109)  and  (110)  into  (108)  there  is  obtained  for  the  left  side  of 
Equation  (107) 


Jp^  A  UA  dCj  /i_4 


(111) 
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Assuming  that  the  streamline  is  close  to  the  body  sur&ce  so  that  the  Mach  lines 
and  Jig  basically  straight  and  the  pressure  varies  linearly  along  the  left  running 
characteristic,  then  the  second  integral  of  Eqiiation  (107)  may  be  approximated  by 
use  of  the  trapezoidal  rule  as 


J4  Acosp  Aj  cospj  V  dCj /j  A^cosp^^^C/j 

The  last  integral,  of  course,  integrates  directly  to 

f  d0  =  e.-0, 

14  = 

Again  expanding  in  a  Taylor  series  expansion  yields 


(s.-s  ) 
o  4 


0=0  +  +■•• 
or 

Substituting  (111),  (112),  and  (113)  into  (107)  one  obtains 


^  V  A  ^1-4  *-  I  A  \  dCj  /  5-2  I  AjCos  V  dC^  /j 


Pg 


f— ) 

\  8C.  A 


} 


+  b. 


/  80  \  /  50  \ 

:l8Cj/2 


(112) 


(113) 


(114) 


Note  that  62  -  61  cancels  on  each  side  of  Equation  (114).  Now  divide  through  by  bi 
and  take  the  limit  as  bi  -♦  0  to  obtain 


\  A  8C.  /I  V  bj  /  \  A  8C.  Jz  i  Ajcos  p^  V  aC^  A 


,  ‘  P”  H  f  (115) 

A^cos  p^  '8Cj/2  I  2bj  b.  laC^A  '  5C,  A 

Note  that  b2  /  bi  and  (S5  -  S4)  /  bi  are  not  infinite  but  in  fact  are  finite  as  will  be  shown 
shortly.  Now  Equations  (105)  and  (106)  can  be  used  to  evaluate 
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Substituting  for  these  relations  and  carrying  out  the  algebra,  Equation  (115)  becomes 


2  cos 


dp  36 

—  -  A  — 
ds  3s 


sin  (pj)  ^  (6^) 


ba  j  200SP2  I  3p  30  j  sin  (p^)  sin (e^) 

X  ' 

ds  ds  12  T 


b,  I  A, 


.  I  ^.U^, {!£.,£?)  1 

I  Aj^cos  <Pj‘  t  1  \  3s  3s  /j  I 


f  dp  30  A  ]  ]  ®5  ®4 
cos(pJ  —  -  A  — 

^  '  3s  3s  /2  ]  J  2  b, 


A^cos  (p^) 


(116) 


Now  except  for  the  matchpoint  between  a  blunt  nose  and  a  conical  frustrum 

(-)  =(-)  =  0 
'  3s  /I  '3s  1 2 

since  0  is  constant  along  conical  segments.  At  the  matchpoint 

/dex  1 

(  ~  R 


(117) 


where  Rjf  is  the  radius  of  the  spherical  nose  tip.  Using  the  relation  (117)  and  solving 
for 


(si 


one  obtains  for  Equation  (116) 


\  3s  /2 


/  1  3p  A  -^5  I  1  I  \ 

I - )  -  +4co3(u  )  j  +  -  jan(p  )sin(9  ) -(  —  ).sin{p  )sn(8  ' 

'  2  A  ds  /I  I  b,  ^  ]  T  {  ^  "  \h,  ’  ^  ^ 

1  i 


(118) 


/b,'2ccs(p2)  Sj-s^ 

'  ^  ‘'2 
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From  geometry,  referring  to  Figure  13  and  the  sketch  below,  one  can  write 
approximately  for  small  angles  e 


(s.  -  =  ^2  ~  *^^^2 


This  can  be  written  as 


(119) 


From  conservation  of  mass,  one  can  write 

PjVjbjSin  (Pj)  =  p^V^b^sin  (p,? 

or 

l>2 

l>j  PgVjSin  (p^l  (120) 

Using  (120)  in  (119)  there  is  obtained 
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«5  ~  ^  _  \2 


PjV,s;n  (p.) 
Pj^aSin  (Pa^ 


(121) 


Note  that  with  Equations  (120)  and  (121),  Equation  (118)  can  be  calculated 
independent  of  the  distances  hi  and  82.  The  only  distance  that  enters  the  equation  is 
r,  which  is  the  local  body  radius  at  a  point  in  question.  Defining 


and  ^2~  ^ 


b 


1 


where  these  expressions  are  determined  by  Equations  (120)  and  (121),  and 


F 


3 


(S).J 


‘'*=7 


sin  ( pj)  an  (0  j)  -  sin  (p^  sin  (G^) 


F 


5 


2  cos 


(Pa) 


allows  E-quation  (118)  to  be  written  in  an  abbreviated  form  as 


/  £p  \  _  ^3“^  ^4 

Ids  4  “  F5-F^/2.\2 


(122) 


This  is  our  final  expression  we  desired  to  relate  the  pressure  gradient  behind  an 
expansion  point  for  use  in  the  SOSET. 
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Note  that  Equation  (122)  gives  a  first-order  estimate  for 


\  ds  > 


ds  ^2 


based  on  all  known  quantities  once  the  Prandtl-Meyer  integration  of  Appendix  D  has 
taken  place  through  each  turn.  This  equation  is  valid  for  real  gases  since  no  perfect 
gas  assumption  has  been  made. 

Returning  now  to  Equations  (35)  and  i36),  use  of  Equation  (122)  in  conjunction 
vihth  the  Prandtl-Meyer  solution  of  Appendix  D  now  allows  one  to  estimate  the 
pressure  ail  along  the  body  surface  to  second-order  accuracy  for  real  gases.  With 
entropy  constant  along  the  body  siirface,  other  thermodynamic  and  flow  field 
quantities  can  be  calculated  from  Equations  {D-27)  and  (D-28)  along  with 

T  =  T(p,S) 


M  =  V/a 


(123) 


y  = 


o 

a  p 


Z  = 


P 

pRT 


Previous  applications  of  SOSET  have  applied  the  pressure  computed  from  Equation 
(35)  at  the  midpoint  of  all  body  sections  (see  Figure  6)  for  force  and  moment 
computations.  While  this  is  the  correct  location  for  2-D  siu-faces  and  cylindrical 
portions  of  the  body,  it  is  not  the  most  appropriate  location  for  conical  frustums.  For 
conical  frustums,  both  on  the  nose  and  t  ttail,  the  centroid  of  the  surface  area  is 


9  (x .  .  -  X  .Mr?  -  rh 

^  1+1  :  i+l  I 

3  (r.  -r.)^(r.^  +r) 

1+1  1  1+1  I 


(124) 


Here  x;  + 1,  r;  + 1  are  the  x  and  r  coordinates  of  the  aft  portion  of  the  frustum  and  xj,  r, 
the  X  and  r  coordinates  of  the  forward  portion.  Hence,  for  cylindrical  segments  and 
two-dimensional  sections  (such  as  wings),  the  pressure  computed  from  the  Shock 
Expansion  is  applied  at  the  midpoint,  whereas  for  conical  frustums,  it  is  applied  at 
the  point  given  by  Equation  (124). 
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2.3.5.2  Matching  Point  Between  MNT  and  SOSET 

Several  alternatives  have  been  used  in  past  literature  to  match  the  MNT  to 
SOSET  on  the  blunt  nose  part  of  the  body.  The  first  alternative  was  that  of  Jackson, 
et  al.20  where  the  body  slope  was  assumed  to  be  that  for  the  maximum  wedge  angle 
allowed  for  an  attached  shock  wave.  Their  results  showed  good  agreement  with 
experimental  data  for  2.3  ^  M®  but  only  fair  agreement  for  lower  supersonic  Mach 
numbers.  DeJamette,  et  ai.35  developed  an  empirical  equation  for  the  Mach  number 
to  match  MNT  to  SOSET.  This  approach  gave  improv^  results  over  Refe*-ence  20, 
particularly  at  low  supersonic  Mach  numbers.  Reference  2  foimd  that  using  a 
constant  value  of  Mj.  =  1.1  was  about  as  good  as  the  results  of  Reference  20. 

Appendix  B,  which  contains  a  derivation  c  the  new  methodology  by  DeJamette 
for  pressures  on  spheres  in  hypersonic  flows  found  a  constant  value  of  a  match  point 
of  Gfn  =  25.95°.  The  pressure  at  this  point  is  defined  by 


>•  M 


P,'P=  =  1 


-  C 

2 


(125) 


where  Cpj  is  given  by  Equation  (41).  Differentiating  Equation  (125)  there  is  obtained 


d  p  y  M*  p  ^  ^P,  d  6 

d  s  2  d  5  d  s 

eq 


But,  for  a  sphere 


(126) 


i/R. 


(127) 


Also,  from  Appendix  B 


d.s 

eq 


dC  j 


P  I  2 


i  1  c  +  — 

•  ^  Y  *  mn 


(128) 


C 

p 


=  C 


sin 


b 


(129) 


Using  the  relations  (126)  through  (129)  and  recognizing  that  5eq  =  25  95°  at  the 
match  point.  Equation  (116)  just  past  the  match  point  becomes 
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'  ds  /l 


A 


R. 


1 


P,V 


2 

1 


1 

1.05 M“p  (0.i676C^  -  1.251  C  )“  ! 


(130) 


Pi,  Vi,  Ml,  and  Cp^  of  Squaiion  (130)  are  computed  from  Equation  (79)  and  the  MNT 
of  section  2.3.4.2. 


Thus,  in  summary,  at  the  match  point  the  pressure  is  computed  from  Equation 
(41).  Downstream  of  this  point,  the  solution  is  continued  by  use  of  Equations  (35), 
(36),  (122),  ano  (130)  where  p^  is  the  value  of  pressure  computed  by  a  tangent  cone 
approach  (given  freestraam  Mach  number  and  local  slope,  given  entropy  along  the 
surface,  compute  the  value  of  p^.  by  the  approach  of  Appendix  C).  The  value  of  p, 
comes  from  the  Prandtl-Meyer  integration  in  Appendix  B. 

2.3.6  Angle-of-attack  Solution  for  B(«dv  Alone 

The  augie-of-attack  solution  will  follow  the  general  approach  of  References  2 
and  23  for  including  angle-of-attack  effects  in  the  pressure  distribution.  However, 
those  references  were  not  concerned  with  real  gas  effects  or  with  estimating 
tempei’atures  for  aeroheating  and  thermostructural  analyses.  Hence,  while  the 
general  approao>''<-  of  References  2  and  23  will  be  used,  modifications  and  additions 
will  be  added  ixe  the  approach  more  relevant  to  the  problems  at  hand. 

2.3.6. 1  hointed-cone  Pressure  Distribution  ata  >  0 

To  start  SOSET  on  a  pointed  body  or  to  compute  the  lc*cai  properties  on  a  conical 
frustum,  a  pointed-ccne  solution  for  a  >  0  is  needed.  One  would  desire  an 
approximate  solution,  if  good  accuracy  can  be  obtained,  as  opposed  to  a  numerical 
solution  of  partial  differential  equatiens^S  modified  to  include  real  gas  effects.  Tiie 
approach  of  Reference  23  for  perfect  gases ’S  given  by  Equation  (40).  Equation  (40a) 
is  independent  of  the  type  of  gas  and  is  therefore  appropriate  in  its  present  form  for 
computing  the  ACj,(  a,  0,  4>,  Ma,)  foi  a  real  gas  as  well.  However,  the  approximate 
relation  (40b),  for  the  pressure  coefficient  on  cones  at  zero  angle  of  attack,  is  replaced 
by  the  pressure  computed  in  Appendix  C  for  cones  in  real  gas  How.  ftudgins^? 
presents  an  approximate  formula  for  the  cone  surface  pressure  for  real  gas 
computations  based  on  a  fairly  broad  range  of  parameters  including  freestream 
velocity,  altitude,  and  cone  half-angles.  The  pressure  formula  is 


p 

«- 


1.4932{?4  1.3017 

*  c 


for  S3JI 8^  >  5 .5 


(131) 


or  since 
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/  • 


V  » 


^  .1 


•  • 

•c. 


* 

I  ♦  ' 


2  p 

-  =  - ^  ^ 

^0=0  y  M-  I  P, 


(C  )  =  -^—  i.4932  (  M  sine,.  -r  0.301 ' 

P  \  *  c ; 


(132) 


Reference  37  indicates  accuracies  of  better  than  1.5  percent  for  Equation  (131).  It 
should  be  pointed  out  that  for  values  of  Mo  sin  9c  ^  1.5,  real  gas  effects  on  pressure 
are  negligible  and  £k}uation  {<0b)  can  \e  used  instead  of  Equation  (132).  For  alues 
of  Ma>  sin  0c  >  1.5,  Equation  (132)  replaces  Equation  (40b). 

Reference  23  used  Equations  (40)  and  (40a)  to  estimate  angle-of-attack  effects 
on  the  cone.  A  slightly  more  accurate  method  will  he  used  here.  The  method  uses  the 
Equations  (40)  and  (40a)  but  combines  these  relations  in  a  Taylor  series  expansion  in 
angle  of  attack. 

Assume  small  a.ngles  of  attack  so  that  the  pressure  coefficient  can  be  written  in  a 
Taylor  Series  Expansion  about  a  =  0  as  foliov/s: 

I  1  I  J  j  3"C  j  q2 

C  lp,0ix,r),4>  =  |C  (6)j  +  -  n  o  - -  — 

PI  !  I  ?  I  do  ;i„2  a=0  2: 


.  ax  5  ; 

+  i - P  -  +  ... 

1  ac^  3^ 


(133) 


Now,  if  we  lise  the  methodology  for  the  cone  at  angle  of  attack  to  evaluate  the 
derivatives  of  Equation  (133),  then  referring  to  Equations  (40)  and  (40a) 

I  ] 

C  —  !C  i6  )!  ,  —  sin 2a  san20  costb 

=  {pc  jo-O  c  ^ 

-i-  sir:‘c.«-0^  M  2  -  -  ji  I  -  uin-0^  )  -  (  2  +  “  )  sir.*^4,|  (134) 

Oifferentiazing  Equation  (134)  and  evaluating  at  a=0,  one  obtains 


1  'S  I 

1  -  ‘  .  =  -  2  sn  (26)  cosA) 

1  ia  5"='' 


(135) 


NAVSWCTR  90-683 


Here,  the  value  of  6  is  taken  to  be  that  that  exists  on  the  conical  segment  in  question. 
In  essence,  use  of  Equations  (133)  and  (134)  allows  an  analogous  but  more  accurate 
technique  than  a  tangent  cone  approach  to  solve  for  the  pressure  coefBcient  on  a  body 
at  angle  of  attack  without  recomputing  local  cone  solutions  for  each  point  around  and 
along  the  body.  Differentiating  Equation  (134)  a  second  and  third  time,  and 
evaluating  at  a  =  0,  one  obtains 


where 


2F  cos”  6 


and 


(136) 


(137) 


Substituting  (135),  (136),  (147)  into  (133),  one  obtains 


C  (o,6{x,r;,4>)  =  IC  (0) 

P  1  P 


-  —  2a  sin  (26}cos4> 


+  (F  cus 


■0)  ^  — sin20cos4>j 


(138) 


Equation  (138)  indicates  that  the  local  angle-of-attack  effects  of  pressure  can  be 
estimated  with  the  zero  angle-of-attack  si.lution  already  computed,  the  local  value  of 
body  slope  and  roll  position,  and  the  freestream  Mach  number  and  angle  of  attack. 


The  other  question  is  how  dees  one  compute  other  flow  properties  at  the  cone 
surface?  Consider  ti.  o  cases.  First  is  the  case  where  the  cone  solution  is  being  used  as 
the  starting  solution  foi  the  pointed  body.  For  that  case,  since  entropy  behind  the 
shock  is  not  known,  the  flow  field  for  the  cone  at  a  =  0  must  be  computed  as  in 
Appendix  B.  Since  tho  flow  properties  are  ccr-  puted  numerically  to  the  accuracy  of 
the  thermodynamic  curve  fits  of  Reference  35,  the  pressure  on  the  cone  surface  is  also 
given  by  this  solution  i-ather  than  Equation  (132)  or  (40b).  Now,  since  Equation  (138) 
gives  the  Cp  due  to  a  for  various  values  of  6,  M,  then  a  new  Cp  for  4>  =  180“ 
(windward  plane)  is  computed  with  this  eqtmtion.  A  new  value  of  an  effective  cone 
angle  corresponding  to  the  windward  plane  is  then  given  by  solving  Equation  (132) 
for  sin  6  with  Co  replaced  by  the  value  of  p  corre^onding  to  C^  _  Then  the  cone 
solution  of  Appendix  B  is  repeated  with  a  new  cone  angle.  This  gives  the  value  of 
entropy,  temperatiaa,  and  other  properties  as  vreli.  Since  the  flow  between  the  body 
and  shock  is  isentropic,  tliis  value  of  entropy  is  held  constant  for  the  rest  of  the 
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computations  downstream.  Once  the  value  of  entropy  is  known  and  Cp  at  other 
positions  around  the  body  is  computed  from  Equation  (138),  then 


T(e,<j>,o)  =  T  j  p  (9,^,0  ,  S 


can  be  determined  from  the  thermofit  equations. 


(139) 


Other  flow  properties  can  be  computed  from  the  thermofit  equations  or  energy 
equations  as  follows: 


P  -  p 


p(0.4>,a)  .S 


a{9,4>,a)  =  a  I  p{6,4),a)  ,s|  ^ 


e(e,4»,a)  =  e|  ?(C,4».q) 


J 


h  =  e  +  p/p 


V  = 


2(H-h) 


M  =  V/a 


a^p 


r  = 


Z  =  p/{pBT.> 


I 

J 


(140) 


(141) 


The  aeccnd  case  »,f  interest  is  when  we  desire  to  compute  the  flow  properties  for 
a  cone  at  angle  of  attack  for  use  in  a  tangent  cone  approach  where  entropy  has  been 
defrned  by  the  properties  across  either  a  pointed  body  or  blunt  body.  In  both  of  these 
cases,  the  entropy  is  fixed  at  either  the  value  in  the  windward  plane  of  the  starting 
cone  solution  or  that  on  the  stagnation  streamline  behind  a  normal  shock.  For  tliis 
case,  since  entropy  is  known,  the  local  value  of  p  can  be  determined  directly  from 
Equations  (138);  (40),  t40a),  and  {40b)  if  Mco  sin  8c  ^  1-5.  Then,  the  values  of 
temperature,  density,  speed  of  sound,  and  internal  energy  can  be  computed  via 
Equations  (139)  and  (140)  as  previously  discussed. 
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2.3.6.2  Body  Angle-of-attack  Effects 

References  2  and  23  used  so-called  loading  functions  to  include  angle-of-attack 
effects  at  a  given  point  along  the  body  The  essence  ef  this  approach  was  Equations 
(40)  and  (40a).  A  similar  approach  will  be  used  here  except  the  slightly  improved 
pressure  prediction  method  given  by  Equation  (138)  ‘.viH  be  used.  In  essence, 
Elquation  ‘138;  can  be  used  in  a  tangent  cone  approach  to  compute  pressure  at  any 
point  on  the  body  surface  given  an  accurate  value  of  pressure  on  the  given  body  at 
zero  angle  of  attack.  That  is,  at  a  given  angle  of  attack,  local  body  slope,  position 
around  the  body,  and  freestream  Mach  number,  pressures  ail  around  the  body  can  be 
computed  from  Equation  (138).  Knowing  pressure  and  the  correct  value  of  entropy. 
Equations  (139)  through  (141)  can  be  used  to  determine  other  flow  properties  at  the 
body  surface. 

In  using  Equation  (138),  it  was  found  that  for  blunt  bodies  at  moderate  angles  of 
attack,  results  in  the  leeward  plane  v/ere  not  as  good  as  desired.  To  illustrate  tb^s 
point,  consider  Figure  14,  which  shows  the  perfect  gas  pressure  computations  on  a  20- 
percent  blunt  Von  Karman  ogive  at  M®  =  15  and  a  =  10“.  Figure  i4(a)  shows  the 
pressure  coefficients  predicted  by  Equation  (138)  in  several  planes  around  the  body, 
and  Figures  14(b)  through  (d)  compare  the  results  in  the  4)  =  180°,  90“,  and  0“  planes 
to  the  exact  Euler  code.  ZELS.41,4.  Jjj  examining  Figures  14(a)  through  (d),  two 
points  are  worth  noting:  (1)  the  present  method  agrees  reasonably  well  with  ZEUS 
computations  in  the  d)  =  180“  and  90“  planes,  and  (2)  while  Equation  (138)  works 
well  in  the  windward  plane  area  for  cones,  it  exhibits  unacceptable  behavior  in  the 
leeward  plane  for  blunted  ogives. 

The  unacceptable  behavior  exhibited  by  Equation  (138)  in  the  leeward  plane  of 
configu.'-ations  other  than  cones  comes  from  two  phenomena:  (1)  for  a  blunt  body,  the 
angle  0  is  fairly  large  in  the  nose  tip  region  causing  the  second  term  of  Equation  (138) 
to  be  larger  than  desired  >n  this  region,  and  (2)  on  ogives  as  0  goes  to  zero  near  the 
shoulder,  the  second  and  fourth  terms  of  Equation  (138)  go  to  zero  whereas  the  third 
term  reaches  its  maximum  positive  value.  This  causes  the  increase  in  pressure  past  x 
=  1  seen  in  Figures  14(a)  and  i4(d  <.  It  was  found  that  a  way  to  remedy  both  of  these 
problems  in  the  4>  ^  30°  planes  was  to  modify  Equation  (138)  as  follows: 

2a  sin  i20)cos(}> 

- ; -  (142) 


C  J  a.  0,  <|>  1  =  C  iQ' 

P  p 


Equation  (142)  is  used  for  configurations  other  than  sharp  cones  in  the  leeward  plane 
area.  The  results  of  this  modification  are  shown  in  Figure  15(a)  for  pressure  and 
Figure  10(b)  for  temperature.  Note  that  significant  improvement  in  pressure 
coefficient  compared  to  the  ZEUS  code  is  obtained.  The  ZEUS  code  data  in  the 
various  planes  are  not  showm  on  Figure  15  for  clarification  piu-poses.  To  show  the 
improved  comparison,  one  can  simply  use  the  leeward  plane  curve  (4)  =  0°)  of  Figure 
15  and  compare  this  to  Figure  14(d).  Since  Equation  (138)  is  still  used  in  the  4)  = 
180°  and  4)  =  90°  planes.  Figures  14(b)  and  (c)  remain  the  same.  Also,  use  of 
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AXiAL  LOCAuON.  X  (CAL) 

FIGURE  14<b).  COMPARISON  OF  APPROXBUTE  ^EQUATION  138)  AI^D  EXACT  PRESSURES 
ONA20-PERCENTBLUNTVONKARMANOGIVE  {M^  =  15,o  =  10',  4.  =  ISC') 
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FIGURE  14(c).  COMPARISON  OF  APPROXIMATE  (EQUATION  13S)  AND  EXACT  PRESSURES 
ON  A  20-PERCENT  BLUNT  VON  KARMAN  OGIVE  (M=  =  15,  o  =  10“,  &  =  90“) 


AXIAL  LOCAUON,  X  (CAL) 

FIGURE  14(d).  COMPARKON  OF  APPROXIMATE  (EQUATION  138)  AND  EXACT  PRESSURE 
ON  A  20-PERCENT  BLUNT  VON  K.\RMAN  OGIVE  (M.  =  15.  a  =  10*,  <>  =  O') 


COEFFICIENT  OF  PRESSURE 
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AXIAL  LOCATION.  X  (CAL) 

FIGURE  15fb).  APPROXIMATE  TEMPERATURE  PP^DICTIOX  ON  A  20-PERCENT  BLU-NT 

VONKARMANOGn'E  =  15.  a  =  10*^) 
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Equation  (142)  versus  (138)  in  the  leeward  plane  area  presents  little  change  to  the 
axial  force  but  improves  the  normal  force  and  center-of-pressure  predictions  by  about 
5  and  6  percent,  respectively.  Compared  to  ZEUS  computations,  the  errors  in  the 
aerodynamic  coefficients  for  wave  drag,  normal  force,  and  center  of  pressure  using 
the  present  approach  of  Equations  (138)  and  (142)  is  5,  11,  and  4  percent, 
respectively.  This  is  considered  to  be  quite  acceptable  for  an  engineering  type 
aeroprediction  code. 

Past  the  nose  portion  of  the  body,  standard  Modified  Newtonian  Theory  is  used 
to  predict  the  pressuies  with  the  pressure  coefficient  assumed  to  be  zero  in  the 
shadowed  region.  That  is 

C  =  C  sin* <5  )  for  5  SO 

?  Pg  «! 

C  =  0  for  6  <  0  (143) 

P 

Since  entropy  is  constant  on  the  body  surface,  values  of  temperature  and  other  local 
properties  are  computed  with  Equations  (139)  through  (141)  as  before. 


2.3.7  Wing  and  Interference  Aerodynamics 

For  high  Mach  numbers  (Ma  typically  greater  than  6),  the  method  of  References 
2  and  21  will  be  used  for  real  gas  computations  except  that  the  pressures  and  other 
properties  will  be  computed  using  real  gas  as  opposed  to  perfect  gas  computations. 
References  2  and  21  used  local  slope  techniques  for  the  wing  aerodynamic 
computations.  That  is,  MNT  on  the  leading  edge  of  the  wing,  tangent  wedge  from  the 
match  point  rearward  on  a  compression  surface,  and  PME  on  expansion  surfaces.  The 
methods  of  section  2.3.4,  Appendix  C,  and  .Appendix  D,  respectively,  are  used  to 
compute  the  local  pressure  coefficient  and  temperature  at  the  outer  edge  of  the 
boundary  layer.  The  match  point  between  MNT  and  the  tangent  wedge  is  taken  to  be 
0M  =  25.95“.  Since  this  is  a  first-order  as  opposed  to  a  second-order  approach,  there 
could  be  a  slight  discontinuity  in  the  pressure  coefficient.  This  is  due  io  the  fact  that 
there  is  no  pressure  derivative  to  blend  the  pressure  coefficients  together  as  with 
SOSET. 

Use  of  a  first-order  as  opposed  to  a  second-order  solution  for  the  w'ings  and  other 
lifting  surfaces  can  be  justified  by  the  fact  that  these  surfaces  generally  contribute 
less  than  25  percent  of  the  axial  force.  Experience^O  has  shown  that  the  axial  force  is 
where  second-order  accuracy  is  required  for  pressure  coefficients.  On  the  other  hand, 
lift  can  be  reasonably  well  estimated  by  a  first-order  theory.  Since  the  win^ 
contribute  significantly  to  the  lift  force,  a  first-order  theory'  is  acceptable. 

It  was  found  in  References  2  and  21  that  when  local  slope  techniques  were  used 
in  a  2-D  “strip  theory”  approach  (wing  is  considered  to  be  made  up  of  several 
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streamwise  2-D  surfaces  and  total  forces  and  moments  summed  from  these  surfaces), 
no  interference  aerodynamics  were  needed.  This  was  due  to  the  fortuitous 
cancellation  of  errors  and  the  fact  that  2-D  pressures  were  higher  than  3-D  pressures 
on  most  wings.  At  the  wing  root,  the  pressure  was  underpredicted  due  to  use  of 
freestream  versus  local  conditions,  and  at  the  wing  tip  the  pressure  was 
overpredicted  due  to  the  same  reason.  The  two  effects  canceled  each  other  in  an 
approximate  sense.  A  summary  of  the  approaches  for  computing  real  gas  properties 
on  bodies  and  wing  bodies  using  SOSET  is  given  in  Figure  16. 


I.  BOOYATa  =  Q 

1.  For  pointed  body,  use  cone  solution  for  real  gases  of  Appendix  C  to  start  solution.  For 
blunt  bodies,  use  MNT  of  section  23.4  to  start  solution. 

2.  Use  Equations  (35),  (36),  and  (122)  to  calculate  pressure  on  conical  segments  along  body 
downstream  of  tip. 

3.  Use  Equations  (133).  (140),  {l4l)to  calculate  otherflow  properties  on  each  conical 
segment  given  the  entropy  from  step  1  and  pressure  from  step  2. 

ii.  POINTED  BODY  ATa  >  0 

1 .  Computa  (Cp)a from  Equation  (1 38}  and  (Cp)a = o  f rom  Equation  (1 32)  (or  Equation  (40b} 
if  Mx  sin  St  £  1 .5).  This  gives  Cp  on  cone  at  a  >  0. 

2.  Compute  effective  Gtfrom  Equation  (131)  where  pt  of  Equation  (131)  is  ••eplaced  by 
thatfor*})  =  180®  of  step  1. 

3.  Use  this  value  of  6  to  obtain  the  cone  solution  for  real  gases  and  obtain  the  value  of 
entropy  between  the  shock  and  body. 

4.  Compute  T,  a.  e.  h.  V.  y,  Zfrom  Equations  (1 33),  (1 54).  and  (135)  around  conicai  tip. 

5.  Use  Equations  (138),  (1 42).  (143),  and  Appendix  D  to  compute  pressures  downstream 
on  body  at  a  >  0. 

6.  Using  Cp  (a.6,4>)  from  step  5,  tiie  value  of  entropy  from  step  3.  compute  otiier  flow 
properties  from  Equations  (139),  (140),  and  (141 }  at  desired  points  along  and  around 
body. 


FIGURE  16.  SUMM.4RY  OF  SOSET  FOR  RE.\L  G.4SES  ON  BODIES 
.4ND  WLNG  BODY  CO.MBIN.ATIONS 
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HI.  BLUNT  BODY  ATa  >  0 

1 .  Compute  conditions  ocoss  noimaf  shock  and  at  stagnation  point  from  Appendix  C  and 
section  23.4. 

2.  Compute  conditions  up  to  match  point  from  section  23.4. 

3.  Compute  body  solution  at  a  =  9  as  described  above 

4.  Compute  C^{olO,4>}  from  equations  (138).  (142).  and  (143). 

5.  Using  the  value  of  entropy  from  step  1.  pressure  from  step  4.  compute  other  properties 
from  Equations  (139).  (140).  and  (141). 


iV.  WING 

1.  For  blunt  leading  edge,  computeflow  properties  at  each  point  based  on  section  23.4. 
For  sharp  leading  edges,  use  wedge  solution  for  real  gases  of  Appendix  C  Angies  used 
are  local  angles  that  include  angle  of  attack,  local  slope,  and  sweepback. 

2.  Use  PME  and  tangent  wedge  aft  of  the  leading  edge  usin^  local  slopes  and  freestream 
conditions.  This  gives  all  flow  properties  directly  based  on  the  entropy  of  step  1  and  the 
pressure  of  PME  and  tangent  wedge. 


HGURE 16.  SUMMARY  OF  SOSET  FOR  REAL  GASES  ON  BODIES 
AND  WING  BODY  COMBIN.ATONS  (CONTINUED) 


If  one  is  most  interested  in  correct  values  of  wing  lift  and  axial  force  and  not  as 
concerned  with  high  temperature  effects,  the  3-D  thin  wing  theory  of  Reference  31 
combined  v/ith  interference  factors  is  more  accurate.  This  approach  is  still  allowed  as 
an  option  to  the  user  in  the  aeroprediction  code.2-4 


2.3.8  Viscous  Considerations 

The  main  intent  of  this  report  was  to  extend  SOSET  to  include  real  gas  effects 
and  hence  allow  computation  of  forces,  moments,  and  surface  temperatures  at  Mach 
numbers  above  about  6.  With  the  exception  of  surface  temperature,  these  objectives 
have  now  been  accomplished.  However,  some  additional  comments  are  in  order  as  to 
how  to  include  the  viscous  effects  that  occur  on  the  wing  and  body.  Ultimately,  the 
current  technique  will  be  transitioned  into  the  latest  version  of  the  NAVSWC 
aeroprediction  code2-4  and  an  updated  version  will  be  generated. 
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The  skin  friction  and  body  base  pressure  drag  can  be  estimated  by  the  same 
techniques  of  Reference  30,  Van  Driest  11,39  and  empirical  methods,  respectively. 
The  2-D  base  pressure  drag  on  the  rear  of  wing  surfaces  is  again  empirical  and  the 
same  approach  as  that  of  Reference  31  can  be  used.  The  nonlinear  effects  of  lift  and 
drag  that  occur  due  to  higher  angles  of  attack  can  be  accounted  for  by  the  Allen- 
Perkins  Viscous  Crossflow  Theory .40 

The  other  viscous  consideration  is  that  of  temperature.  Recall  that  one  of  the 
main  objectives  of  the  present  work  was  to  calculate  real  gas  temperatures  at  the  wall 
for  use  in  heat  transfer  analysis.  The  temperatures  calculated  by  the  present 
analysis  are  those  that  exist  at  the  outer  edge  of  the  boundary  layer  and  thus  must  be 
carried  to  the  wall  and  heat  transfer  analysis  performed  to  get  the  true  wall 
temperature.  This  can  be  done  with  one  of  several  stata-of-the-art44  engineering 
techniques.  However,  without  the  more  accurate  real  gas  versus  perfect  gas 
temperatures  computed  at  the  outer  edge  of  the  boundary  layer,  these  estimates  of 
W’all  temperature  will  be  considerably  in  error  at  higher  Mach  numbers.  This  in  turn 
could  cause  costly  and  possibly  unnecessary  design  decisions.  Hence,  while  the 
present  report  is  primarily  focused  on  providing  improved  estimates  of  temperatures 
at  the  outer  edge  of  the  boundary  layer,  those  temperatures  will  make  local  wall 
temperature  and  heat  transfer  computations  (which  are  dependent  on  material 
properties,  insulation  assumptions,  etc.)  much  more  accurate. 

The  present  inviscid  analysis  gives  values  of  all  the  properties  at  the  outer  edge 
of  the  boundary  layer.  They  could  thus  be  input  to  a  boundary  layer  code  and  surface 
properties  calculated  for  heat  transfer  analysis.  To  aid  in  heat  transfer  analysis,  an 
approximate  technique  will  be  given  for  estimating  two  of  the  parameters  needed  as 
inputs;  the  adiabatic  wall  temperature  or  enthalpy  and  heat  transfer  coefficients. 
While  the  technique  is  not  as  elaborate  as  a  numerical  solution  of  the  axisymmetric 
or  three-dimensional  boundary  layer  equations,  it  does  give  an  engineering 
approximation.  The  results  can  then  be  compared  to  more  accurate  computations. 

The  aerodynamic  heating  to  the  wall  is  defined  by 

-h,)  (144) 

w 

Here  the  subscript  l  refers  to  inviscid  computations  at  the  outer  edge  of  the  boundary 
layer  so  these  quantities  are  known.  Thus,  by  computing  the  heat  transfer  coefficient 
Ch  and  adiabatic  wall  enthalpy,  the  heat  transfer  to  the  wall  for  a  given  wall 
enthalpy  can  be  computed. 

The  adiabatic  wall  enthalpy  means  the  change  in  enthalpy  normal  to  the  wall  is  zero. 
That  is 


(145) 
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where  y  is  normal  to  the  wail.  Solution  of  the  boundary  layer  equations  with  this 
boundary  condition  will  give  the  true  adiabatic  wall  enthalpy.  However,  most 
engineering  computations  are  performed  with  ha.^  related  to  the  recovery  factor  Tc 
as44 


=  h,  +  r 
L  c 


2 


For  a  turbulent  boimdary  layer,  r  has  been  shown  to  be 


(146) 


For  a  real  gas,  Prandtl  number  is  a  function  of  other  thermodynamic  properties  and 
must  be  foimd  from  the  thermodjmamic  fit  equations  of  Keferences  34  and  35.  That  is 


Pr  = 


(147) 


Thus,  since  Pj^,  Tj^,  hj^,  are  all  known  fii-om  the  inviscid  computations,  ha.^  can  be 
computed  from  Equation  (146).  The  adiabatic  wall  temperature  is  also  defined  in 
many  references  by 


2 

1  +  r  — -  .Mf 

C  2  ^ 


(148) 


The  perfect  gas  assumption  can  be  somewhat  compensated  for  in  Ekauation  (152)  by 
using  the  local  value  of  as  opposed  to  1.4. 

To  compute  the  heat  transfer  coefficient,  resort  will  be  made  to  the  reference- 
temperature  method45,46  and  Reynolds  analogy .33  The  reference  temperature,  T*,  is 
given  byll3 


T 

T 


L 


1  +  0.032  + 


Then  the  local  skin  friction  coefficient  is33 


0.0592 


0.2 


(149) 


(150) 


Here 
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.  » 

/• 


^  1 


,  C  • 


*\ 


^-  = 


p  V 


(151) 


p*  and  Uj*  are  evaluated  at  the  local  pressure  and  reference  temperature  defined  by 
Equation  (149)  fi’om  the  real  gas  thermofit  eaaations.S4,35  That  is 


p  •=  D  (T  .p.) 


p,  =  Pj  (T  ,p  } 


(152) 


(153) 


Reynolds  analogy  then  gives 


-  — (p  )  ^ 

2  ^ 


(154) 


Pr*  is  also  evaluated  from  the  thermofit  equations  where 


P.  -  ,p  ) 


(155) 


Finally,  if  one  desires  enthalpy  as  opposed  to,  or  in  addition  to,  temperature,  then 


h  =  h  j 


h  =  h  (pj  ,  T  ) 


{156a) 


(156b) 


can  be  computed  once  again  from  the  thermofit  equations  since  pressure  is  assumed 
ocnstant  across  the  boundary  layer.  That  is 


Equations  (150)  and  (154)  are  valid  for  flat  plates  in  a  turbulent  compressible 
boundarj'  layer.  They  are  tlius  the  appropriate  equations  for  the  wings  and  other 
lifting  surfaces  on  the  configuration  of  interest.  For  the  body,  Reference  33  indicates 
that  Elquations  (150)  and  (154)  should  be  multiplied  by  the  Mangier  fraction,  V  3. 
Hence 


^  _  0.103 


(157) 
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(158) 


Again  it  should  be  emphasized  that  the  reference  enthalpy  or  temperature  and 
Reynolds  analogy  are  approximate  engineering  techniques.  They  have  been  used 
extensively  and  have  a  theoretical  foundation.  The  more  rigorous  alternative  of 
using  the  inviscid  properties  as  inputs  to  a  chemically  reacting  turbulent  boundary 
layer  computation  is  considered  beyond  the  scope  of  the  present  work. 


3.  RESULTS  AND  DISCUSSION 


3.1  STARTING  SOLUTIONS 

It  is  of  interest  first  tn  check  the  validity  and  accuracy  uf  the  starting  solutions 
for  both  the  pointed  and  blunted  body  cases.  Figures  17  and  18  compare  ♦he  perfect 
gas  pressure  and  temperature  predicted  by  the  present  methed  to  the  exact  cone 
solution  of  JonesS3  fn»  several  conditions.  The  first  of  these  conditions,  M»  =  10.6, 6c 
=  15°,  a  =  10°,  is  the  only  one  of  interest  as  far  as  real  gas  effects  are  concerned. 
However,  the  other  lower  Mach  number  and  higher  angle-of-attack  cases  are  shown 
to  indicate  the  robustness  of  the  present  method  for  computing  pressure  and 
temperature.  Note  that  for  ah  four  condition?  of  Figures  17  and  18,  good  agreement 
un  pressure  and  temperature  is  obtained  compared  to  the  exact  solution,  in  most 
cases,  accuraci*;3  of  3  percent  or  better  are  found  for  the  conditions  investigated  to 
date.  An  exception  is  'Jie  pressure  predicrion  at  the  bottom  of  Figiu-e  18  where  errors 
in  pressure  of  up  to  1 5  percent  are  encountered.  Comparing  Figures  17  and  18  shows 
that  the  method  gives  slightly  iniproved  accuracy  at  the  higher  Mach  numbers. 

Figure  16  compares  the  present  conical  starting  solution  with  the  exact  results 
from  the  ZEUo  tode-*^  for  M®  =  15,  a  =  10°,  6c  =  Roth  real  and  perfect  gas 
results  are  shown.  Note  that  tlie  approximate  and  ZEUS,  real  ard  perfect  gas 
pressures  are  nearly  identical  (Figure  i9(a».  This  is  as  anticipated  frora  knowledge 
of  compressioii  processes  with  respect  to  real  gas  effects.  On  the  other  hand.  Figure 
19(b)  shows  the  lar,je  differences  in  temperature  that  can  occur  due  to  real  gas  effects, 
lo  comparing  approximate  and  exact  temperature  calculations  in  Figure  19(b),  it  is 
seen  that  the  perfect  gas  results  are  within  a  couple  cf  percent  whereas  the  real  gas 
deviates  by  up  to  10  percent.  The  ZEUS  real  gas  temperatures  are  actually  higher 
than  the  perfect  gas  results  in  the  leeward  plane,  which  is  not  correct.  The  reason  for 
this  is  that  the  ZEUS  code  uses  the  real  gas  thermofit  equations  to  obtain  entropy 
directly.  On  tl'a  other  hand,  referring  to  Figure  C-15  and  the  discussion  that  goes 
w  itf-  the  figure,  it  is  seen  that  obtaining  entropy  in  this  manner  gives  temperatures 
about  250°  too  high  for  cone  aagies  cf  15°  and  25°.  Referring  to  Figure  19(b),  this  is 
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the  approximate  value  of  the  difference  betwet*;  the  ZEUS  code  and  the  present 
computations.  This  250“  reduction  would  apply  fjl  aroimd  the  body  since  this 
incorrect  value  of  entropy  is  held  constant.  This  explains  why  there  is  a  crossover  of 
the  real  and  perfect  grs  temperatures  with  the  ZEUS  computations  in  Figure  19(b). 
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EXACT  COME  SOLN  (REF  38) 
APPilOXiMATE  METHOD 


•RfiSSURl 
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Figures  20(a)  through  (d)  illustrate  the  accuracy  of  the  Improved  Modified 
Newtonian  Theory  (IMNT)  over  MNT.  The  figures  show  the  IMNT,  MNT,  and  the 
exact  results  from  Reference  48.  Note  that  the  IMNT  gives  virtually  identical  results 
to  the  exact  computations.  Although  not  shown,  this  excellent  agreement  in  pressure 
prediction  holds  true  all  the  way  to  M®  =  30.  In  the  present  method,  the  match  point 
between  SOSET  and  IMNT  is  at  8eq  =  25.95°  so  that  IMNT  is  used  for  8eq  ^  25.95° 
and  SOSET  for  values  of  8eq  <  25.95°.  Knowing  accurate  values  of  pressure  and  a 
good  estimate  of  entropy  should  allow  good  estimates  of  other  properties  at  the  body 
surface. 


90  80  70  60  50  40  30  20  10  0 
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FIGURE  20ia).  SURF.4CE  PRESSURE  DISTRIBUTION  0\T:R 
A  HEMISPHERICAL  FOREBODY  .AT  M*  =  5 
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FIGURE  20fb).  SURFACE  PRESSURE  DISTRIBUTION  0\^ 

A  HEMISPHERICAL  FOREBODY  AT  M«  =  10 
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HGURE  20(c).  SURFACE  PRESSURE  DISTRIBUTION  OVER 
A  HEMISPHERICAL  FOREBODY  AT  M*  =  15 
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3.2  PRESSURE  GRADIENT  USED  IN  SOSET 

Before  computing  pressures  and  temperatimes  on  general  body  shapes  using  the 
section  3.1  starting  solutions  combined  with  SOSET,  it  is  of  interest  to  examine 
Equation  (36),  which  defines  q.  Recall  that  q  is  the  exponent  containing 


2 


which  is  used  in  Equation  (35)  in  the  SOSET  to  determine  the  blend  of  Tangent  Cone 
Theory  (TCT)  and  Generalized  Shock  Expansion  Theory  (GSET).  To  do  this,  a  simple 
biconic  nose  shape  was  considered  as  indicated  in  Figure  21.  This  allowed  a  simple 
starting  solution  and  then  a  computation  of  the  properties  around  a  turn  of  a  given 
A8  =  02  -  01  as  shown  in  Figure  1.  Results  were  obtained  for  the  parameter 


imder  various  conditions  using  the  approximate  method  of  Syvertson  and  Dennis, 19 
the  exact  prediction  by  DeJamette,23  and  the  approximate  present  technique  given 
by  Equation  (122).  Since  the  Reference  19  and  23  results  were  for  perfect  gas  only, 
the  cases  considered  for  comparison  were  for  perfect  gas.  The  results  of  this  study  are 
presented  in  Table  1,  which  shows  M®  varying  from  3  to  15,  cone  half  angle  from  5  to 
20,  and  A0  2“  and  5°.  Note  first  of  all  that  the  present  approximation  for 


agrees  very  closely  with  the  exact  integration  of  Reference  23.  The  only  difference 
between  the  two  techniques  is  the  fact  that  the  present  method  approximates  the 
integral  of  Equation  (112)  iising  the  Trapezoidal  Rule  whereas  Reference  23 
integrates  the  entire  equation  numerically.  On  the  other  hand,  the  Reference  19 
results  neglect  the  integral  of  Equation  (112)  altogether.  As  seen  in  Table  1, 
neglecting  the  integral  gives  reasonable  results  (within  10  percent)  in  most  cases 
shown  where  A0  is  small. 


The  most  important  aspect  of  the  Table  1  results  is  the  fact  that  q  becomes 
negative  quite  often  for  Mach  numbers  of  5  and  greater.  This  means  that  q  must  be 
defined  and  SOSEIT  of  Equation  (35)  reverts  back  to  GSET  or  TCT;  that  is,  p  =  p2  or 
p  =  Pc.  As  a  result  of  this  dilemma,  we  investigated  a  different  approach  that  would 
allow  us  to  take  advantage  of  both  pc  and  P2.  This  approach  defined  p  in  the  SOSET 
3S 


p  =  p^-fp^-p2)nj 


(159) 
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=  CONE  HiOf  ANGLE  =  PRESSURE  BEHIND  AN  EXPANSION 

a,  =  SHOCK  ANGLE  P^  =  PRESSURE  FAR  DOWNSTREAM  ON  A  ERUSTRUM 

a9  =  e^-  82 


FIGURE  21.  BICONIC  CONFIGURATION  GEOMETRY  FOR  INVESTIGATING 
PRESSURE  GRADIENT  AROU'ND  A  CORNER  (see  Table  1 ) 


TABLE  1.  ESTIMATES  OF  (ra  ./p.)  ( 8p  /  as )  2  FOR  VARIOUS  CONDITIONS 


M. 

Oc 

(deg) 

AO 

Exact  23 

Syv& 
Dennis  19 

Present 

(ra/p*)  {dp  Ids)  2 1 
(Pc-P2) 

P2 

Pc 

3 

15 

2 

0.0216 

0.0209 

0.0216 

0.657 

1.86 

1.89 

5 

0.0478 

0.0441 

0.0478 

0.622 

1.50 

1.58 

5 

20 

2 

-0.0148 

-0.0142 

-0.0147 

-0.143 

4.68 

5 

-0.0129 

-0.0109 

-0.0120 

-0.077 

3.56 

mm 

5 

10 

2 

0.0153 

0.0147 

0.0154 

0.389 

1.88 

1.92 

5 

0.0341 

0.0306 

0.0341 

0.359 

1.34 

1.43 

10 

15 

-0.0879 

-0.0822 

-0.0872 

-0.209 

8.29 

8.71 

5 

-0.0952 

-0.0785 

-0.0896 

-0.188 

5.20 

5.71 

10 

10 

2 

-0.0144 

-0.0131 

-0.0141 

-0.097 

3.96 

4.11 

5 

0.0014 

0.0030 

0.0034 

0.010 

2.20 

2.35 

15 

15 

2 

-0.2602 

-0.2580 

-0.206 

16.73 

5 

-0.2840 

-0.2672 

-0.193 

Ba 

15 

10 

2 

-0.0719 

-0.0649 

-0.0706 

-0.141 

7.17 

7.68 

5 

-0.0420 

-0.0294 

-0.0351 

-0.122 

3  43 

3.77 

15 

5 

2 

0  0061 

0.0056 

0.0063 

0.068 

2.04 

2.13 

5 

0.0165 

0.0135 

0.0167 

0.051 

0.72 
_ ! 
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where  iij  =  0  gives  p  =  Pc  and  ili  =  1  gives  p  =  P2  and  0  <  <  1  gives  a  blend  of  the 
two.  Thds  approach  would  allow  a  direct  input  of  q,  as  opposed  to  a  computation  as 
done  in  the  past. 

To  investigate  the  Equation  (159)  approach,  two  cases  were  considered.  The 
first  of  these  is  a  20-percent  blunt  Von  Karman  ogive  at  Moo  =  15.  Results  of  the 
pressure  computation  are  shown  in  Figure  22.  In  Figure  22,  IMNT  +  P2  is  when 
=  1,  IMNT  +  Pc  is  where  qj  =  0,  old  AP  is  the  current  aeroprediction  code,  and  ZEUS 
is  the  exact  computation.  In  the  current  aeroprediction  code,  SOSET  is  implemented 
with  Equation  (35)  and  when  q  becomes  negative,  p  is  set  to  Pc.  As  a  result,  note  in 
Figure  20  that  the  present  results  for  q^  =  1  are  slightly  better  than  the  old  AP 
compaxed  to  the  ZEUS  computations.  This  is  due  solely  to  the  IMNT.  However,  if  one 
were  to  use  q^  =  1,  significantly  worse  results  would  be  obtained. 


-X- 

IMNT  +  P2 
IMNT  +  PC 
OLD  AP 


ZEUS 


FIGURE  22.  COMPARISO.V  OF  V.4RIOUS  PRESSURE  PREDICTION  TECHNIQUES 
ON  A  20-PERCENT  BLU.NT  VON  KAR.MAN  OGIVE  .M.  =  15.  a  =  0”,  h  =  100k  ft; 
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The  second  case  investigated  to  examine  the  Equation  (159)  approach  is  a 
simple  hemisphere  forebody.  These  results  are  shown  in  Figure  23.  Note  here  that 
the  IMNT  +  Pc  and  old  AP  results  are  inferior  to  the  IMNT  +  p2  results  compared  to 
ZEUS.  This  is  even  more  reason  for  allowing  to  be  input  as  opposed  to  computed. 

To  summarize  the  pressure  prediction  using  SOSET,  it  was  found  that  for  Mach 
numbers  greater  than  5,  q  becomes  negative  quite  often  requiring  one  to  revert  back 
to  either  TCT  or  GSET.  Since  this  is  the  case,  it  is  believed  that  the  Equation  (159) 
approach  is  a  better  way  to  implement  SOSET.  A  value  of  q^  =  0  is  recommended  for 
hemispheres  or  near  hemispheres.  For  other  configurations,  a  value  of  q,  =  1.0  is 
recommended.  This  is  the  approach  that  is  being  implemented  in  the  aeroprediction 
code  for  Mach  numbers  greater  than  about  6.  For  Mach  numbers  between  about  2.5 
and  6,  conventional  SOSET  is  used,  and  for  Mach  numbers  1.2  to  about  2.5,  second- 
order  Van  Dyke  Theory 30  is  used. 


3.3  CONFIGURATION  RESULTS 

To  thoroughly  investigate  the  current  new  methods  for  calculating  real  gas 
effects,  several  configurations  were  selected  that  covered  the  typical  range  of  missile- 
type  configurations.  These  include  a  20-percent  blunt  Von  Karman  ogive,  a  20- 
percent  blunt  Von  Karman  ogive-cylind'^r,  a  blunt  cone  nose  shape,  and  a  sharp  cone- 
cylinder-flare  configuration  (see  Figure  24).  Exact  solutions  were  also  generated  as  a 
basis  for  comparison  of  the  present  approximate  engineering  methods. 

Figure  15  showed  the  results  of  using  the  IMNT  +  SOSET  for  the  zero  angle-of- 
attack  solution,  applying  Equations  (138)  and  (142)  to  get  angle-of-attack  pressure 
effects,  and  knowing  entropy,  applying  Equation  (139)  to  get  temperature.  Figure  16 
was  for  perfect  gas.  Figure  25  further  illustrates  the  present  methodology  for  the  Von 
Karman  ogive  by  comparing  the  pressure  and  temperature  in  the  windward  plane 
(perfect  gas)  to  conventional  TCT  and  ZEUS  results.  Conventional  TCT  is  where  one 
computes  pressure  at  a  given  point  on  a  body  surface  with  the  freestream  Mach 
number  and  the  local  slope  with  respect  to  the  freestream  velocity  vector.  Tne 
standard  Taylor  McCaiUS  cone  solution  can  be  solved  for  perfect  gases,  Hudgins3'' 
solution  for  real  gases  or  approximations  to  these  solutions  such  as  Equations  (40b)  or 
(132).  The  IMNT  was  used  on  the  blunt  nose  portion.  Note  that  in  comparing  the 
three  theoretical  approaches  for  both  pressure  and  temperature,  the  present  method 
agrees  closer  with  the  exact  solution  than  does  the  TCT. 

Figure  26  presents  the  results  for  the  same  case  as  Figure  25  except  here  the  gas 
is  equilibrium  chemically  reacting.  Note  that  in  comparing  Figure  26(a)  to  25(a), 
little  difference  in  pressxire  is  noted  (and  indeed  in  the  aerodynamic  coefficients); 
however,  comparing  Figure  26(b)  to  Figure  25(b),  real  gas  temperatures  are  lower  by 
as  much  as  a  factor  of  2.  Figure  26(c)  compares  the  present  predictions  to  those  of  the 
ZEUS  code  in  three  planes.  Agreement  on  temperature  in  the  critical  windward 
plane  is  within  4  percent. 
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FIGURE  23.  COMPARISON  OF  VARIOUS  PRESSURE  PREDICTION  t'ECHNIQUES 
ON  A  HEMISPHERE  FOREBODY  (M*  =  15.  o  =  O',  h  =  100k  ft) 
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FIv.-URE  24.  CONFIGURATIONS  USED  AS  TEST  CASES  FOR  NEW  THEORY 
fDDklENSIONS  IN  CALIBERS' 
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TEST  CASE  4:  20-PERCEMT  BLUNT  VON  KARMAN 
CGJVE— CYUNDER  WITH  CRUCIFORM  TAILS 


EV;  THEORY 


FIGURE  24.  CONFIGURATIONS  USED  AS  TEST  CASES  FOR  N 
< DIMENSIONS ’NCALiBEPS?  ‘CONTINUED; 
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FIGURE  2e(a).  APPROXIMATE  PRESSURE  PREDiCTIOK  ON  A  20-PERCENT  BLUNT  VON 
KARXLAN  OGi\=r:  FOR  A  REAL  GAS  iM*  =  15,  o  =  10®,  h  =  100k  fi> 
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FIGURE  26>c!.  C0MP.4RIS0N  OF  .4PPR0XIM.ATE  .4ND  EX.\CT  TEMPER.4TURE  PREDICTION 
ON  A  2&-PERCEN7  BLUNT  VON  K.\R.M.AN  OGI\T  FOR  .A  RE.AL  G.AS 
<M=  =  15,0  =  ir,h  =  100k  fu 
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The  second  test  case  considered  is  the  blunt  cone  shown  in  Figure  24.  While  the 
bluntness  is  the  same  as  the  Von  Karman  ogive,  the  aft  part  of  the  nose  has  a 
constant  angle  so  it  is  of  interest  to  see  how  the  results  of  the  approximate  method 
compare  w'ith  ZEUS  computations.  These  comparisons  are  presented  in  Figure  27. 
Figures  27{a)  and  27(b)  present  the  windward  plane  pressure  coefficient  and 
temperature  conparisons  between  the  ZEUS  and  approximate  computations.  In 
examining  Figure  27(a),  ZEUS  computations  do  show  some  slight  differences  in 
pressure  coefficient  between  real  and  perfect  gas  on  the  overexpansion  region  of  the 
blunt  nose  whereas  the  approximate  results  give  basically  identical  pressures.  Also 
note  that  good  agreement  between  the  exact  and  approximate  technique  is  evident 
except  in  the  overexpansion  region.  In  examining  Figure  27(b),  good  agreement  in 
temperature  predictions  between  the  tw’o  theoretical  methods  is  noted.  Even  in  the 
overexpansion  region  where  the  disagreement  is  largest,  the  maximum  deviation 
from  the  exact  results  is  only  about  8  percent.  On  the  maximum  temperature 
portions  of  the  nose,  the  deviation  is  less  than  2  percent.  Also  worthy  of  note  is  the 
almost  5000“  R  lower  temperature  of  real  versus  perfect  gas  computations. 


!  •  i 

0  J _ ! 
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FIGURE  27<.t  I  CONfP.ARISON  OF  .APPROX LM.ATE  .A.VD  EX.ACT  PRESSURE  COEFFICIE.VTS 
IN  THE  WL\*D\VARDPL.A.\'E  OF  .A  20-PERCENT  BLUNT  CO.VE  (NU  =  15,o  =  10*) 
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FIGURE  27(b).  COMPARISON  OF  APPROXIMATE  AND  EXACT  TEMPERATURE  IN 
THE  WINDWARD  PLANE  OF  A  20-PERCENT  BLUNT  CONE  (M.  =  15.  a  =  10“) 
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FIGURE  2ncl  COMPARISON  OF  APPROXIMATE  AND  EXACT  PRESSURE  COEFFICIENTS 
ON  A  20-PERCENT  BLUNT  CONE  {M«  =  15,a  =  10'.  REAL  GAS) 
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nCURE  27(d).  COMPARISON  OF  APPROXIMATE  AND  EXACT  TEMPERATURES 
ON  A  20-PERCENT  BLUNT  CONE  (M.=  15,a  =  10'.  REAL  GAS) 
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Figures  27{c)  and  27(d)  present  the  real  gas  pressure  coefficient  and 
temperature  comparisons  in  the  =  180°,  90°,  and  0°  planes.  The  worst  case 
comparisons  on  temperature  are  in  the  leeward  plane  where  the  deviation  is  as  high 
as  10  percent.  However,  as  alreadj'  noted,  the  windward  plane  maximum 
temperature  deviation  is  only  2  percent  or  less.  Figures  28(a)  and  28(b)  present  the 
blunt  cone  real  and  perfect  gas  pressures  and  temperatures,  respectively,  as  a 
function  of  Mach  number  for  a  =  0  and  15°  using  the  approximate  method.  The 
results  presented  are  those  near  the  base  of  the  cone,  x  =  1.5,  and  in  the  windward 
plane.  Several  points  are  worthy  of  note.  First  of  ail,  Figime  28(a)  reiterates  the 
negligible  effect  of  real  gas  conditions  on  pressure  at  all  Mach  numbers.  It  also 
illustrates  the  strong  pressure  increase  as  a  function  of  angle  of  attack  and  Mach 
number.  Figure  28(b)  illustrates  how  the  real  gas  affects  temperature  as  Mach 
number  increases.  It  is  also  interesting  to  note  that  while  temperature  differs  by  up 
to  several  thousands  of  degrees  due  to  angle-of-attack  effects  for  a  perfect  gas,  the 
real  gas  difference  is  a  maximum  of  1500°  at  Moo  =  20  and  a  =  15°  and  a  =  0°. 
Finally,  as  already  discussed,  the  present  approximate  code  gives  very  good 
agreement  for  inviscid  surface  properties  with  the  exact  Euler  solver,  ZEUS. 


Alp  =  0,  Perfect 

— ♦ — 

Alp  =  0,  Real 
Alp  =  15,  Perfect 
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Alp  =  15,  Real 
ZEUS 


FIGURE  28<a  t.  PRESSURES  PREDICTED  BY  .\PPROXIM.'\TE  METHOD  O.V  A  20-PERCENT 
BLUNT  CONE  .4S. A  FUNCTION  OF  .M.4CH  NUMBER  =  180\x  =  1. .Y  h  =  100k  fi  = 
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FIGURE  28(bi.  TENIPERATURES  PREDICTED  BY  APPROXIMATE  METHOD 
ON  A  20-PERCENT  BLUNT  CONE  AS  A  FUNCTION  OF  NL\CH  NUMBER 
«I>  =  180”,  X  =  1.5,  h  =  100k  ft) 
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The  third  case  considered  is  a  cone-cylinder-flare  (see  Figure  24).  Results  of  the 
pressures  and  temperatures  on  this  configuration  are  shown  in  Figures  29.  Figiu*es 
29(a)  and  29(b)  present  the  pressure  and  temperature  in  the  windward  plane 
respectively  as  a  function  of  distance  along  the  body  surface.  Note  the  present  theory 
shows  no  overexpansion  behind  the  shoulder  due  to  the  fact  that  MNT  is  used  to 
estimate  pressures  around  the  surface.  Also  note  in  Figure  29(b)  that  for  the  conical 
surfaces,  real  gas  teirperatures  are  lower  by  only  about  10  to  15  percent,  whereas  on 
blunt  nose  configurations  at  the  same  Mach  number,  temperatures  can  be  as  much  ais 
a  factor  of  2  lower  for  real  compared  to  perfect  gases.  Figure  29(c)  presents  the 
pressures  around  the  cone-cylinder-fiare  configuration  at  three  x  stations 
corresponding  to  a  point  on  the  cone,  cylinder,  and  flare.  Figures  29(d),  (e),  and  (f) 
present  the  temperature  comparisons  for  the  same  x  stations  of  Figure  29(c).  Note 
that  in  Figures  29(e)  and  29(f),  values  of  temperature  are  constant  in  the  leeward 
plane  area  due  to  the  Newtonian  assumption  that  Cp  =  0  in  shadowed  regions. 
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FIGURE  29(c!.  PRESSURE  DISTRIBUTION  ON  A  lO-CONS-CYI-INDER-FLAHr. 

(M=  =  is.  a  =  !G^  X  =  2.  i.  10: 
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FIGURE  29(0.  TE^fFERATUSE  DISTRIBUTION  AROUND  THE  FLARE 
OF  A  KrCONE-CYUI'rDER-FLARE 
(M*  =  10,  a  =  i0‘.  X  =  iO) 
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The  final  case  considered  is  the  20-percent  blunt  Von  Karman  ogive-cylinder- 
fin  configuration  shown  in  Figure  24.  The  angle-of-attack  0°  pressure  coefficient  and 
temperature  for  real  and  perfect  gases  are  shown  in  Figure  30  for  Moc  =  15.  Note 
again  the  excellent  agreement  up  to  the  location  of  the  tail  fins.  This  figure  was 
primarily  shown  to  indicate  that  while  the  present  approximate  engineering  code  is 
quite  applicable  for  preliminary  design,  it  does  not  have  the  physics  included  for 
detailed  interaction  effects.  These  effects  could  include  bow  shock  waves  intersecting 
fins,  or  fin  shock  waves  intersecting  the  body.  These  interaction  effects  can  cause 
local  “hot”  spots  and  more  detailed  analysis  codes  such  as  ZEUS  or  Navier  Stokes 
solvers  are  required.  The  fin  interaction  effects  are  shown  by  the  ZEUS  results  at 
X  =  9  calibers  to  the  end  of  the  body. 
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nOURE  30(b).  COMPARISON  OF  TEMPER.ATURE  PREDICTIONS  ON  A  20-PERCENT 
BLUNT  VON  KARMAN  OGI\'E-J‘',YUNDER-FIN  CONFIGURATION  (NU  =  15,  o  =  0“) 
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The  emphasis  to  this  point  has  been  on  accurate  values  of  inviscid  surface 
temperature  to  aliow  accurate  heat  transfer  analysis.  In  getting  accurate  values  cf 
heat  transfer,  accurate  values  of  pressure  prediction  were  lequired.  Tliese  accurate 
values  of  pressure  prediction  also  give  good  force  and  moment  predictions.  Tlie  final 
Figure  31  illustrates  force  and  moment  predictions  as  a  function  of  Mach  number  and 
angle  of  attack  compared  to  the  ZEUS  code.  Note  that  the  axial  force  only  includes 
wave  drag  since  ZEUS  at  present  is  an  inviscid  code.  Accuracies  on  wave  drag  and 
normal  force  are  within  10  percent,  and  center  of  pressure  near  angle  of  attack  zero 
within  8  percent  of  the  body  length  compared  to  the  full  Euler  code.  About  half  of  the 
error  in  wave  drag  is  due  to  the  high  fin  alone  predictions  using  the  present  strip 
theory  approach.  Figure  31(a).  Note  that  the  theory  does  predict  some  slight  changes 
in  forces  and  moments  due  to  real  gas  effects.  However,  except  for  veiry  specialized 
problems.  It  is  believed  that  these  effects  can  be  neglected,  particularly  in  an 
engineering  code  such  as  the  aeroprediction.  However,  as  already  noted  many  times, 
heat  transfer  analysis  definitely  needs  to  consider  the  real  gas  effects  when  flying  at 
any  appreciable  time  above  M  oe  —  6to8. 


NAVSWCTR  90-683 


AP  -  Perfect 
AP  -  Real 


ZEUS  -  Perfect 
ZEUS -Real 


ANGLE  OF  ATTACK 


FIGURE  31(b).  NORMAL  FORCE  COEFnCHENT  OF  A  2C^PERCENT  BLUNT  VON  KARMAN 
OGF/E-CYLINDER-FIN  CONFIGURATION  (^^«  =  15) 
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FIGURE  3i(c}.  CENTER  OF  PRESSURE  OF  A  20-PERCENT  BLUNT  VON  KARMAN 
OGIVE-CYUNDER-FINCONnGURATIONCM.  =  15) 
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FIGUBE  31  (e).  NORMAL  FORCE  COEFFICIENT  FOR  A  20-PERCENT  BLUNT  VON  KARMAN 
OGIVE-CYLINDER-FIN  CONFIGURATION  (a  =  5") 
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4.  SUMMARY 


New  methods  have  been  developed  to  compute  inviscid  surface  pressures  and 
temperatures  for  both  perfect  and  equiUbrium  chemically  reacting  flows  on  both 
pointed  and  blunt  bodies  of  revolution.  These  new  methods  include  an  improved 
Shock-expansion  Theory,  an  improved  MNT,  and  an  improved  method  for  angle-of- 
attack  effects.  Comparison  of  these  approximate  engineering  techniques  to  exact 
inviscid  computations  using  a  full  Euler  code  showed  the  following: 

1.  Agreement  on  the  critical  windward  plane  inviscid  temperatures  generally 
of  4  percent  or  better. 

2.  Agreement  of  inviscid  surface  temperatures  of  10  percent  or  better. 

3.  Agreement  of  axial  wave  drag  and  normal  force  of  10  percent,  and  center  of 
pressure  8  percent  of  body  length.  For  the  body  alone,  wave  drag  estimates 
are  generally  within  5  percent  of  exact  Euler  computations. 

A  new  real  gas  formulation  for  pressure  gradient  behind  a  comer  was  derived. 
However,  in  implementing  this  into  the  SOSET  it  was  found  to  be  of  little  value.  This 
is  because  at  high  Mach  numbers  the  exponential  decay  term  used  in  the  SOSET 
becomes  positive  requiring  SOSET  to  revert  back  to  ei^er  GSET  or  TCT.  It  was 
shewn  that  neither  of  these  theories  was  best  for  all  cases  and,  as  a  result,  a  user 
input  to  allow  a  choice  of  which  method  tc  use  was  considered  the  best  alternative  for 
use  of  the  traditional  SOSET. 

With  the  nev/  technology  developed,  the  NaVSWC  aeroprediction  code  can  now 
be  used  to  give  engineering  estimates  of  inviscid  surface  temperature  for  any  Mach 
number  of  interest.  These  approximate  temperatures  can  then  be  used  as  inputs  for 
more  detailed  heat  transfer  analysis. 


98 


NAVSWCTR 90-683 


S.  REFERENCES 


1.  Moore,  F.  G.,  CompiUational  Aerodynamics  at  NAVSWC:  Past,  Present,  and 
Future,  NAVSWC  TR  90-569,  October  1990, 

2.  Devan,  L.,  Aerodynam  'irs  of  Tactical  Weapons  lo  Mack  Number  8  and  Angle  of 
Attack  18(P:  Parti.  Thet^ry  and  Application,  NSWC  TR  30-346.  October  1980. 

3.  Devan,  L.,  and  Mason ,  L.,  Aerodynamics  of  Tactical  Weapons  to  Mack  Number  8 
and  Angle  of  Attack  i30®r  Part  II,  Computer  Program  and  Users  Guide,  NSWC 
TR  81-358,  September  1981, 

4.  Devan,  L.,  Mason,  and  Moore,  F.  G-,  Aerodynamics  of  Tactical  Weapons  to 
Mach  Number  3  and  Angle  of  Attack  I80\  AIAA  Paper  No.  52-0250,  20th 
Aero^jace  Sciences  Meeting,  Orlando,  FL, 

5.  Vincenti,  W.  G.,  and  Kruger,  C.  H.,  Introduction  to  Physical  Gas  Dynamics, 
John  Wiley  and  Sons,  Inc.,  New  York,  1965. 

6.  Jones,  J.  B..  and  Hawkins,  G.  A.,  Engineering  Thennoifynamics,  John  Wiley 
and  Sons,  Inc.,  New  York,  1960. 

7.  Liepmann,  H.  W.,  and  Roshko,  A.,  Elements  of  Gasdynamics,  John  Wil^  and 
Sons,  Inc,,  New  York,  1957. 

8.  Andersoii,  J.  D.,  Modem  Compressible  Flow,  McGraw-Hill  Book  Company, 
1982. 

9.  Truitt,  R.  W,,  Hype:^sonic  Aerodynamics,  The  Ronald  Paress  Company,  New 
York,  1959. 

10.  Hayes,  W.  D.,  and  Probsteia,  R.  F.,  Hypersonic  Flow  Tlieory,  Academic  Press, 
New  York,  1966. 

11.  Hilsenrath,  J.,  Klein,  M.,  and  Woiiey,  H.  W.,  Tables  cf  TTienriodynamic 
Prcnerti^  of  Air  Including  Dissociation  and  Ionization  from  1500°K  to 
150b0”K,  A£DC  TR-59-20. 

12.  Lees,  Lester,  Hypersonic  Flow,  Inst.  Aero.  Science,  Preprint  No.  554, 1955. 

13.  Ames  Research  Staff,  Equations,  Tables,  and  Charts  for  Compressible  Flow, 
NACATR 1135, 1953. 

14.  Eggers,  A.  J.,  Syvertson,  C.  A.,  and  Kraus,  S.,  A  Study  of  Inviscid  Flow  About 
Airfoils  at  High  Supersonic  Speeds,  NACA  Report  1123, 1953. 


99 


NAVSWCTR  90-683 


16.  Eggers,  A.  J.,  and  Savin,  it  C-,  A  Unified  Tivo~Dimensional  Approach  to  the 
Calculaiion  c/  Three-Dimensional  Hypersotiic  Flows,  with  Application  io  Bodies 
of  Revolution,  NAC  A  Report  1249, 1955. 

16.  Eggers,  A.  J.,  and  Savin.  R.  C.,  Approximate  Methods  for  Calculating  the  Flow 
About  Nonlifting  Bodies  of  Reitolution  at  High  Supersonic  Airspeeds,  N  ACA  TN 
2579, 1951. 


17.  Savi»i,  R.  C.,  Application  of  the  Generalized  Shock  &^iansian  Method  to  Inclined 
Bodies  ofReodUiiian  Ttaoeling  at  High  Supersonic  Airspeeds,  NACA  TN  3349, 
1955. 

18.  Taylor,  G.  I.,  and  MacCall,  J.  W.,  The  Air  Pressure  on  a  Cone  Moving  at  High 
Speed,  Proc-  Royal  Society  (Londtsi)  Ser.  A.,  Vol.  139, 1933,  pp.  278-311. 

19.  Syvertson,  C.  A.,  and  INanis,  D.  H.,  A  Second-Order  Shock  Expansion  Method 
Applicable  to  Bodies  of  Reoolution  Near  Zero  Lift,  NACA  TR 1323, 1957. 

20.  Jackscn,  C.  M..  Sawyer,  W  C.,  and  Smith,  R.  S.,  A  Method  for  Determining 
Surface  Pressures  on  Blunt  Bodies  of  Revolution  at  Small  Angles  of  Attack  in 
Supersonic  Flow,  NASA  TN  P-4865, 1968. 

21.  DeJamette,  F.  R.,  and  Jones,  K.  M.,  Development  of  a  Computer  Program  to  Cal¬ 
culate  Aerodynamic  Characieristics  of  Bodies  and  Wing-Body  Combinations, 
NSWC/DLTR-3829,  April  1978. 

22.  DeJarnatte,  F.  R.,  and  Ford,  C.  P.,  A  New  Method  for  Calculating  Surface 
Pressures  on  Bodies  at  an  Angle  of  Attack  in  Supersonic  Flow,  Proceeding  of  the 
11th  Navy  Svmposium  os  Aerobeliistics,  Vol.  il,  .August  22-24, 1978,  pp.  199- 
231. 

23.  Ds-Jamette,  F.  R.,  Ford,  C.  P.,  and  Young,  D.  E.,  A  New  Method  for  Calculating 
Surface  Pressures  on  Bodies  at  an  Angle  of  Attack  in  Siqrersonic  Flow,  AIAA 
Paper  No.  79-1552,  AIAA  12th  Fluid  and  Plasma  Dynamics  Conference, 
Williamsburg,  V  A,  July  1974. 

24.  WittiifiT,  Charles  E.,  and  Curtis,  James  T..  Normal  Shock  Wave  Parameters  in 
Equilibrium  Air,  Cornell  Aero  Lab.  Report  CAL-III,  November  1961. 

25.  Mairone,  P.  V.,  Normal  Shock  Waves  in  Air:  Equilibrium  Composition  omI 
Flow  Parameters  for  Velocides  from  26,000  to  50,000 ft/sec,  CAL  Report  No.  AG- 
1723-.A-2,  August  1962. 

26-  Heims,  Steve,  Prandd-Meyer  Expansion  of  Chemically  Reacting  Gases  in  Local 
Chemical  and  Thermodynamic  BguUibrium,  NACA  TN  4230,  March  1958. 


NAVSWCTR9C-683 


27.  Moeckei,  W.  £.,  Oblique-shock  Relations  at  Hypersonic  Speeds  for  Air  in 
Chemical  Equilibrium,  NACA  TN  3895,  January  1957. 

28.  Wardiaw,  A.  B.,  Baitakis,  F.  P.,  Solomon,  J.  M.,  and  Haekerman,  L.  B.,  An 
Inviscid  Computational  Method  for  Tactical  Missile  Configurations,  NSWC  TR 
81-457,  December  1981. 

29.  Wardiaw,  A.  B.,  1*11010,  F.  J.,  Solomon,  J.  M.,  and  Baitakis,  F.  ?.,  Inviscid 
Multiple  Zone  Calculations  for  Supersonic  Taeticai  Missiles,  AIAA  Paper  84- 
2099, 1984. 

30.  Moore,  F.  G.,  Body  Alone  Aerodynamics  of  Guided  and  Unguided  Projectiles  at 
Subsonic,  Transonic,  and  Supersonic  Mach  Numbers,  NWL  TR-3796,  Novem¬ 
ber  1972. 

31.  Moore,  F.  G.,  and  McKerly,  B.,  Aerodynamics  of  Guided  and  Unguided 
Weapons:  Parti—  Theory  and  Application,  NWLTft-3018,  December  1973. 

32.  Moore.  F,  G.,  and  Swanson,  R.,  Aerodynamics  of  Tactical  Weapons  to  Mach 
Number  3  and  Angle  of  Attack  15  Degrees:  Part  I  —  Theory  and  Appiicaticn, 
NSWC  TR-3584,  February  1977. 

33.  Anderson,  J.  A.,  Hypersonic  and  High  Temperature  Gas  Dynamics,  McGraw- 
Hill  Series  in  Aeronautical  and  Aerospace  Engineering,  1989. 

34.  Tannehili.  J.  C.,  and  Mugge,  P.  H-,  Improved  Curve  Fits  for  the  Thermodynamic 
Properties  of  Equiiibrium  Air  Suitable  for  Numerical  Computation  Using  Time- 
Dependent  or  Shock-Capturing  Methods,  NASA  CR-2470, 1974. 

35.  Srinivasan,  S.,  Tannehili,  J.,  and  Weilmusnster,  K.,  Simplified  Curve  Fits  for 
the  Thermodynamic  Properties  of  Equilibrium  Air,  Report  ISU-ERI-AMES 
86401,  Eng.  Research  Institute,  Iowa  State  University,  hm.es,  L4,  June  1986. 

36.  DeJamette,  F.  R.,  Ford,  C.  P.,  Young,  D.  E.,  Xalcuiation  of  Pressure  os  Bodies 
at  Low  Angies  of  Attack  in  Supersonic  Flow,”  JSR  Aftide  No.  79-1552R,  Vol. 
17.  No.  6,  November-Decsmber  1980. 

37.  Hudgins,  Henry  S.,  Jr.,  Supersonic  Flow  About  Right  Circular  Cones  at  Zero 
Yaw  in  Air  at  Chemical  Equilibrium,  Part  I  —  Correlaticn  of  Flow  Properties, 
TM  1493,  Picatisny  Arsenal,  August  1965. 

38.  Jones,  D.  J.,  Numerical  Solutions  of  the  Flow  Field  for  Conical  Bodies  in  a 
Supersonic  Stream,  National  Research  Council  of  Canada  Report  LR-507,  July 
1968. 


161 


NAVSWCTR  90-683 


39.  Van  Driest,  E.  R.,  “Turbulent  Boundary  Layer  in  Compr^sible  Fluids,”  JAS, 
Vol.  18,  No.  3,  pp.  145-160, 216. 

40.  Alien,  J.  H,,  and  Perkins,  E.  W.,  Characteristics  of  Flow  Over  Inclined  Bodies  of 
Revolution,  NAC  A  RM  A50L07,  Mofifett  Field.  CA,  1965. 

41.  Wardlaw,  Ai^rew,  and  Baitakis,  Frank,  ZEUSBL:  Integral  Boundary  Layer 
Procedure  for  Tactusd  Missiles,  NAVSWC  technical  r^wrt  in  pr^)aration. 

42.  Armistead,  Michael,  and  Mclnviile,  R.  M.,  A  UseFs  Guide  to  the  BLZEUS 
Aeroheating  Analysis  Code,  NAVSWC  t»dinical  note  in  preparation. 

43.  Bender,  D.  R. ,  A  Miniature  Version  of  the  JA  70  Aerodynamic  Heating  Computer 
Program,  H  800  (MINFVER),  McDonnell  Douglas  MDC  G0462,  (COSMIC 
PROGRAM  3  MFS-21951),  June  1970. 

44.  Truit,  R.  W.,  Fundamentals  of  Aerodynamic  Heating,  Tbe  Ronald  Press,  New 
York,  Copyright  19^. 

45.  Rubesin,  M.  W.,  and  Johnson,  A.  K..  “A  Critical  Review  of  Skin-Friction  and 
Heat  Transfer  Solutions  of  the  L-aminar  Boundary  Layer  of  a  Flat  Plate,”  Trans, 
of  the  ASME,  Vol.  71,  No.  4,  pp.  383-388,  May  1949. 

46.  Eckert,  E.  R.  G.,  “Engineering  Relations  for  Heat  Transfer  and  Friction  in  High 
Velocity  Laminar  and  Turbukst  Boundary-Layer  Flow  Over  Surfaces  with 
Constant  Pressure  and  Temperature,”  Trans,  of  the  ASME,  Vol.  78,  No.  6, 
p.  1273,  August  1956. 

47.  Priolo,  F.  J.,  and  Wardlaw,  A.  B.,  Jr.,  “High-temperature  Effects  for  Missile- 
type  Bodies  Using  the  Euler  ^iver,  ZEUS,”  AIAA  Paper  No.  91-3259, 
presented  at  1991  AIAA  Applied  Aerodynamics  Conference,  Baltimore, 
Maryland,  September  1931. 


48.  Morrison,  A.  M-,  Solomon,  J.  M,,  Ciment,  M.,  and  Ferguson,  R.  E.,  “Handbook  of 
Inviscid  Sphere-cone  Flow  Fields  and  Pressure  Distributions:  Volume  I,” 
NSWC-WOUTR  75-45,  December  1975. 


102 


NAVSWCTR 90-683 


a 

b 

C,,C2 

Cf^ 

Ch 

Cp 

Cp,  Cy 
e 
H 
h 

k 

M 

P 

Pr 

q 

qw 

R 

Rn 

s 


6.  SYMBOLS  AND  DEFINITIONS 
speed  of  soiind  (fb'sec) 

distance  between  body  surface  and  streamline 
left  and  right  running  characteristic  coordinates 
local  skin  friction  coefficient 
heat  transfer  coefficient 
pressure  coefficient 

specific  heats  at  constant  pressure  and  volume 

internal  energy  {ft2/3ec2) 

total  enthalpy  (fl2/sec2) 

specific  entlialpy  (ft2/sec2)  and  altitude  (ft) 

adiabatic  wall  enthalpy 

Kelvin  temperature  (d^ 

Mach  number 
pressure  (Ib/ft2) 

Prandtl  number 
heat  flux 
heat  flux  at  wail 

gas  constant  [for  air  R  =  1716  (fl-ib)  /  {si«g-®R)  ] 

Reynolds  number 

recovery  factor 

entropy  (fl-ib)  /  (slug  -®R) 

streamline  coordinate  system 
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T 

T 

T 

*B 

t 

U,.U3 

V 

VnVe 

Vmax 

x,y,z 

z 

a 

P 

y 

8eq 

6 

®Ci  ®s 

A 

y 

Pi 

V 

p 


temperature  (degrees  Kelvin  or  Sankine) 
temperature  of  a  perfect  gas 
temperature  of  a  real  gas 
time 

velocity  (ft/sec) 
total  velocity 

velocity  along  and  normal  to  a  conical  ray 
maximum  velocity 

cylindrical  coordinate  ^stem  with  ^  =  0  =  leeward  plane 
rectangular  coordinate  ^stem 
compressibility  factor  (=  1  for  a  perfect  gas) 
ai^le  of  attack  (deg) 

ratio  of  specific  heats  (Cp  /  C  v) 

angle  used  in  Newtonian  Theory  (=  angle  between  velocity  vector 
and  local  body  slope) 

local  body  slope 

cone  half-angle  and  shock  angle,  req)ectively 

pV2/p 

Mach  angle,  p  =  sin  -l  (1  /  M) 
coefficient  of  viscosity 
Prandtl-Meyer  angle  (deg) 
density  (slugs/!t3) 
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shock  wave  angle  (d^) 

exponent  used  in  Second-order  Shock-expansion  Theory 

Subscripts 


adiabatic  wall 
cone 

dissociation  energy 
electronic  excitation  energy 
frozen  flow 
local  conditions 
match  point 
total  and  reference 
rotational  energy 
stagnation  conditions 
translation  energy 
vibration  energy 
freestream  conditions 

Superscripts 

reference  conditions 
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APPENDIX  A 

DETERMINATION  OF  ANGLE  8  USED  IN  NEWTONIAN  THEORY 
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DETERMINATION  OF  ANGLE  8  USED  IN  NEWTONIAN  THEORY 


Refer  to  the  nomenclature  of  Figure  A-1. 


z 


z 


WHERE:  r2  >  yZ  z2 


FIGURE  A-i.  NGMExN’CLATUEE  USED  FOR  DETERMINATION  OF  ANGLE  5,, 


Note  that  Irom  the  definition  of  the  dot  product  of  two  vectors,  one  can  write 


e  =  {VJ  i  ^  lcos(90 -r  5 


(A-l) 


Also 


A  A 

=  (\^  ^cosQii  +  (V  sina}k 


(A-2) 


and 


A 

e 

B 


A  ,  A  A 

— sio6i  +  (sin^cos6ij  +  cos6}k 


{AS) 
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A  A  /V  A 

=  cos6i  4-  airiOsin^  j  4-  sinBoos^k 


=  cos^j  —  sin4»k 


(A-4) 

(A-5) 


Substituting  Equations  (A-2)  and  (A-3}  into  Equation  (A-1)  and  performing  tiie 
indicated  operations,  one  obtams: 


smb 

eq 


=  sindcosa  —  anacos^tcosd 


(A-6) 


This  is  the  angle  that  is  used  for  a  general  three-dimensional  (3-D)  point  on  a  blunt 
surface  at  angle  of  attack  in  the  Newtonian  Hieory.  Here  ^  =  0  is  the  leeu  ard  plane 
and  ^  =  180  is  the  windward  plane.  If  the  reverse  is  true  (as  many  refe<  ences  use)^ 
the  negative  sign  of  Equation  (A-6)  becomes  positive. 


I 

I 


A-4 


NAVdWCTRSO^ 


APPENDIX  B 

PRESSURE  DISTRIBUTION  ON  A  SPHERE 
IN  HYPERSONIC  FLOW 
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PRESSURE  DISTRIBUTION  ON  A  SPHERE  IN  HYPERSONIC  FLOW 


The  Modified  Newtonian  pressure  distribution  has  been  used  by  numerous 
investigators  to  predict  pressure  distributions  on  blunt-nosed  bodies  in  hypersonic 
flow.  A  close  examination  of  this  method  for  a  sphere  showed  that  for  Mod  >  3.5  it 
overpredicts  pressures  for  x/R,^-  <  0.5625  and  underpredicts  pi;esures  for  most  of  the 
region  beyond  x/R,^  =  0.5625.  Here,  x  is  the  axial  coordinate  with  x=0  at  the 
stagnation  point,  and  is  the  radius  of  the  sphere.  However,  the  Modified 
Newtonian  method  predicts  the  correct  pressure  at  x/R^  =  0.5625  for  all  3.5  <  M®  ^ 
30.  Figure  B-1  shows  tiicse  features  for  Ma>  =  3.5  and  5,  but  they  hold  also  for  the 
higher  Mach  numbers. 

A  second  feature  of  the  Modified  Newtonian  (MN)  method  is  that  at  x/R^^  = 
0.5625  the  ratio 


«q 


1 

2 

MN 


varies  very  little  with  M®.  The  value  of 


dp 


was  obtained  numerically  from  CFD  data  in  Reference  B-1.  The  angle  of  Sgq  is  the 
inclination  of  the  surface  with  respect  to  the  axis  of  the  sphere,  where  8eq  =  n/2  at  the 
stagnation  point.  Table  B-1  shows  the  numerical  values. 


TABLE  B-1 .  VALUES  OF  PRESSURE  GRADIENT  PARANIETER 
AS  FUNCTION  OF  MACH  NU.MBER 


M® 


<Ip 

d£ 


/ 


d  p 

p  — 1  at  x/R„  =  0.5625 

d6  Imn  ^ 


C.5 

1.469 

5.0 

1.519 

10.0 

1.516 

20.0 

1.523 

30.0 

1.524 

Morrison,  A.  M..  Solomon,  J.  M.,  Ciment,  M.,  and  Fer^son,  R.  E.,  Handbook  of  Invisctd  Sphere 
Cone Fiowfields and Pr^sure Distributions:  VoL  /,  NSWCAVOL TR  75-45,  December  1975. 
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o  MOOIHED  NEWTONIAN.  Ma>=3.s 
□  MODIFIED  NEWTONIAN.  Ma>=5.0 
-  CFD.  REFERENCE  B-1 


S/R 


nOURE  B-1.  PRESSURE  ON  A  SPHERE 
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Thus,  the  position  x/Rj^  =  0.5625  will  be  called  the  match  point  since  the  pressiure  is 
Modified  Newtonian  pressure  and  the  pi^sure  gradient  can  be  calculated  fi^om  the 
ratio 


kI\ 


dp 


1 

I- 

i2 


dS  iuM 

eq 


It  has  been  determines  that  using  a  numerical  value  of  1.5  for  this  ratio  for  all  M®  > 
3.5  gives  accurate  results. 


For  the  region  0  ^  x  Rj,  <  0.5625,  a  more  accurate  prediction  of  the  pressiure 
can  be  calculated  firom  the  equation 


AC  =kcos“5  (cos 5  —cos  (5  )„) 

P  tq  tq  «q 


(B-l) 


where 


“p  =  Vn  -  Cp 


and  (6eq  )j|  =  25.95,  which  is  the  value  of  6eq  at  x/R,;  =  0.5625.  Both  k  and  m  are 
constants  to  be  determined.  Note  that  Equation  (B-D  gives  ACp  =  0  at  both  =  0 

and  =  0.5625.  In  addition,  it  gives  d(ACp)  /  d  =  0  at  x/R,,  =  0.  It  now 
remains  to  determine  expr^sions  for  k  and  m. 

The  derivative  of  Equation  (B-l)  is 


d  (AC  )  j  , 

. =  ksan  6  cos*°  *5  ;  {m  +  Dcns  8  —  mens  'b  I 
d6  ■  «q 


(B-2) 


Apply  Equation  (B-2)  at  the  match  point,  =  0.5625 


d<AC  )  f  dC  ^  dC 

?  i 


dS 


- 


.do  /|iN  dS 

«J  eq 


^  =  -ksin(6  ;•  cos“(5  )„ 

M  cfi  n 


(B-3) 


Now  at  the  match  poini,  the  derivative 


dC 
_ ^ 

dS 

eq 

can  be  calculated  firom  the  ratio  given  earlier;  i.e., 
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<ip 

d8 


=  1.5  ? 


«l 


1 

dp  ]2 


dS„  JMti 


which  can  also  be  written  as 


dC  1- 
'  2 


d5 


=  1-5  j(c  +  )  — 2.1 

\\  P  y^2  /  d5„  1 


MN 


(B-4) 


where 


(C  =  C  sin^ 

P  MN  24 


(B-5) 


Equation  (B-4}  can  be  used  in  Equation  (B*3)  but  that  only  gives  one  equation  for  the 
two  unknowns,  k  and  m.  Thus,  an  additional  equation  is  needed. 

A  third  feature  of  Modi^ed  Newtonian  pressures  is  obtained  by  observing  that 
the  data  in  Reference  B>1  show  that  ACp  is  a  maximum  at  x/Rj,  =  0.25  for  all  Ma>  > 
3.5.  Thus,  Equation  {B-2)  can  be  used  to  get 


B.-1  [  I 

— ; -  =0  =  k»xiS  cos*"  *5  i  {ro-i-l}  cos  5  —  iac*»(5  )„  j 

d6  *s  r  «i  "q  M I 

eq 


(B^) 


Apply  this  equation  atx/H^  =  0.25  to  obtain 


m 


=  ( 


oos6 


ods{6^)jj-  os5^ /,fli^=025 


=  2.78 


(B-7) 


Finally,  substitute  Equation  (&4}  into  {B-3)  and  apply  at  the  match  point  (z/R^  = 
0.5625,  (fieq),^  =  25" ,  )  to  get 


k  = 


8in(S  )  c3S®(5  ),, 

cq  H  cq  M 


(B-8) 


Ail  terms  on  the  right  side  can  be  calculated  for  a  given  Moo.  Substituting  for  the 
parameters  af  Equation  (B-8)  at  the  match  point  where  (8eq)|f  =  25.95”,  one  obtains 
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k  = 


—  2.416  C  4.60S  j  0.1SS7  + 


i.l24 


{B-9> 


Wi^  m  and  k  determined,  an  accurate  Cp  can  be  calcuiated  from  Equation  <B-1)  for 
z/Rjj  ^  0.5625. 

This  method  has  been  applied  to  a  sphere  for  3.5  ^  Ma  ^  3C.  The  results  are 
compared  to  CFD  results  from  Reference  B-1  for  0.5625.  At  Ma>  =  3.5  the 

maximum  error  in  pressure  is  3  percent  at  x/Rj^  =  0.25.  For  Ma>  >  3.5,  the  maximum 
error  is  less  than  1  percent  and  in  most  cases  less  than  0.5  percent,  in  confrast,  the 
Modifred  Newtonian  pressures  had  maximum  errors  from  7  to  10  percent  Equation 
(B>1)  can  be  used  for  x/R^  up  to  about  0.7  with  errors  less  than  3  percent.  A  better 
prediction  for  x^Rj^  >  0.5625  can  be  obtained  by  using  the  second-order  shock- 
expansion  method. 

Ihe  new  method  gives  a  much  mors  accurate  prediction  of  surface  pressures  on 
a  sphere  in  hypersonic  flow  than  Modified  Newtonian  Theory.  The  additional 
computatkmai  effort  is  negligible  for  even  the  smallest  PCs.  Results  of  this  new 
technique  are  shown  in  the  Results  and  Discussion  section  of  this  report. 
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APPENDIX  C 

COBIPUTATION  OP  PROPERTIES  ACROSS  SHOCK  WAVES 
IN  FROZEN  AND  EQUIUBRIUM  CHEMICAIXY 
REACTING  FLOWS 
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COMPUTATION  OF  PROPERTIES  ACROSS  SHOCK  WAVES 
IN  FROZEN  AND  EQUIUBRIUM  CHEMICALLY 
REACTING  FLOWS 


To  compute  the  Sow  over  2-D  or  3-D  configurations  using  the  Shock-E^ansion 
Theory  requires  an  initial  starting  soiuticn.  This  solution  is  provided  by  computing 
the  Sow  field  variables  acr<»s  a  normal  (in  the  case  of  a  blunt-nose  body  or  blunt- 
leading-edge  wing)  or  oblique  ^ock  for  ^mrp-nosed  bodies  or  wing  leading  edges.  If 
the  flow  is  low  enough  in  Mach  number  so  that  real  gas  effects  are  n^ligible  or  if  one 
is  only  inter^ted  in  approximate  prediction  of  forces  and  moments,  a  perfect  gas  can 
be  assumed  and  the  shock  wave  relations  of  Reference  C-1  can  be  used  directly.  This 
last  situation  results  firom  the  fad  that  real  gas  effects  have  a  fairly  smaU  effect  on 
pressure  (usually  less  than  10  percent)  but  can  have  a  large  effect  on  dem-ity  and 
temperature.  Since  we  are  interested  in  temperature  profiles  along  the  body  as 
izq>ut&  for  structural  analysis,  we  must  consider  both  normal  and  oblique  shock  wave 
computations  in  real  gases. 


NORMAL  SHOCK  WAVES 

The  conservation  of  mass,  momentum,  and  energy  —  Equations  (43),  (44),  and 
(45)  of  file  main  text  of  this  report—for  steady  (d/3t  =  0),  adiabatic  (q  =  0)  fiow  with 
no  body  forces  (f  =  0),  redu<»  to  the  following  for  flow  across  a  normal  shock  wave  (see 
Figure  C-1): 


(C-1) 

p,+p,u;=p,+P3u* 

(C-2) 

u-  u- 

h  +  —  =  h-  +  — 

(C-3) 

1  2  2  2 

Solving  Equations  (C-2)  and  (C-3)  for  p2  and  h2,  while  using  Equation  (C-1),  there  is 
obtained 


P2  =  Pl  +  PlUf(i-Pj/p2) 

{C-4) 

h2=hj  +  Uf/211  -(pj/p/j 

C-5) 

C-3 
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Note  that  aii  terms  on  the  right-hand  ^e  of  Equations  (C-4)  and  (C~5)  are  known 
except  for  g^.  To  solve  Equations  (C-4}  and  (C-5),  one  guesses  a  value  of  /  p2-  A 
value  of  0  or  0.1  is  adequate.  Equations  (C-4)  and  (C-5)  are  then  solved  for  P2  and  h2. 
Knowing  and  a  new  value  of  is  computed  from  the  eqiiilibrium  dmmistry 
model  discussed  in  the  Simplified  Procedures  for  Air  section  (section  2.3.2).  New 
values  of  P2  and  K  are  then  computed  and,  once  again,  new  values  of  P2  computed 
from  the  equilibrium  chemistry  model.  Ihis  process  is  repeated  until  the  change  in  p2 
is  small  and  within  the  desired  accuracy  for  convergence.  This  then  defines  the 
correct  values  of  P2>  h2,  and  p2,  behind  the  ^ock.  T2,  a2  can  be  computed  from  section 
2.3.2  once  02  is  computed  from  Equation  (67).  The  compressibility  factor  Z  is  defined 
as: 


2  = 


P 

pRT 


(C-6) 


For  a  perfect  or  thermally  perfect  gas,  Z  is  one  but  for  a  real  gas  that  undergoes 
chemical  reactions,  Z  represents  the  level  of  dissociation  that  takes  place.  Since  R  is 
the  universal  gas  constant  and  p,  p,  and  T  are  computed  by  Equations  (C-4)  and  (C-5) 
and  Section  2.3.2,  Z  is  known.  Also 


U2  =  (Pi/p2)U, 

(C-7) 

Yj  =  a^P^P 

(C-8) 

M,  =  U,/a, 

(C-9) 
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‘Xlie  Normal  Shock  Solution  for  Equilibrium  Air  is  given  in  References  C-2  and  C-S. 
Figure  C-2  gives  T^/  as  a  function  of  fireestaream  pressure  and  velocity.  Note  the 
lar^  di^erenc^  between  the  perfect  and  real  gas  temperatures  as  temperature 
increases  and  as  atmcepheric  pressures  decrease  (higher  altitude). 


FIGURE  C-2.  INFLUENCE  OF  FREESTREAM  PRESSURE  AND  VELOCITY  ON 
NORMAL  SHOCK  TEMPERATURE  FOR  EQUILIBRIUM  AIR  (T,  =  225‘K) 
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Fcr  the  flow  field  over  a  blunt  body,  the  properties  vary  between  the  shock  and 
body.  Hence,  frozen  flow  would  consist  of  computing  real  gas  properties  immediately 
behind  the  shock  and  then  treating  the  flow  field  as  a  perfect  gas  thereafter,  with 
constant  values  of  y  =  yf  and  compressibility  factor  Z  =  Zf.  Equilibrium  chemically 
reacting  flow  would  allow  different  values  of  these  parameters  as  a  frmction  of 
distance  between  the  body  and  shock.  However,  if  one  is  only  interested  in  properties 
along  the  body  surface  from  a  method  such  as  Modified  Newtonian  combincKl  with 
3CSET,  only  the  stagnation  streamline  that  wets  the  body  surface  needs  to  be 
examined.  Reference  C-2  showed  that,  along  the  stagnation  streamline,  y  and  Z 
change  very  little.  Hence,  effectively  the  flow  is  frozen  at  values  that  exist  behind  the 
shock  (within  1  or  2  percent)  between  the  shock  and  body  stagnation  point.  This 
assumption  will  be  made  in  the  present  work. 


OBUQUE  SHOCK  WAVES:  2-D  OR  WEDGE  FLOWS 

Examining  Figure  C-3,  it  soon  becomes  clear  that  one  can  use  the  previous 
normal  shock  relations  (C-l)  through  (C-5)  to  solve  for  flew  parameters  behind  an 
oblique  shock  if  the  and  components  of  flow  velocity  normal  to  the  shock  wave 
are  replaced  with  their  equivalent  components  normal  to  the  oblique  shock  wave. 
That  is 


U  =U,-sino  (C-10) 

1 

U  =U,sin(o-e)  (C-11) 

2 

Here,  6  is  the  flow  direction  which,  for  attached  flow  on  a  wedge,  is  constant  and 
equal  to  the  wedge  angle.  Equations  (C-l),  (C-4),  and  (C-5)  now  become,  with  these 
substitutions 

U ,  sin  o  =  p,  U,  -sLn  <o  -  0}  (C-12) 


+  PjUjSin^otl  -  9^!?^  (C-13) 

U'sin^o 

hj  =  +  - 2 - ‘  ^  ‘*2^'  ^  (C-14) 


Also,  since  =  Uij  from  the  conservation  of  momentum,  then 
UjOoso  =  U,cos(o  —  0) 


cos(o  —  0) 
U,  cx>s(a) 

O 


(C-15) 
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Combining  Equation  (C-15)  with  Equation  (C-12),  there  is  obtained 


_  tan  (o  —  6) 


(C-16) 


Ejqianding  tan  (a  -  8)  using  a  higoncmetric  identity,  going  through  some  algebraic 
manqiulation,  and  taking  ihe  weak  shock  solution,  one  obtains 


/  /  Po  \  1  Po 

tano  =  <p  /p  —  l)/{2tan6)—  V  (  i - ) - -  (0-17) 

Pl'  4Un^e  Pi 


To  solve  for  the  How  field  properties,  a  value  of  pi/p2  is  once  again  guessed.  The  shock 
angle  is  computed  from  ^nation  (C-17).  With  &is  value  of  o  and  pi/p2  known,  P2 
and  h2  are  computed  firom  Equations  (C-13)  and  (C-14).  A  new  value  of  p2  is 
computed  from  the  process  of  Section  2.3.2  and  an  iteration  takes  place  until 
successive  changes  in  P2  are  within  some  error  bound.  This  defines  the  values  of  o, 
P2,  P2>  and  h2  behind  the  shock.  T2  and  "a”  are  then  computed  firom  Section  2.3.2  and 
Z,  U2  from  Equations  (C-€},  <C-7)  and  (C-15).  Local  Mach  number  and  isentropic 
exponent  y  are  then  computed  from  Equations  (C-8)  and  (C-9). 
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It  should  be  pointed  out  that  the  iaentrt^e  aqxment,  defined  as 

V= 

la(inp)l, 

is  the  same  as  the  gamma  of  Equation  (C~8).  This  is  different  than  the  “equivalent” 
gamma  defined  in  Reference  C4  as 


h 

Y  =  - 
e 


Here  the  word  equivalent  in  quotes  is  used  because  y  =  y  for  a  perfect  gas.  However, 
for  a  real  gas  where  dissociation  is  present,  the  isentropic  exponent  used  in  true  speed 
of  sound  computations  is  different  than  the  “equivalent”  gamma. 

Since  the  body  is  a  wedge,  ffie  flow  properties  behind  the  straight  shock  are 
constant  throughout  the  flow  field.  Hence,  the  pressure  P2  is  also  the  pressure  on  the 
surface  which  can  be  used  to  compute  forces  and  moments.  Also,  Frozen  and 
Equilibrium  Chemically  Reacting  Flows  give  the  same  results  since  all  properties  are 
constant  between  the  body  and  shock.  Figure  C4  gives  the  results  ofHie  shock  angle 
versus  wedge  half  angle  at  an  altitude  of  lOOK  ft  for  two  velocities  for  a  perfect  and 
real  gas  (either  firozen  or  equilibrium  chemically  reacting).  Note  the  small  change 
between  the  perfect  gas  and  real  gas  at  small  flow  deflections  and  Mach  numbers. 
Figure  0-5  is  an  example  of  the  shock  waves  for  a  particular  case  where  the 
conditions  are  such  to  show  a  large  difference  between  pmTect  and  real  gas  computa¬ 
tions.  This  same  simplified  flow  field  does  not  hold  true  if  the  body  is  axisymmetric 
as  will  be  discussed  shortly . 

Figures  0-6  through  08  show  the  temperature,  density,  and  pressure  for  the 
same  lOOK  altitude  condition  at  two  velodties  as  a  function  of  wedge  angle.  Note  the 
major  effect  of  real  gases  is  in  the  temperature  and  density  at  higher  velocities  and 
large  wedge  angles.  The  pressure  is  only  affected  slightly. 
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ANGLE 

(•) 


FIGURE  C-4.  DEFLECTION  ANGLE/WAVE  ANGLE  VELOCITY  DIAGRAM 
FOR  OBUQUE  SHOCKS  (ALTITUDE  =  100k  R) 
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FIGURE  C-S.  COMPARISON  OF  OBUQUE  SHOCK  WAVES  FORA  CALORICALLY 
PERFECT  GAS  VERSUS  AN  EQUIUBRIUM  CHEBOCALLY  REACTING  GAS 
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FIGURE  C-€.  TEMPERATURE  IN  THE  SHOCK  LAYER  FORA  WEDGE 

(ALTITUDE  =  100k  R) 
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FIGURE  C-7.  DENSITY  IN  THE  SH(X:E  LAYER  FOR  A  WEDGE 
(ALTITUDE  =  100k  R) 
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FIGUREC-8.  PRESSURE INTHESHOCKLATERFORA WEDGE 
(ALTITUDE  =  100k  ft 


OBUQUE  SHOCK  WAVES:  AXISYMMETRIC  CONICAL  FLOWS 


The  flow  field  between  fiie  shock  and  conical  sur&ce  is  not  constant  here  as  was 
the  case  with  the  wedge.  However,  fiow  properties  are  constant  along  rays 
emanating  firom  the  conical  tip  for  a  perfect,  frozen,  or  equilibrium  chemically 
reacting  gas.  Furthermore,  the  entropy  is  constant  throu^out  the  flow  field  between 
the  shock  and  body  for  these  flow  types.  The  continuity  and  momentum  equations  in 
spherical  coordinates  for  this  type  of  How  bea}meU~^ 


(C-18) 


dV 

— I=V*  (C-19) 

de  « 


dV 


8 


de 


+V-C0t6- 

0 


(C-20) 
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Figure  C-9  shows  the  geometry  and  nomenclature  used.  These  equations  represent  a 
set  of  three  firstnrder  ordinary  differential  equations  for  five  unknowns.  'Hie 
remaining  two  equations  come  firom  hi^-temperature  thermodynamic  properties 


p  =  p(p.S)  (0-21) 

a  =  a  (p,  S)  (C-22) 


Perfect  Gas 

For  a  perfect  gas»  Equations  ((^>21)  and  (022)  are  replaced  by 


a 


2 


((^23) 


P  =  ^  (C-24) 

a 

Substituting  Equations  (C-23)  and  (C-24)  into  Equations  (C-18)  through  ((%20)  and 
performing  the  algebra,  one  arrives  at  Equation  (30)  of  section  2.2.2.  The  discussion 
on  the  solution  of  this  equation  was  aiso  given. 


nGUHEC-9.  OBUQUE  SHOCK  GEOMETRY  FOR  A  CONE 
(SPHERICAL  COORDINATESJ 
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Frozen  Flow 

For  frozen  flow,  the  gas  is  allowed  to  change  chemical  composition  across  the 
shock  im,tantaneously,  and  tiie  values  of  y  and  Z  are  no  longer  constant  across  the 
shock.  Also,  since  the  flow  angle  behind  the  shock  is  not  constant  and  equal  to  the 
cone  angle,  as  was  the  case  for  a  wedge.  Equation  (C-17)  is  not  the  most  appropriate 
equation  to  use  in  the  soiction  process  since  this  involves  three  unknowns.  Using 
Eq^tiaas  (C'lS)  ?md  (C-14)  instead  requires  only  a  double  iteration  to  solve 
Equations  iC-18)  ttuoiigh  \C-20)  frt)m  the  body  to  the  shock.  Tlie  solution  proceeds  as 
follows: 

1.  Guess  values  of  o  and  pi/p2.  Perfect  gas  values  are  adequate. 

2.  Solve  Equations  (C-13)  and  (C-14)  for  p2,  h2. 

3.  Solve  the  appropriate  form  of  Equation  (73)  for  a  new  value  of  p2  [i.e.,  p2  = 
P2(P2.h2)orp2  =  p2  (P2»S)ifS  is  given  firom  blunt  b'-iy  flow!. 

4.  Repeat  steps  2  and  3  until  succesmve  changes  in  p2  within  some 
acceptable  tolerance.  11110  defines  the  values  of  p2,  P2,  h2  behind  the  shock 
of  guessed  angle  a. 

5.  Calculate  values  of  “a"  and  T  behind  the  shock  from  the  appropriate 
thermodynamic  curve  fit  equations  and  a  value  of  Z  =  Zfandy  =  yf  from 


and 


P2 

p,RT3 


The  values  of  y^  and  Zf  are  fixed  at  these  values  for  the  solution  cf  Equations 
(C-18)  through  (C-20). 

(Recall  from  the  wedge  Sow  computations  that  the  values  cf  y  and  y  are 
different.  The  question  arises  as  to  which  to  use  for  frozen  flow.  The 
answer  is  that  ei^er  can  be  used  or,  for  that  matter,  other  options  for  frozen 
flow  could  be  used.  However,  to  be  consistent  with  the  use  of  total 
temperature,  pressure,  and  density  relationships  along  the  body  where  the 
flow  is  in  fact  isentropic  and  the  isentropic  exponent  y  is  used,  the  isentropic 
exponent  y  is  the  one  that  will  be  used  for  frozen  flow  computations.) 

6.  Define  initial  values  cf  Vf  and  Vg  behind  Ihr  shock  frt>m  the  momentum 
and  energy  equations.  That  is 
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(V^)2  =  VjCoso 

1 

<V<,>2=-[vL-2'-s-<v4' 

7.  Numerically  int^^te  Equations  (C-18)  through  (C-20)  from  the  shock  to 
the  body  where  6  =  Oc  using  Fourth-Order  Runge  Kutta.  During  the  in¬ 
tegration,  new  values  of  p,  Vr,  Ve  are  computed  at  each  A6  interval.  At 
each  interval,  new  values  of  a^  and  p  must  be  computed  based  on  the  values 
ofp,Ve,Vr.  These  new  values  are 


(C-25) 

(C-26) 


2 


a 


(C-27) 


p  = 


(C-28) 


8.  At  the  body,  check  the  value  of  Ve.  If  it  is  not  zero,  guess  a  new  value  of  o. 

9.  Repeat  steps  1  throi^h  8  until  Ve  is  zero  to  some  acceptable  tolerance. 
Once  this  has  occurred,  the  entire  flow  field  between  the  shock  and  body  is 
defined.  Only  the  values  of  properties  at  the  body  surface  are  saved  as  they 
are  all  that  will  be  used  in  the  Shock  Expansion  Theory. 

10.  Calculate  the  entropy  behind  the  shock  from  the  appropriate  thermofit 
equation.  This  value  of  entropy  will  be  used  later  in  the  shock  expansion 
theory.  If  entropy  values  are  known  firom  the  blunt  body  stagnation 
solution,  this  step  can  be  omitted. 


Equilibrium  Chemically  Reacting  Flow 

For  equilibrium  (hemically  reacting  flows,  the  ^imputation  process  illustrated 
for  frozen  flows  is  used  with  two  exceptions: 

1.  The  values  y  and  Z  are  allowed  to  vary  between  the  body  and  shock. 

2-  The  valves  of  “a"  and  “p”  at  each  point  in  the  computation  process  are 
computed  via  thermofit  Equation  (73)  based  on  the  latest  computed  value  of 
*‘p”  and  the  value  of  which  is  constant  behind  the  shock.  That  is 


a  =  a  (p,  S) 
P  =  p{p.S) 


(C-29) 


C-15 
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Both  of  these  relations  are  given  in  Reference  C-6  and  they  replace  the  equations  for  p 
and  a  in  the  frozen  flow  computations  discussed  previously. 

Knowing  the  values  of  Ve,  Vr,  S,  p,  a,  and  p  during  the  solution  process  of 
Equations  (C-18)  through  (C-20),  the  remaining  thermodynamic  properties  can  be 
computed  as  follows: 


> 

1 

« 

> 

1 

II 

(C-30) 

e  =  h  —  p/p 

(C-31) 

[vjtvjp 

M=  - - - 

a 

(C-32) 

T  =  T(p.S) 

(C-33) 

Z=  ** 
pRT 

(C-34) 

a"p 

y  =  h/e  ;  y  =  - 

(C-35) 

p 


To  illustrate  the  effects  of  real  gas  computations  on  conical  shapes.  Figures  C-10 
throu^  C-13  have  been  prepared.  These  figures  duplicate  the  results  of  Reference 
C-5  within  numerical  accuraqr.  Results  here  are  similar  to  those  of  the  wedge  caty 
presented  previously  except  the  frozen  and  chemical  equilibrium  cases  are  not  the 
same.  Hie  trends  of  decreased  shock  angle,  higher  tenqieratures,  lower  densities, 
and  near  constant  pressures  of  real  gases  compared  to  perfect  gases  holds  true  for  the 
cone  also.  These  trends  will  not  be  true  for  PME  flows  considered  next  in  Appendix  C. 

Before  proceeding  to  discuss  PME  flows,  discussion  of  a  problem  that  occurred 
in  the  solution  of  equililn’ium  chemically  reacting  flows  over  cones  using  the  curve 
fits  of  Reference  C-€  will  be  giv^.  Since  entropy  is  constant  between  the  shock  and 
body  for  a  wedge  or  cone  and  is  constant  along  streamlines  for  equilibrium  flows,  it 
seemed  natural  to  use  tills  as  one  of  the  independent  variables  [see  Equation  (C-29)l. 
However,  when  attempting  to  oimpute  the  entiopy  firom  internal  energy  and  density, 
it  was  foui^  that  slight  discontinuities  (or  jumps)  occurred  in  the  entropy  values 
between  various  segments  of  the  curve  fits  in  Reference  C-6.  These  jumps  and 
attempts  at  smoothing  the  solution  are  illustrated  in  Figure  C-14.  Unfortunately, 
when  both  the  actual  and  smoothed  values  of  entropy  were  used  in  the  solution  of  the 
equations  for  a  cone  \  (C-18),  {C-19),  (C-20),  and  {C-29)  1,  discontinuities  occurred  in 
the  density  and  temperature  at  the  cone  sur&ce.  An  example  of  the  cone  siufacs 
density  is  given  in  Figure  C-15  (curve  labeled  p  =  p  (p,  S).  Note  that  the  smoothed 
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CONE  HALF-ANGLE  (•) 


FIGURE  C-12.  EFFECT  OF  FREESTREAM  VELOCITY  AND  CONE  ANGLE 
ON  CONE  SURFACE  TEMPERATURE  (H  =  100k  R) 


CONE  HALF-ANGLE  (•) 


FIGURE  C-13.  EFFECT  OF  FREESTREAM  \^LOCITY  AND  CONE  ANGLE 
ON  CONE  SURFACE  DENSITY  (H  =  100k  R) 
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values  show  some  improvement  over  the  unsmoothed  entropy  values  but  both  results 
deviate  too  much  from  the  solution  of  Reference  C-5  (up  to  10  percent). 

To  get  around  this  problem,  several  other  alternatives  v/ere  investigated.  The  first  of 
these  was  to  solve  the  conical  equations  using  the  total  enthalpy  and  locally 
computed  values  of  p,  Vr,  and  Ve  from  Equations  (C-18)  to  (C-20)  to  solve  for 
intermediate  values  of  enthalpy  from  Equation  (C-30).  Internal  energy  could  be 
computed  using  Equation  C-31  and  then  density  and  speed  of  sound  could  be  updated 
with  p  (h,  p)  and  a  (e,  p).  Using  this  approach,  the  results  labeled  p  (h,  p)  of  Figure  C- 
15  were  obtained.  These  results  agree  almost  identically  to  those  of  Hudgins.C-S 
However,  they  still  do  not  use  one  of  the  key  parameters  that  will  be  needed 
downstream  in  the  computation  process:  the  constancy  of  entropy  along  the  body 
surface. 


Q  t  1  l - 1 - 1 - 1 - 1  I  * 
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nCURE  C-14.  ENTROPY  VERSUS  ENERGY  FOR  CONSTANT 
VALUES  OF  DEXSTTY 

A  final  attempt  was  made  to  rise  entropy.  The  approach  was  to  solve  the  conical 
flow  equations  (C-18)  through  (C-20)  with  conditions  behind  the  shock  determined  as 
previously  discussed. 

Behind  the  shock,  internal  energy  was  found  from  Equation  (C-31).  Entropy 
was  found  from  Equation  73  where  S  =  S  (e,  p).  Using  this  value  of  S  and  p,  a  new 
value  of  e  was  computed  from  e  =  e  (p,  S).  This  value  of  e  was  compared  to  the 
previous  value  of  e.  If  the  two  did  not  agree  to  some  desired  accuracy,  a  new  value  of 
S  was  guessed  and  an  iteration  took  place  until  the  two  values  of  e  agreed.  This  then 
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was  the  value  of  S  that  was  held  constant  in  the  wedge  or  cone  flows  and  along  the 
body  surface  (S  =  Si).  Using  this  approach  to  determine  Si,  Equations  (C-18)  through 
(C-20)  were  solved  in  conjunction  with  (C-29).  These  results  are  shown  in  Figure  C- 
15  as  p  (p.  Si)  where  Si  refers  to  the  value  of  entropy  found  from  iteration,  lliese 
results  are  very  close  to  those  where  density  is  computed  from  p  =  p  (h,  p)  up  to  cone 
half-angles  of  20°.  The  largest  error  up  to  cone  half-angles  of  30°  is  about  4  percent 
and  for  cone  angles  less  than  20°  is  less  than  1  percent.  This  then  is  the  approach 
chosen  to  compute  entropy  that  will  be  used  throughout  the  flow  field  computations. 
Figures  C-16  through  C-18  show  the  comparison  of  the  various  cone  solutions  on 
surface  pressure,  temperature,  and  density  for  two  velocities.  Note  that  if  one  were 
interested  in  pressure  only,  any  of  the  approaches  are  adequate.  However,  to  get  the 
most  correct  value  of  temperature  that  is  needed  for  aerothermal  analysis,  the 
approach  of  iterating  to  find  entropy  appears  to  give  results  closer  to  the  correct 
solution. 

The  only  explanation  that  can  be  offered  for  this  discrepancy  is  the  fact  that  the 
curve  fits  are  stated  to  be  accurate  to  within  2  percent.  The  curve  fits  are  formulated 
basically  independent  of  each  other.  On  the  other  hand,  by  using  results  from  other 
curve  fits  to  force  a  more  correct  value  of  entropy  to  be  computed  requires  more 
consistency  in  the  conservation  laws  and  thermodynamic  properties.  The  approach  of 
iterating  to  find  entropy  takes  a  little  more  computer  time.  However,  for  the  work 
herein  where  entropy  will  be  held  constant  along  the  body  surface  and  hence  only 
computed  once,  the  computational  time  is  insignificant. 


DENSITY 

(SLUGSm3*103) 


FIGURE  C-15.  COMPARISON  OF  DENSITY  ATTHE  CONE  SURFACE  USING  VARIOUS 
WAYS  OF  UPDATING  DENSITY  IN  FLOW  FIELD  SOLUTION  {h  =  100  k  ft) 
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FIGURE  C-16.  CONE  SURFACE  PRESSURE  COMPUTED  WITH  VARIOUS  WAYS 
OF  IMPLEMENTING  REAL  GAS  EFFECTS  (h  =  100k  R) 
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FIGURE  C-17.  CONE  SURFACE  TEMPERATURE  COMPUTED  WITH  VARIOUS  WAYS 
OF  IMPLEMENTING  REAL  GAS  EFFECTS  (h  =  100k  R) 
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nGUBEC-18.  CONESUKFACEDENaTYCOMPUTEDWITHVABIOUSWAYS 
OF  QfPUSMENTtNG  REAL  GAS  EFFECTS  (h  =  100k  ft) 
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APPENDIX  D 

COMPUTATION  OF  PROPERTIES  ACROSS  EXPANSION  WAVES 
IN  FROZEN  AND  EQUIUBRIUM  CHEMICALLY 
REACTING  FLOWS 
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COMPUTATION  OF  PROPERTIES  ACROSS  EXPANSION  WAVES 
IN  FROZEN  AND  EQUIUBRIUM  CHEMICALLY 
REACTING  FLOWS 


Prandtl-Meyer  Expansion  (PME)  flow  is  standard  in  the  aerodynamics 
literature  (see,  for  example.  References  D-1  thrcugh  D-5).  For  convenience,  a  brief 
summary  of  the  theory  and  results  is  presented  in  this  appendix. 

The  PME  consists  of  an  isentropic  turning  of  a  streamline  from  its  initial 
conditions  at  state  1  to  its  final  conditions  at  state  2  (see  Flgiure  D-1).  Note  that, 
while  this  is  an  isentropic  process  (for  equilibrium  air),  this  process  applies  from 
streamline  to  streamline,  and  each  streamline  can  have  different  conditions. 
Referring  to  Figure  D-1,  this  means  that  we  do  not  necessarily  have  constant 
conditions  in  the  entire  flow  field  leading  up  to  the  turn.  Indeed,  in  general,  they  are 
not.  However,  on  the  body  surface,  since  it  is  considered  a  streamline,  the  flow  is 
isentropic  so  long  as  nonequilibrium  conditions  do  not  occur.  The  Mach  waves  are 
actually  characteristic  lines  along  which  no  flow  information  is  transmitted 
upstream.  Of  course,  as  conditions  from  streamline  to  streamline  change,  the  Mach 
lines  or  characteristics  will  curve  to  accommodate  these  changes.  If  conditions  ahead 
of  the  turn  are  constant  throughout  the  flow  field,  the  characteristics  will  be  straight. 
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Another  important  feature  is  that  the  velocity  parallel  or  iangentiJl  to  the 
Mach  line  will  be  constant  on  either  side.  This  comes  hrom  the  conservation  of  mass 
and  momentum  equations  used  in  derivation  of  the  equations  for  flow  across  shock 
waves.  This  result  is  independent  of  type  of  flow  (i.e.,  shock  or  expansion  wave). 
Looking  at  an  incremental  change  in  the  velocity  as  the  flow  turns  around  a  comer, 
we  can  derive  the  equations  for  P.M.  flows.  Referring  to  the  bottom  ^etch  of  Figiare 
D-2,  we  can  write  the  foUowii^  from  geometry: 


dV^cosp  =  (V  +dV)(-d  0)  =  -Vd9 

(D-1) 

dV 

sinu  = 

d  V 

(D-2) 

n 


The  figure  shows  but  the  process  applies  to  any  portion  of  the  expansion,  so  the 
subscripts  have  been  dropped  in  the  analysis.  Substituting  (D-2)  into  (D-1),  we  get 


dV 


cosp 


sin|i 


-Vd9 


or 


dV 

-d0  =  ootp  — 


(D-3) 


But,  since 


then 


anp  = 


1_ 

M’ 


cotp 


-  1 


Equation  (D-3)  becomes 

-de=  1  ^ 

From  the  adiabatic  energy  equation,  one  can  write 

H  =  Const  =  h  +  - 
2 

Differentiating  (D-5),  solving  for  dV/V,  wfc  obtain 


(D-4) 


(D-5) 
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Also 


V  _  V2>:H-fa) 
a  a 


(D-T) 


Substituting  Equations  (D-€)  and  (D-7)  into  (D-4),  the  general  equation  for  PME 
flows  is  obtained 


-de  = 


[2{H-h> 


2(H-h) 


(D^) 


Tliis  equation  is  quite  general  as  the  only  assumption  is  iseniropic  flow  along  a 
streamline.  It,  therefore,  could  sqmly  to  an  isentropic  expansion  or  compression  and 
to  perfect,  frozen,  or  equilibrium  chemically  reacting  flows. 
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PESFECTGASFLOW 

Ifweaarame  a  perfisct  gas,  then  the  q>8cific  heats  are  cGnstant(Cp,Cv)  and  Z  = 
1.  Thu8,h  =  CpTacda^  =  yRT.  With  these  8ubstitutiocs,£qaat?cnCD>8)  becomes 


-de=  - 


2(T, 


("t  “  ^ r** 

,-T)i(y-l)\T  /  1 


fD-9) 


But  since 


T  2 


Equation  (D-9)  can  be  rewritten  in  the  form 

VyT^l  dM 


-de  = 


y  —  1  •  M 

1  +  - - M* 


(D-IO) 


In  int^ral  form,  th:  omes 

f  ft 


“2 


M  y  ~  t 
“1  1 +  - — 


L-  i!! 

M 


(D-ll) 


If  6  —  is  deSned  as  the  flow  angle  that  aaresiMnds  to  M  =  1,  Equatimi  (D-ll)  can 

be  rewritten  as 

rIL  rW  rM.  rl 

-  de-  de=  I  ?KM)6M+  Fiy.M)dM  {D-i2) 

Je* 


where 


?(y.M)  = 


1  +  M*  M 
2 
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The  integrais  of  Equation  (D-12)  can  be  migrated  by  the  method  of  partiai  &actions 
to  obtain: 


v{M)*  f  FG«.y)dM=  Ip  (D-13) 

ii  /-I  ty  +  1  J  ;  J 

-Substituting  into  Equation  (D-12),  there  is  obtained 

-02  +  0,=Ae= -HvOljJ-vOli)  {D-i4) 

AO  and  are  known  quantities  so  v(M,)  can  be  computed  firom  Equation  <D-13}. 
▼(M2)  is  then  computed  from  Equation  (D-14),  from  which  M2  is  solved  fur  in  a 
numerical  itoratian  process  by  Equation  (D-13).  Knowing  M2  and  the  fact  that  the 
flow  along  the  streamline  is  isentropic,  all  the  other  flow  field  properties  can  be 
computed  by  the  tcta!  pressure,  denmty.  snd  temperature  relationships  which  are 
constant.  That  is 


a?-i5) 

(D-16) 

P2=P2/RT2 

(D-17) 

e  adiabatic  'energy  relations  we  have 

(D-18) 

*2=S'^2 

(D-19) 

1^2 

V  ^ 

(D-20) 

H  =  C  T, 

(D-21) 

9  0 

(D-22) 
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FROZEN  FLOW 

If  the  flow  h  frozen,  the  condiUons  behind  the  shock  for  a  ncimal  shock  on  a 
blunt  body  will  he  taken.  These  conditions  come  from  the  normal  shock  real  gas 
commutations. 

For  wedge  cr  cone  flows  (sharp  or  pointed  bodies),  the  conditions  behind  an 
oblique  shock  can  be  used.  Again  these  results  come  from  real  gas  computations  for 
oblique  shock  wav&'.  :a  Appendix  B.  Hence,  the  initial  conditions  for  frozen  flow  on 
blunt  or  pointed  bodies  or  blunt  or  sharp  win^  are  known.  This  gives  values  of 


which  are  prx^iertiesofihe  gas  and  its  chemical  states. 

The  rest  of  the  flow  can  then  be  solved  similar  to  a  perfect  gas.  That  is, 
Equaiiocs  (D-13)  through  {D-16>  are  solved  where  fr(»en  values  of  yr  are  used  for  y. 

Equations  (D-17)  through  {l)-22)  are  replaced  with 


P2  =  P2 


»2=-' 


“2 


'^2=^3  3, 


^2  =  \  -  Pi'Pz 


EQUIUBPJUM  CHEMICALLY  REACTING  FLOW 

For  the  case  of  chemically  reacting  flows  where  the  ratio  cf  specifle  heats  and 
compressibility  factor  can  change  as  the  flow  expands  aroimd  the  surface,  one  can  no 
longer  assume  a  perfect  gas  and  int^rate  Equation  (D-d)  in  dosed  form.  Equatitm 
(D^)  must  be  solved  through  numerical  integration.  In  differential  equation  form. 
Equation  (Et^}  becomes 
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d(K-h) 

de 


+ 


2(H-h> 

f2(H-h) 


(D-23) 


Once  again,  Fnurth-order  Runge  Kutta  or  other  appropriate  numerical  integration 
schemes  can  be  used.  The  A6  interval  (O2-81)  is  divided  into  a  given  number  of 
equaliy-^>aced  intervals.  The  initial  values  of  H,  h,  61, 82  (defined j,  a,  M,  p,  p,  S,  etc. 
are  known.  These  all  come  from  the  conditions  behind  a  normal  shock  or  oblique 
^ock.  Of  c(*iirse,  the  conditions  behind  the  shock  must  be  taken  to  the  surface  of  the 
body  and,  in  the  case  of  a  blunt  body  or  leading  edge,  taken  to  the  point  where  the 
PME  process  begins.  For  a  blunt  body,  this  process  is  discussed  in  the  Mixiifisd 
Newtonian  Theory  section  of  the  test  and  in  the  discussion  of  wedges  and  cones  is 
Aj^ndix  B. 

One  positive  aq>ect  about  using  the  Fomrth-order  Runge  Kutta  integration 
routine  (as  opposed  to  Simpson’s  rule  or  some  other  integration  process  where  h2  is 
required)  is  the  fact  that  no  initial  guess  of  h2  is  necessary  nor  any  follow-on  iteration 
to  correct  Uiis  guess  based  on  the  boundary  conditions.  Unfortunately,  while 
Equation  (D-23)  is  given  ia  a  very  straightforward  manner  for  integration,  a(h,  S)  is 
not -available  from  either  Reference  I>-6  or  D-7.  Informal  communication  with 
Tannehill  also  revealed  that  this  particular'  thermofit  equation  still  had  not  been 
developed.  As  a  result,  either  a  resort  to  tablesD-S  must  be  used  for  a{h,  S)  or 
E-quation  (D-23)  must  be  put  in  a  form  compatible  with  the  thermofit  equations  of 
References  D-fi  and  D-7.  Since  the  thermofit  equation  process  is  much  faster  Uian 
table  look-up  procedures,  and  at  only  a  small  sacrifice  in  accuracy,  two  different 
alternatives  were  investigated  for  solviner  Equation  (D-23}. 

The  first  approach  uses  the  Prardtl-Meyer  Equation  D-4;  that  is 


^  _  -V 

d8  y/‘y2_  j 


24) 


From  the  one-dimensional  momentum  equation 


dV  =  - 


dp 

pV 


(B-25) 


so  that  substituting  Equation  (D-25)  into  (D-24}  one  obtains 


D-8 
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dp  pV^ 


£ 


(D-26) 


To  solve  Equation  CD-26),  jaie  knc'Jtrs  initiai  values  of  p,  8,  A0,  p,  a,  and  V. 
Equation  (D-26)  is  integraced  using  again  Fourth-order  Rimge  Kutia  to  obtain  a 
value  of  p  at  the  next  integration  interval.  Knowing  p  and  the  fad  that  entropy  is 
constant,  updated  values  of 

p  =  p(p.S) 

a  =  a(p,S}  ^  (D-2f) 

e  =  e(p.S)  J 

can  be  computed  firom  the  thermoflt  equations  of  Reference  D-7.  Then 
h  =  e  4-  p/p 

and  (D-28) 

V  =  V  2{H-h) 

(D-27)  and  (D-28}  are  then  used  along  with  (D-26)  in  the  integration. 

A  second  method  investigated  involved  the  simultaneous  solution  of  E<{uations 
(D-24)  and  (D-28).  That  is,  knowing  initial  values  of  V  and  a.  Equation  (D-24)  was 
solved  to  obtain  an  updated  value  of  V.  This  new  value  of  V  plus  initial  values  of  p,  p, 
and  a  were  used  to  solve  for  p  at  the  next  integration  step.  Using  p  and  S,  p  and  a 
were  updated  usii^  Equation  {D-27).  Equations  (D-24)  and  {D-26)  were  then  solved 
at  the  next  integration  step.  Figures  D-3  and  D-4  give  the  results  of  these  two 
approaches  for  pressure  and  temperature  (referred  to  on  the  figures  as  methods  1  and 
2)  and  compares  these  results  to  Heimsl^-s  for  equilibrium  ilcw.  It  should  be  pointed 
out  that  Heims^  used  table  lock-up  procedures  for  the  real  gas  properties  so  his 
results  should  be  close  to  exact  Note  that  the  method  1  approach  of  only  integrating 
one  difierentia!  equation  ^ves  more  accurate  results  compared  to  Heims. Hence, 
this  is  the  af^roach  used  herein  for  integrating  the  PME  for  equilibrium  flew. 

Figures  D-5  through  D-8  give  pressure,  temperature,  density,  and  Mach 
number  fer  perfect,  equilibrium,  and  frozen  flow  cases.  The  initial  conditions  for  the 
expansion  were  the  same  as  those  on  Figures  D-3  and  D-4;  that  is,  Ti  =  6140",  M:  = 
1,  and  Pi  =  1.2  atm.  It  is  quite  interesung  to  note  that  while  pressure  is  not  affected 
to  any  appreciable  extent  in  a  compression  process  compared  to  perfect  gas 
computations,  it  increases  m  an  expansion  process  of  equilibrium  versus  perfect  gas 
computations  (see  Figaro  D-5).  Also,  the  equilibrium  temperature  increases  in  an 
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FIGURE  D-4.  TEMPERATURE  COMPAHLSON  OF  INTEGRATION  TECHNIQUES 
FORPME  OF  EQUILIBRIUM  FLOW  (Ml  =  i.O.Ti  =  6i4(rK.  pi  -  i.2aUn) 


FIGURE  ns.  PERFECT.  FROZEN,  SQUILffiRIUM  PME;  TEMPERATURE  \^RSUS 
FLOWDEFLECilON  ANGLE  iMi  =  l.O.Ti  =  614(nLpi  =  l-2atin) 


FIGURED-?.  PERFECT,  FROZEN,  EQUILIBRIUM  PME:  DENSITY  VERSUS 
FLOW  DEFLECTION  ANGLE  (Mj  =  1.0,  Ti  =  6140'K,Pi  =  1.2  atm) 
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FIGUHED-S.  PERFECT,  FROZEN,  EQUILIBRIUM  PME:  MACH  NUMBER  VERSUS 
FLOW  DEFLECTION  ANGLE  (Mi  =  1.0,  Ti  =  6140*K,pi  =  1.2  aUn) 
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expansion  as  opposed  to  decreasing  in  a  compression  compared  to  perfect  gas 
computations  (see  Figure  D-6).  Note  that  density  of  equilibrium  Sows  versus  perfect 
gas  is  nearly  constant  for  an  expansion  (Figure  D-7),  whereas  it  increases  for  an 
equilibrium  compression  compared  to  a  perfect  gas.  Finally,  Mach  number  decreaees 
for  an  equilibrium  expansion  compared  to  a  perfect  gas  (Figure  D-8).  Frozen  Qow 
cases  are  also  shown  on  Figures  D-5  through  D-8  where  the  gas  conditions  that  exist 
at  the  initial  condition  state  are  held  constant  throughout  the  expansion  process. 
Two  values  of  the  initial  state  are  held  fixed  corresponding  to  the  isentropic  exponent 
(  y  =  a2p  /  p)  and  the  so-called  “equivalent”  gamma  (  y  =  h/e).  Note  that  for  the 
expansion  process  there  is  very  little  difference  between  the  two,  which  indicates  that 
for  this  set  of  initial  conditions  both  values  of  gamma  are  about  the  same.  This  may 
not  necessarily  be  the  case,  however,  had  the  fi^ozen  flow  values  of  y  and  y  been  these 
that  existed  behind  a  normal  shock  on  a  blunt  body. 
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APPENDIXE 

DERIVATION  OF  CHARACTERISTie;  COBfPATIBlLlTT 
CONDITIONS  IN  STREAMLINE  COORDINATES 


E-l/'E-2 
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DERIVATION  OF  CHARACTERISTIC  COMPAUBH  JTY 
CONDITIONS  IN  STREAMLINE  COORDINATES 


Our  equations  of  motion  in  vector  form  are  as  follows: 
Continuity 


V  (p  V)  =  0 


Momentunt 


P 


D  V 
Dt 


=  -Vp 


or  for  steady  flow 


p  V  -  V  V  =  Vp 


Energy 


—  -r  k  =  H.  =  conscaot 
9.  t 


State 


(E-l) 


{E-2) 


(E-3) 


h  =  h(p,p} 
a  =  a(p,p} 


(E-4) 


would  like  to  6rst  transform  these  equations  to  the  intrinsic  or  streamline 
coordinate  system.  In  streamline  coordinates,  s  is  along  the  body,  which  is  a 
streamline,  and  n  is  normal  to  it.  Also 


V 


(E-5) 


Substituting  (E-5)  into  (E-1)  and  (£-2)  there  is  obtained 
(pV)  V  -  V(pV)  =  0 


E-3 


(E-6) 
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add 


-  v)v  =  -  ^ 

S  t  t  S  p 


(E-7) 


Take  the  dot  product  of  Equation  (E-7)  wilhe^  to  obtain 
-  |{^  -  V)  e  I  -  Vp 

n  [  •  a\  p  n 

Using  the  vector  identities 


e  •v)e  =  vxc  xe 

•  I  a  • 


-  (ori^  I  )  =  01x1^  -  T.^  ) 

n  sa  s  ^  o 

andthe&cttfaat^x^  ='^,  Equation  E-8  become 


(E-8) 


^  <  V  X  ^  )  +  — —  ^  -  ?  p 

‘  '  pV2  “ 


(E-9) 


Also,  for  equilibrium  Cow,  entropy  is  conserved  along  streamlines,  which 
requires 


•  V I 


(E-10) 


Substituting  Equation  (£-10)  into  the  continuity  Equation  (6)  and  expanding  there  is 
obtained 


t'S  -  V)V  -  -  -  Vp  -  V(V 

•  pa 


(E-11) 


But  Com  the  momentum  Equation  (E-7),  after  taking  the  dot  product  with  one 
obtains 


.  vv  =  -  ~  -  Vp)  (E.12) 

Substituting  (E-12)  into  (E-11/  and  performing  the  algebra,  one  obtains 
-2 

-  ^-Vp  +  V-  e  — 0  (E-13) 

PV^  • 


E-4 
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where 


-  1 


The  next  step  in  the  transformation  of  the  equations  of  motion  (E-l)  and  (£-2) 
into  streamline  coordinates  is  to  relate  the  iinit  vectors  e,  and  'en  to  those  in  the 
cylindrical  coordinates  and  To  do  this,  simpfy  rotate  the  e^,  ^  axis  through  an 
angle  0  to  the  e,,  e^  plane.  This  gives 


^  =  oos9'^  +  sinO'e 

•  z  r 


(E-14) 


^  =(  — 8in0)^  +  0086'^ 

o  z  r 


(E-15) 


The  vector  operators  in  OTiindrical  coordinates  are 


ax  *  ar 


(E-16) 


^  1  f  a(rA  ) 

V  -  A  =  - - 

r  I  dx 


a(rA^) 


(E-17) 


Vx  A  =  - 
r 


i 


For  axisymmetric  flow 


so  that 


—  =-  A-  =  0 
8$  ^ 


Vx  A  I  — 
v\  dx 


as-18) 


E-5 
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Nowingenei-alx  =  x(s,n>andr  =  r(s,n).  Thus,  using  the  chain  rule,  one  can 
write 


f  ±] 

=  -i 

''  ds  \ 

a  /  dn\ 

dn  \  dx  / 

\  dxj 

ds  • 

‘  dx/ 

\  dr  / 

II 

as\ 

arl 

+ 

d  / 

dn  V  dr  / 

Using  £qiiations  (£-14)  and  (£-15),  Equations  (£-19)  and  (£-20)  become 


—  =  cosS —  —  sm6  — 
dx  ds  dn 


a  a  a 

—  =  sinG  —  +  oosG  — 
dr  ds  dn 


(E-19) 


(£-20) 


(E-21) 

(£-22) 


Now  substituting  Equations  (£i-14),  (£-16),  and  (£-17)  into  the  continuity  equation 
(E-13),  one  obtains 


il 

PV^ 


[  [(cosG)^^  +  11  p 


la  „  Id 
+  —  —  (roosG)  +  —  —  (rsinG)  =  0 
r  dx  r  dr 


(E-23) 


Now  substituting  Equations  (E-21)  and  (£-22)  into  (£-23),  performing  the  algebra, 
there  is  (ditained 


dp  dG  _  — sinG 

pY^  dn  r  (£24) 

This  then  is  the  continuity  equation  in  streamline  coordinates  in  terms  of  the 
dependent  variables  p  and  6. 

For  the  momentum  equation  (£9),  where 


NAVSWCTRSO^ 


substitute  Equations  (E>18),  {£-15),  (E-16)  and  perform  the  dot  product  to  obtain 


3  3  1  ^  3o  dp  I  ^ 

—  lm0)  —  —  (0080)  +  — -  j  — Bm0  —  +  OOS0  —  1=0  (E-25) 

ax  ax  pV^  l  3x  dr  I 


Finally,  substitute  the  relations  (E-21)  and  (£-22)  into  (E-25)  and  perform  the 
algebra  to  get  the  final  form  of  the  n  momentum  equation 


1  dp  36 

— 0 
P  y2  dn  3s 


(E-26) 


Equations  (£-24),  (E-26),  and  (E-10)  are  the  continuity  and  momentum  equations  in 
streamline  coordinates.  Hey,  in  conjunction  with  (E-d)  and  (£-4)  allow  a  complete 
solution  of  the  flow  field  for  ardsymmetric  bodies  at  zero  angle  of  attack  in  steady, 
invisdd,  adiabatic  flow  of  an  equilibrium  chemically  reacting  fluid. 

The  next  step  in  the  derivation  of  the  charactmistic  compatibility  relations  is  to 
put  the  equations  in  a  form  for  solution  of  the  partial  derivatives.  It  will  then  become 
envious  as  to  why  the  equations  were  put  in  the  form  of  streamline  coordinates  before 
the  cocmatibility  equations  were  derived.  In  essence,  this  allowed  reduction  of  one 
viriable  since  in  streamline  coordinates  the  velocity  is  along  the  streamline,  whereas 
in  rectangular  or  (ylindrical  coordinates  the  velocity  has  components  along  the 
req)ective  axes.  Also,  the  entropy  is  constant  along  the  streamline,  simplifymg  the 
energy  equation. 

Repeating  the  equatiems  for  our  mathematical  model  for  convenience,  we  have 


^  dp  30  _  — sinO 

pV^  3s  3n  r  (E-24) 


1  3?  30 

— -  —  -H  —  =  0 
py2  dn  3s 


Also,  we  can  write 


dp  dp 

dp  =  —  ds  +  —  da 
ds  3n 


(E-26) 


(E-27) 


30  .  36 

d0  =  —  ds  —  dn 

3s  dc 


(£-28) 


Putting  Equations  (£-24),  (£-26),  and  (£-28)  in  symbolic  determinant  form  one  has 


E-7 
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dp 

dp 

dO 

de  ; 

ds 

dn 

ds 

dn 

f 

0 

0 

1 

— sin6 

pV^ 

r 

0 

1 

1 

0 

j 

0 

pV^ 

ds 

d  n 

0 

0 

dp 

0 

0 

ds 

d  n 

de 

One  can  then  solve  for  any  one  of  the  partial  derivatives  through  a  simple  deter¬ 
minant  process.  For  example 


au  _  N 
ds  D 


where  N  = 


0 


— san6 


r 


1 


0 


0  0 


d  s  d  n 


dp  0 


0 


dO 


(E-29> 
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and 


D  = 


0  0  i 

PV^ 


0 


1 


0 


ds  dn  0  0 

G  0  ds  dn 


(E-30) 


The  ci^rscteristk  curves  are  de0s^  as  those  Hnes  along  vdiich  the  derivatives  of 
fluid  properties  mi^t  be  discontinuous.  To  determine  these  curves  requires  the 
denominator,  Equation  (E-30),  to  go  to  zero.  Performing  this  (iteration  and  carrying 
outthe  algrina  one  obtains  the  result 


dn  _  ^ i 
ds  ■  "8 


±tanp 


(E-31) 


l^is  says  that  along  lines,  called  characteristic  lines  (which  are  also  equal  to  Mach 
lin»  because  the  lines  make  the  angle  p  with  respect  to  the  streamline),  the 
derivative  Sow  jn'tqiertie  may  be  discontinuous. 

Now  in  order  for  the  derivatives  to  exist,  if  the  denominator  is  zero,  then  the 
numerator  mmt  also  be  equal  tc  zero.  Setting  Equation  (E-29}  to  zero  and  can3nng 
out  the  math,  rme  obtains  the  soKralled  compatibility  condition 


d9  ds  +  ds 


+ 


da  ds 


=  0 


(E-32) 


Substituting  (£-31)  into  (E-32),  realizing  that  the  projections  (Stto  the  s  axes  of  the 
Cl  and  C2  coordinates  are 

d 5  =  dC,  CHS  -i-  dC.cosa 
1  -  2  - 


then  one  obtains 


S-S 
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iz-  + 

de 

—  ^Osnip 

pV* 

dC, 

dc,  - 

r 

(E-33) 

p 

dp 

de 

—  smOmp 

pV* 

dC, 

dc^- 

r 

(E-34) 

Here  Ci  is  called  the  left  running  and  C2  the  right  running  characteristic 
corresponding  to  the  plus  and  minus  signs  of  Equation  (E-31),  respectively. 

Note  that  along  the  characteristic  lines,  the  partial  differential  equations  (E-1) 
and  (E-2)  have  been  reduced  to  two  first-order  ordinary  differential  equations.  Thisis 
the  major  advantage  of  the  method  of  characteristics  as  it  always  reduces  the 
equations  of  motion  by  one  dimensicm  whmi  solved  along  the  characteristic  lines  or 
surfiices. 
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